AD-  78  4  1  91 


DEVELOPMENT  AND  EXPERIMENTAL  VERI¬ 
FICATION  OF  PROCEDURES  TO  DETERMINE 
NONLINEAR  LOAD- DEFLECTION  CHARAC¬ 
TERISTICS  OF  HELICOPTER  SUBSTRUCTURES 
SUBJECTED  TO  CRASH  FORCES.  VOLUME  I. 
DEVELOPMENT  OF  SIMPLIFIED  ANALYTICAL 
TECHNIQUES  TO  PREDICT  TYPICAL  HELI¬ 
COPTER  AIRFRAME  CRUSHING  CHARAC¬ 
TERISTICS  AND  THE  FORMULATION  OF 
DESIGN  PROCEDURES 

G.  Wittlin,  et  al 


Lock  heed- Califor  ni  a  Company 


Prep  ared  f  or: 

Army  Air  Mobility  Research  and 
Development  Laboratory 


May  1974 


DISTRIBUTED  BY: 


National  Technical  information  Service 
U.  S.  DEPARTMENT  OF  COMMERCE 

5285  Port  Royal  Road,  Springfield  Va.  22151. 


_ Unclassified _ 

SECURITY  CLARIFICATION  OF  THIS  PAPE  (Ww  Dm*  Imm« 


< 


f 


REPORT  DOCUMENTATION  PAGE 


1.  USF5X7  numSKR 
USAAMRD J-TR-74-12A 


REAS  mSTRUCnOMt 
BEFORE  COMPLETWO  FORM 


12.  OOVT  ACCESSION  MO. 


2.  RECIPIENT'S  CATALOOJIJJMAER 

ah 


I  CATALOG MUMS  Eft 

7$T/9/ 

EPORT  A  PERIOD  COVERED 


«.  TITLE  (m4  tmBHIlt) 

DEVELOPMENT  AND  EXPERIMENTAL  VERIFICATION  OF  PRO¬ 
CEDURES  TO  DETERMINE  NONLINEAR  LOAD-DEFLECTION 
CHARACTERISTICS  OF  HELICOPTER  SUBSTRUCTURES 
SUBJECTED  TO  CRASH  FORCES  -  VOLUME  I 


TYPE  or  REPORT  A  PERIOD  COVERED 

Final  report 

A.  PERFORMING  ORC.  REPORT  NUMBER 


7.  AUTHORS 

Go  Nittlin 
K.  C.  Park 


N.  CONTRACT  OR  GRANT  NUMAERfaJ 

Contract  DAAJ02-73-C-00S1 


A.  PERFORMING  ORGANIZATION  NAME  ANO  ADDRESS 

Lockheed-Califomia  Coarptny 
Burbank,  California  91503 


10.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  A  WORK  UNIT  NUMAERS 

Task  1F162205AH88 


It.  CONTROLLING  OFFICE  NAME  ANO  ADDRESS 

Eustis  Directorate 

U.  S.  Any  Air  Mobility  R$D  Laboratory 
Fort  Eustis,  Virginia  23604 

T4.  MONITORING  AGENCY  NAME  A  ADDRESSfH  ditto  not  horn  Controlling  Olllco) 


12.  report  date 

May  1974 


IS-  NUMAER  OF  PAGES 


IU.  SECURITY  CLASS.  {ot  thlo  ropo n) 

Unclassified 

IS*.  OEC^ASSIFICATION/OOWNORAOING 


IS.  ChSTRIAUTION  STATEMENT  (of  thlt  Boporl) 


Approved  for  public  release;  distribution  unlimited. 


IT.  OISTRIAUTION  STATEMENT  (ol  U i*  ototrocl  ontorod  In  Block  20,  It  dlttotdnl  Bom  ftaporf) 


IS.  SUPPLEMENTARY  NOTES 

Volune  I  of  a  two-volume  report 


>•  KEY  WORDS  (Continue  on  It  ntcenary  and  identity  My  btoek  mmbor) 


Aircraft 

Crash  resistance 
Dynamics 

Structural  properties 


Tolerances  (physiology) 
Helicopters 
Energy  absorbers 
Mathematical . models 


I.  ABSTRACT  ( Contlm*  an  ravtraa  a/da  It  n«c««*«ry  and  Identify  by  block  mmbot) 


The  results  of  a  study  to  develop  simplified  design  proceduies  for  deter¬ 
mining  nonlinear  load-deflection  characteristics  for  helicopter  airframe 
structures  are  presented.  A  literature  survey  irr  performed  in  which  60 
technical  reports  and  papers  are  evaluated  with  regard  to  their  applica¬ 
bility  to  crashworthiness  analysis  and  designer-oriented  procedures. 


DD 


I  JAN *75  1473  EDITION  OF  I  NOV  «5  IS  OBSOLETE 

R*'p'  yrtlirAtl  hy  • 


NATIONAL  TFTHNlf  Al  { 

INFORMATION  '  f  RVK  T 


_  Unclassified _ 

SECURITY  CLASSIFICATION  OT  THIS  PAGE  (Vftwi  D»f.  Knloftit) 


|»  „n*  nf  ' 


_ Unclassified _ 

tKCUWTY  CLAgMTtCATMM  OF  THI»  *AM(9km  Dm* 

4.  Continued. 

DEVELOPMENT  OF  SIMPLIFIED  ANALYTICAL  TECHNIQUES  TO  PREDICT  TYPICAL 

HELICOPTER  AIRFRAME  CRUSHING  CHARACTERISTICS  AND  THE  FORMULATION 

OF  DESIGN  PROCEDURES 

20.  Continued. 
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deflection  curves  are  compared  to  the  test  load-deflection  curves  and  show 
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cally. 

Program  KRASH  is  refined  to  facilitate  its  use  by  designers.  In  particu¬ 
lar,  the  input  data  is  reordered,  some  inputs  are  standardized  and  more 
general  load-deflection  curve  characteristics  are  possible.  The  capucity 
of  the  program  is  increased  to  (0  lumped  masses,  100  internal  beams  and  120 
load-deflection  tables.  The  refined  program  was  r’;n  to  demonstrate  capabil¬ 
ity  to  treat  a  three-dimensional  impact  velocity,  mass  penetration  into  an 
occupiable  space,  and  simplified  rotor  blade  contact.  Specimen  test  data 
is  also  iSel  to  refine  tine  31  mass  UH-1H  model. 

IV. e  analytical  techniques  developed  herein  are  presentee  in  the  form  of 
design  ciiarts,  nomographs,  curves,  tables  and  equations  and  form  the  basis 
of  a  structural  crashworthiness  design  manual.  The  besirn  procedures  are 
outline  1  in  a  step-by-step  process  including  examples.  Included  in  the 
design  procedures  are  design  modeling  guidelines  for  such  typical  structural 
elements  as  columns,  beams,  honeycomb  panels  an'1  frangible  tubes. 

Volume  II  contains  supporting  analytical  and  test  data  and  a  literature 
matrix  categorization.  A  descripti  >n  of  refined  program  KRASH  is  provided 
which  S!.  "-’o  the  new  input-output  format,  listing,  and  ..•ample  problems. 
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INTRODUCTION 


BACKGROUND 


To  date,  research  into  the  determination  of  load -deflect ion  characteris¬ 
tics  for  structural  elements  has  been  concentrated  on  the  elastic  region* 
Literally  thousands  of  tests  have  been  run  on  as  many  specimens  vhich 
have  varied  in  material,  size,  shape,  loading,  and  end  fixity.  Semiempir- 
ical  procedures  for  determining  yield  inputs,  elastic  curve  and  types  of 
failure  have  been  generated  for  all  kinds  of  conceivable  structures.  By 
comparison,  there  is  a  sparsity  of  data  available  for  structural  behavior 
beyond  the  elastic  regime.  This,  of  course,  is  not  to  say  that  no  useful 
information  has  been  developed.  A  development  program  could  be  undertaken 
to  obtain  the  same  type  of  data  to  support  analysis  for  nonlinear  behavior 
as  has  already  been  accomplished  on  linear  behavior  of  structural  elements. 
However,  before  such  a  program  is  undertaken,  the  following  factors  should 
be  considered: 


The  cost  in  manpower,  material  and  time  would  be  very  great. 


•  The  data  can  be  obtained  by  performance  of  proof  of  crash- 
worthy  design  tests  as  a  modest  addition  to  the  standard 
laboratory  structural  tests  that  are  conducted  to  qualify  an 
airframe . 

Reference  (l)  provides  a  program  by  which  an  iterative  process  is  used  to 
incorporate  improved  crashworthy  characteristics  into  helicopters  during 
the  early  design  phase.  Since  large  deformations  are  involved  (beyond  the 
elastic  range)  in  determining  the  energy  absorbed  by  the  structure  during 
a  crash,  the  nonlinear  response  of  the  various  significant  structural 
parts  comprising  the  airframe  is  needed.  In  addition,  because  of  the 
large  number  of  different  structures  used  in  airframe  construction,  the 
designer  needs  the  data  in  some  orderly  form  that  categorizes  the  response 
characteristics. 

The  Crash  Survival  Design  Guide  (Reference  2)  provides  valuable  informa¬ 
tion  which,  if  applied  properly,  will  enhance  occupant  survivability  in  a 
severe  yet  potentially  survivable  crash.  The  design  guide,  as  currently 
constituted,  presents  in  detail  aircraft  crash  kinematics,  survivable  im¬ 
pact  envelopes,  alx-frame  crashworthiness  principles,  seat  design,  restraint 
system  design,  occupant  environment  design,  and  additional  applicable  data. 
However,  the  design  guide  does  not  describe  the  following  significant  items: 

•  procedures  by  which  crashworthiness  capability  can  be 
evaluated 


•  design  criteria  as  a  function  of  the  aircraft  design  con¬ 
figuration  and  operational  usage 
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•  load  deflection  determination  procedures  in  the  nonlinear 
regime  for  various  types  of  structural  elements  used  in  air¬ 
frame  design 

In  the  last  decade,  there  have  been  significant  efforts  to  develop  crash- 
worthy  seats  (pilot,  troop  and  gunner).  The  capabilities  of  these  seats 
can  be  greatly  enhanced  if  the  structure  below  the  seat,  namely,  the  land¬ 
ing  gear  and  airframe,  can  be  designed  to  absorb  additional  energy.  An¬ 
other  pertihent  factor  that  must  be  considered  in  crashworthy  design, 
whether  it  be  the  airframe,  seats  or  major  mass  items,  is  that  aircraft 
differ  in  design  configurations,  in  usage  and,  consequently,  in  impact 
environment  and  dynamic  response  characteristics. 

Until  recently  it  has  been  difficult  to  assess  aircraft  crashworthiness 
capability.  However,  the  development  of  digital  computer  programs  (refer¬ 
ences  1  and  82)  has  served  to  provide  basic  analytical  tools.  Program 
KRASH  (reference  l),  in  particular,  has  proved  to  be  capable  of  predict¬ 
ing  structural  and  occupant  response  during  combined  velocity (vertical- 
lateral  and  longitudinal)  impact.  Program  KRASH  provides  a  foundation 
from  which  the  practical  application  of  a  consistent  crashworthiness  ap¬ 
proach  can  evolve.  The  recent  study  described  in  reference  1  established 
that  an  analytical  method  is  available  which  is  a  useful  tool  for  the  de¬ 
signer,  if  guidelines  for  the  determination  of  nonlinear  load-deflection 
characteristics  of  structural  elements  used  in  aircraft  design  are  pro¬ 
vided.  The  development  of  simplified  methods  to  obtain  this  type  of  data 
will  result  in  the  establishment  of  procedures  which  can  readily  be  applied 
to  various  airframes  to  evaluate  their  crashworthiness  capability.  Further¬ 
more,  designers  of  aircraft  will  have  the  means  to  improve  airframe  crash- 
worthiness  in  an  orderly  manner  to  a  desired  level,  with  a  minimum  incre¬ 
mental  weight  and  cost  increase.  Since  the  development  of  procedures 
applicable  to  crash  analysis  is  primarily  intended  for  designers,  with 
special  emphasis  on  the  preliminary  design,  special  consideration  is  given 
to  simplified  expressions  which  have  general  application  as  opposed  to 
rigorous  approaches  which  are  oriented  toward  unique  design  or  loading 
conditions. 

There  are  many  different  types  of  structural  elements  used  in  the  various 
sections  of  a  helicopter.  The  structure  is  designed  to  take  loads  acting 
in  different  directions;  consequently,  different  modes  of  failure  occur 
such  as  buckling,  tension,  compression,  bearing  and  shear.  Table  I  gives 
a  sampling  of  the  types  of  structural  elements  and  their  locations  in  a 
UH-1H  helicopter.  Some  useful  load-deflection  information,  both  linear  and 
nonlinear,  has  been  developed  over  a  period  of  time  for  many  of  the 
different  types  of  structural  elements  used  in  a  helicopter.  In  particular, 
much  experimental  and  analytical  effort  has  been  expended  to  determine 
elastic  responses  of  columns,  struts,  stringers,  longerons,  honeycomb 
panels  and  beams  to  dynamic  and  static  loads. 
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While  research  has  been  oriented  predominately  in  the  linear  regime,  the 
results  provide  valuable  information  (i.e.,  the  yield  points,  the  elastic 
curve  and  buckling  behavior)  for  establishing  the  characteristics  of 
certain  structures  when  the  nonlinear  region  is  encountered. 


taw.v  t  ttw-iw  stpu 

CTURAL  ELEMENT  TYPES  AND  LOCATIONS 

1.  Honeycomb  panel 

Main  beams,  bulkheads,  skin,  fuel 

cell,  door,  etc. 

2.  Torque  box  (beams) 

Fuselage  mid  section 

3.  Simply  supported  beam 

Lift  link  (transmission) 

4.  Cantilevered  beam 

Fwd.  of  f.s.  155 »  aPt  of  f.s.  211 

5 .  Tension  rod 

Fuselage  (f.s.  211  to  243) 

6.  Flanged  webs 

Bulkheads,  main  beam 

7.  Frames 

Cabin,  tail  boom 

8.  Stanchions,  columns 

Cabin 

9.  Strut,  bipod,  tripod 

Engine  mount 

10.  Bulkheads 

Fuselage,  tail  boom 

11.  Stringers 

Tail  boom 

12.  Longerons 

Tail  boom 

13.  R^b  spar 

Vertical  fin 

l4.  Torsion  tube 

Horiz.  stabilizer,  engine-transmission 

15.  Panel  (flat,  curved) 

Fuselage 

l6.  &cin  stiffeners 

Fuselage,  tail  boom 

17.  Tubular  (curved,  flat) 

Landing  gear,  troop  seats 

18.  Rubber  mounts 

Transmission 

19.  Beam  stiffeners 

Main  beam 

An  important  consideration  in  defining  the  energy  absorption  characteristics 
of  .structural  elements  in  the  nonlinear  region  is  to  establish  the  signifi¬ 
cance  of  rate  of  loading.  To  this  end, the  energy  absorption  characteristics 
of  a  typical  aircraft  substructure  (P2V-4  fuselage  bumper)  were  obtained 
under  static  and  dynamic  loading  conditions  (Reference  l).  The  results  of 
these  tests  are  compared  in  Figure  1.  The  dynamic  test  results  show  good 


agreement  in  the  deformation  of  the  shell  structure  when  compared  to  the 
static  test  data.  The  results  obtained  from  the  tests  of  the  fuselage 
bumper  are  significant  in  that  the  dynamic  test  of  the  bumper  indicates 
that  for  some  typical  aircraft  structure,  a  simple  static  test  will 
suffice . 


ENERGY, 

LOAD,  LB 

IP. -LB 

DEFLECTION  UNDER  LOAD 


Figure  I.  Load-Deflection  and  Energy-Absorbed  Results  for  Static 
and -Dynamic  Fuselage  Bumper  Tests  (Reference  l). 


Souse  structures,  on  the  other  hand,  will  behave  differently  under  dynamic 
loading  conditions  than  under  a  static  load  environment.  However,  if 
static  tests  can  be  performed  in  lieu  of  dynamic  tests  to  obtain  load- 
deflection  characteristics,  then  significant  savings  can  be  realized  in 
future  programs.  What  is  needed  is  an  approach  which  recognizes  that 
existing  linear  and  nonlinear  load  deformation  data  can  formulate  a  basis 
for  establishing  guidelines  which  will  help  designers,  during  the  prelimi¬ 
nary  design  stage  of  a  vehicle's  development,  to  determine  the  amount  of 
energy  absorption  necessary  bo  achieve  a  crashworthy  design.  The  exact 
load-deflection  curve  of  the  element  which  deforms  plastically,  although 
desirable,  is  not  essential  to  a  crashworthiness  analysis.  More  important 
is  the  understanding  of  the  gross  behavior  in  the  plastic  region,  i.e., 
constant  load-carrying  capacity  as  compared  to  instantaneous  loss  of  load¬ 
carrying  capability. 

From  an  economic  point  of  view,  it  is  more  efficient  to  optimize  for  the 
incorporation  of  crashworthiness  features  early  in  the  design  stage.  How¬ 
ever,  during  preliminary  design  it  is  doubtful  if  sufficient  design  details, 
which  influence  the  exact  plastic  deformation  shape  of  structural  elements, 
will  be  available.  The  availability  of  simple  procedures  to  predict  energy 
absorption  and  load-deformation  characteristics  will  allow  the  designer  to 
initiate  valuable  cost,  weight  and  geometry  trade-off  studies.  The  develop¬ 
ment  of  these  procedures  will  require  some  testing  of  typical  specimens. 

This  testing  should,  as  a  minimum,  verify  the  validity  of  proposed 
procedures  for  providing  pertinent  nonlinear  load-deformation  data. 

PROGRAM  OBJECTIVES 

The  objectives  of  the  program  are  as  follows: 

•  To  develop  simplified  analytical  techniques  by  which  nonlinear 
load-deflection  characteristics  of  typical  structural  elements 
can  be  obtained. 

•  To  verify  with  tests  that  these  simple  analytical  techniques 
are  valid  for  crashworthiness  analysis. 

•  To  provide  designers  with  simplified  techniques  with  which  to 
perform  crashworthiness  studies  during  the  preliminary  design 
phase . 

•  To  simplify  data  input  requirements  of  computer  program  KRASH. 

To  achieve  the  objectives,  the  results  of  the  study  will  determine: 

(l)  The  pertinent  structural  element  load-deflection  characteristics 
that  are  to  be  modeled  in  order  to  obtain  acceptably  accurate 
structural  and  occupant  responses. 


5 


‘j.X:*.^y  y  1  ■ 


(HKsessw* 


inw#.  m*mt*-~ 


(2)  The  accuracy  with  which  typical  airframe  structure  load-deflection 
characteristics  can  be  predicted  using  simplified  analytical  pro¬ 
cedures. 


BEPQRT  FORMAT 


The  table  of  contents  conveys  the  general  plan  of  the  report.  Background 
information  is  presented  initially  to  provide  a  proper  perspective  in  re¬ 
lation  to  the  objectives  of  the  program.  Volume  I  is  presented  in  chrono¬ 
logical  order  of  the  program  milestones  and  includes: 


Survey  of  Technical  Publications 
Load  Sensitivity  Studies 
Substructure  Analysis 
Substructure  Test  Program 
Correlation 

Program  KRASH  Refinement 
Design  Procedures 


The  program  results  are  presented  in  summary  form  prior  to  the  conclusions. 


Volume  II  contains  supporting  data  for  the  details  presented  in  Volume  I, 
Volume  II  contains  abstracts  for  the  60  technical  reports  and  publica¬ 
tions  reviewed  during  the  program.  Included  in  the  literature  survey 
section  is  a  matrix  categorization  of  the  reports  by  subject  and  applic¬ 
able  area  of  interest.  Us  test  data  is  presented  in  its  entirety,  in¬ 
cluding  accelerometer,  load  cell,  strain  gage  deflection,  and  time  (or 
scan)  histories  for  each  of  the  twelve  tests.  Supporting  data  for  anal¬ 
ysis  and  design  procedures  is  also  presented.  Program  KRASH  refinements 
which  include  standardization  of  input,  modified  input  format,  and  user's 
guide  are  included. 


.  i 


im 


SURVEY  OP  TECHNICAL  PUBLICATIONS 


GENERAL  DISCUSSION 


Technical  publications  were  reviewed  and  their  contents  categorized  in  accor¬ 
dance  with  the  various  aspects  of  developing  improved  crashworthiness 
criteria  and  concepts  to  assist  in  future  reference  as  needed.  A  review  and 
evaluation  of  60  publications  are  presented  herein.  A  complete  listing  of 
publications  which  comprise  the  literature  survey  is  presented  in  Volume  II. 
The  literature  is  reviewed  with  regard 'to  the  impact  on  developing  designer- 
oriented  crashworthiness  design  procedures.  In  particular,  emphasis  is 
placed  in  the  following  areas: 


•  Simplified  analytical  techniques 

•  Experimental  data 

•  Design  procedures  and  guidelines 


•  Load-deflection  data 


Energy  absorption 


•  Methods  of  analysis 


9  Strain-rate  sensitivity 


•  Inertia  effects 


EVALUATION 


All  documents  reviewed  are  evaluated  for  their  applicability  to  each  of  the 
subject  areas  listed  above,  the  contribution  of  each  is  briefly  stated,  and 
a  composite  summary  of  the  pertinent  aspects  is  recorded. 


Volume  II,  which  provides  additional  literature  survey  information,  con¬ 
sists  of  a  summary  or  abstract  for  each  report  and  a  literature  survey 
subject  index  which  contains  a  matrix  categorization  of  the  contents  of 
the  technical  reports.  The  matrix  categorization  lists  the  reports  by 
number  and  associates  each  report  with  an  area  of  specific  content  or 
applicability.  All  report  numbers  in  this  section  refer  to  the  reference 
numbers  as  they  are  listed  under  Literature  Cited. 


Simplified 


The  determination  of  load  deformation  characteristics  of  aircraft  structure 
will  enhance  the  ability  of  the  designer  to  predict  structure  and  occupant 
responses  during  severe  yet  suivivable  crashes.  However,  it  is  important 
that  this  data  is  presented  in  a  manner  which  will  aid  the  designer  in 
developing  the  desired  level  of  crashworthiness  in  the  airframe  structure. 
The  data  presented  in  the  literature  provides  analytical  procedures  and/or 
empirical  data  which  is  applicable  co  a  great  many  types  of  structural  ele- 
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ments.  Unfortunately,  much  of  the  data  in  its  present  form  is  not  direct¬ 
ly  applicable  for  the  following  reasons: 

(1)  The  analytical  techniques  require  computer  programs  in  order  to 
formulate  solutions. 

(2)  The  empirical  data  is  not  related  to  analytical  procedures. 

(3)  The  data  is  for  structural  elements  under  a  unique  loading  con¬ 
dition  or  boundary  condition. 

(4)  The  data  neglects  effects  of  inertia  loads,  strain  rate,  wave 
propagation  and/or  geometry  changes. 

There  are  several  references,  most  notably  5,  6,  7>  9 ,  12,  18,  24  and  30, 
■which  provide  information  that  can  be  incorporated  into  procedures  that 
could  serve  as  guidelines  for  designers.  Several  of  the  aforementioned 
references  are  texts  (6,  7,  8,  9,  24,  30)  which,  although  primarily  limit¬ 
ed  to  elastic  and  plastic  behavior,  cover  a  wide  range  of  structural  ele¬ 
ments  such  as  beams,  columns,  plates,  rings,  arches  and  composite  structure. 
The  data  in  these  references  provide  basic  information  regarding  yield 
point  loads  and  methods  by  which  the  effects  of  plasticity,  crippling  and 
buckling  can  be  taken  into  account. 

Experimental  Data 

References  1,  4,  17,  19,  26,  28,  32,  33  and  35  present  results  of  tests. 

The  data  is  generally  in  the  form  of  load  versus  d  .lection,  force  versus 
time  and  stress  versus  strain.  Reference  1  static  and  dynamic  test  re¬ 
sults  shoved  that  for  some  typical  aircraft  structure,  the  load  deflection 
characteristics  are  the  same  for  both  static  and  dynamic  loading  condi¬ 
tions.  Figure  1  shows  a  comparison  of  the  results  of  the  static  and  dyna¬ 
mic  tests  described  in  Reference  1. 

Ir.  Reference  4,  load  deflection  curves  for  plate  stringer  panel  configura¬ 
tions  are  presented.  The  results  of  these  panel  tests,  illustrated  in 
Figure  2,  show  the  energy  to  area  ratio  of  the  three  configurations.  Of 
particular  importance  is  the  relative  absorption  capabilities  of  the  panels 
and  their  respective  modes  of  failure.  For  example,  the  integrally 
machined  panels  produce  a  reasonably  efficient  energy-area  ratio,  but  the 
load  stroke  performance  is  poor  because  the  mode  of  failure  is  an  explo¬ 
sive  fracture  in  which  the  riser  splits  completely  off  the  skin.  In  addi¬ 
tion,  the  peak  load  of  this  latter  configuration  is  much  higher  than  the 
peak  load  for  the  other  panels,  which  could  result  in  higher  transmitted 
loads.  Reference  (4)  also  presents  load-deflection  data  for  the  lower 
frame  of  a  typical  fuselage  segment  obtained  from  drop  tests.  This  test 
data  is  compared  to  analysis  and  is  shown  to  differ  substantially.  Refer¬ 
ence  10  presents  data  regarding  buckling  characteristics  of  perfect  and 
imperfect  circular  cylindrical  shells  subjected  to  dynamic  axial  loading. 
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Figure  2.  Plate  Stringer  Compression  Panels  Load-Deflection  Test  Data 
(Reference  4). 


Close  agreement  was  found  between  theory  and  experiment  for  both  dynamic 
buckling  strength  and  buckling  mode  shapes. 

Reference  17  presents  data  for  simple  structures  like  a  clamped  beam,  cir¬ 
cular  ring,  circular  plate  and  hemisphere.  The  data  is  limited  to  deflec¬ 
tion  time  histories  and  is  for  a  specific  material,  loading  condition, 
boundary  condition  and  size  of  structure.  In  Reference  19,  a  comparison 
is  made  of  experimental  data  with  the  results  of  a  theoretical  investiga¬ 
tion  of  the  plastic  deformation  of  cantilevered  beams  subjected  to  impul¬ 
sive  loading.  The  test  results  for  this  particular  study  followed  the 
general  trend  predicted  by  rigid-plastic  theory.  A  conclusion  drawn  from 
the  results  is  that  an  increase  in  yield  stress  with  strain  rate  is  the 
primary  cause  of  the  discrepancy  between  theory  and  experiment. 

Reference  26  presents  results  of  an  experimental  investigation  into  the 
use  of  frangible  metal  tubing  as  an  energy  absorber.  Tests  were  performed 
using  a  l/5  scale  model  of  a  proposed  manned  spacecraft,  and  impact  veloci¬ 
ties  reached  30  fps  vertical  and  combined  vertical  (13*5  fps)  and  horizon¬ 
tal  (l8  fps).  Results  are  presented  as  curves  of  force  versus  displace¬ 
ments.  References  27,  32,  33,  34  and  35  present  data  for  a  one-shot 
energy  absorber.  References  32  through  35  provide  data  accumulated  over 
a  period  of  years  for  the  Lockheed  developed  Dynasorb.  The  device  has 
been  tested  for  ground  impact  conditions  as  high  as  112  fps.  Materials 
tested  include  aluminum  alloy  2024-T3,  brass,  copper,  magnesium  AZ-31B-F, 
steel  (1015),  steel  (^130),  and  titanium.  Energy  curves  showing  load  ver¬ 
sus  displacement  are  presented. 

Reference  °9  presents  results  of  experimental  investigations  of  thin 
webbed  plate-girder  beams  to  determine  shear  loads  required  to  cause  web 
rupture  and  flange  rivet  failure.  In  all,  twenty- seven  beams  were  tested. 
Reference  38  discusses  the  impact  response  of  curved  box-beam  columns 
with  large  global  and  local  deformations.  Of  significance  are  load- 
deflection  curves  for  steel  and  aluminum  under  static  and  dynamic  loading 
conditions  (Figure  3)*  This  paper  concludes  that  the  consideration  of 
cross-section  changes  is  necessary  and  important  in  predicting  the  impact 
response  of  beam  columns  with  thin-walled  box  sections  subjected  to  large 
deformation.  The  effects  of  strain  rate  sensitivity  and  strain  wave  prop¬ 
agation  on  the  load-deflection  curve  are  discussed. 

Reference  39  presents  stress-strain  curves  for  aluminum  honeycomb  and 
two  foamed  plastic  structures.  The  characteristics  of  these  three  materials 
are  examined  to  ascertain  how  they  can  be  applied  for  human  protection 
against  accelerations  encountered  at  low  impact  speeds  (30  fps).  Reference 
4l  presents  data  regarding  modes  of  failures  of  multiweb  beams.  The  ulti¬ 
mate  strength  and  buckling  characteristics  of  multi web  beams  have  been  in¬ 
vestigated  both  experimentally  and  theoretically.  The  three  primary  types 
of  instability  that  occur  are  (l)  local  buckling,  (?)  wrinkling,  and  (3)  in¬ 
terrivet  buckling.  It  is  noted  that  beams  of  solid  cross  sections  typically 
ext  Hit  large  plastic  deformations  as  shown  in  Figure  t,  while  suilt-up 
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of  Impact  Speed  on  Dynamic  Force-Deflection 
i  of  Curved  Box- Beam  Columns  (Reference  38). 
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Figure  4.  Structural  Behavior  Typical  of  Solid  Cross-Section 


Beams  of  Ductile  Material  (Reference  39)* 
BUCKLING 


Figure  5.  Structural  Behavior  Typical  of  Built-Up  Beams  (Reference  39) 


II 


sheet  beams  exhibit  structural  behavior  similar  to  that  shown  in  Figure  5» 
Design  Procedures  and  Guidelines 

To  perform  a  crash  analysis  of  an  aircraft  requires  the  use  of  verified 
analytical  tools,  in  addition  to  the  definition  of  the  crash  environment 
and  human  tolerance  considerations.  Potentially,  the  greatest  payoff 
with  regard  to  efficiency  of  design  in ' performing  cbash  analysis  exists 
during  the  preliminary  design  stages.  However,  during  the  initial  stages 
of  design,  little  data  is  available  other  than  basic  concepts,  configura¬ 
tion  and  sizing.  Thus,  the  designer  is  faced  with  the  formidable  task  of 
incorporating  crashworthy  design  features  which  are  compatible  with  the 
crash  environment,  aircraft  operational  requirements,  cost,  weight  and 
geometry  constraints  while  working  with  preliminary  design  data.  Of  use 
to  the  designer  are  several  categories  of  information.  For  example,  de¬ 
sign  principles  which  include  techniques  and  considerations  for  improving 
crashworthiness  capability  provide  one  set  of  guidelines.  References  2, 

3  and  4  are  particularly  useful  in  this  respect.  Reference  1  outlines 
an  approach  which,  based  on  current  available  data  and  techniques,  pro¬ 
vides  an  opportunity  to  assess  the  trade-off  between  the  level  of  crash- 
worthiness  capability  versus  incremental  increase  in  cost  and  weight  to 
achieve  improved  crashworthiness  capability.  It  also  presents  a  systema¬ 
tic  approach  for  ascertaining  crashworthiness  capability  which  can  incor¬ 
porate  the  data  compiled  by  various  agencies  and  studies  over  the  years. 

In  addition.  Reference  1  provides  a  verified  analytical  method  which  can 
form  the  basis  of  a  designer-oriented  tool. 

References  6,  7,  12,  13 ,  29,  30  and  34,  some  of  which  are  texts,  provide 
procedures  by  which  load  deflection  characteristics  of  various  types  of 
structures  in  the  elastic  range  can  be  obtained.  In  addition,  these  refer¬ 
ences  show  how  the  effects  of  plasticity  and  crippling  can  be  incorporated 
These  latter  references  present  data  for  a  variety  of  types  of  aircraft 
structures  and  include  changes  in  geometry  and  boundary  conditions  for  the 
different  structures.  In  general,  the  data  presented  from  the  publica¬ 
tions  in  the  literature  is  not  oriented  toward  design  type  personnel; 
accordingly,  little  in  the  way  of  procedures  and/or  guidelines  is  made 
available . 

Load-Deflection  Data 

Although  much  effort  has  been  expended  to  obtain  structural  characteris¬ 
tics  of  structural  elements,  only  a  small  percentage  of  published  data 
provides  load-deflection  data  that  can  be  directly  incorporated  by  de¬ 
signers.  Table  II  shows  a  cross  section  of  data  available  from  the  refer¬ 
ences  that  are  included  in  the  literature  evaluation.  The  load-deflec¬ 
tion  data  provides  some  useful  data  to  the  designer,  in  that  characteris¬ 
tic  trends  can  be  associated  with  differei  c  types  of  structural  elements. 
Table  III  presents  a  matrix  of  structural  element  load-deflection  cate¬ 
gories  and  load  types.  The  table  indicates,  to  some  degree,  how  the 
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TABLE  II.  LOAD- DEFLECTION  DATA 


)  CATEGORIES  ** Data  In  form  of  stress  versus  strain  and 

Increasing  load  with  increase  in  deflection  load  versus  time 

Constant  load  with  increase  in  deflection  ***  Inadequate  control  ring  design  results  ii 
Decreasing  load  with  increase  in  deflection _ load  increase  versus  deflection _ 


TABLE  in.  MATRIX  0?  STRUCTURAL  ELEMENT  LOAD-DEFLECTION  CATEGORIES  AND 

LOAD  TYPE 

Load 

Type 

Load  Categories  * 

- 1 

2 

3  ** 

Compression 

Skin-Stringer 

Coil  Spring 

Stiffened  Cylinder 

Frangible  Tube 

Telescoping  Tube 

Honeycomb 

Multiple  Cell 
Structure 

Short  Column 

Sheet  Skin 

Flat  and  Curved  Plate 

Axially  Loaded  Cylinder 

Slender  Long  Column 

Stringer  Shell 

Segment 

Beam 

Tension 

Bulkhead 

Stiffener 

Inversion  Tube 

Stainless  Steel 
Strap 

Stranded  Cable 

Short  Elongation  Strap 

Sheet  Skin 

Stringer 

*  LOAD  CATEGORIES 

1.  Increasing  load  with  increase  in  deflection 

2.  Constant  load  with  increase  in  deflection 

3.  Decreasing  load  with  increase  in  deflection 

**  Compression  members  tend  to  fail  as  a  result  of  the  lateral  bending 
induced  by  the  compression  load,  an  action  which  is  commonly  called 
buckling. 
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selection  of  a  structural  element  design  can  influence  the  loads  that  will 
he  experienced  in  an  aircraft.  This  information,  without*  additional  data 
such  as  energy-absorption  efficiency,  linear  load-deflection  curves,  and 
yield  points  for  specific  design  configurations ,  although  a  contribution, 
is  inadequate  for  crashworthiness  design.  References  6,  7  and  30  present 
some  valuable  data  concerning  various  types  of  structural  elements  up  to 
the  point  of  failure  and  including  plasticity  effects.  Of  particular 
importance  in  this  data  is  that  nondimensional  graphs  are  presented;  these 
allow  one  to  easily  determine  load  capability  for  different  geometry  or 
end  constraints.  Reference  6  includes  discussions  of  plastic  bending, 
buckling  of  flanges  and  webs,  lateral  buckling  of  beams,  buckling  of  beams 
in  combined  axial  compression  and  bending,  buckling  of  frames,  and  a  gener¬ 
al  method  for  computing  elastic-plastic  displacements.  Reference  7  pre¬ 
sents  discussions  and  tables  of  data  for  determining  load  and  deflection 
for  beams  and  frames,  including  simultaneous  axial  and  transverse  loading, 
variable  cross  sections  and  curved  sections,  flat  plates  (ultimate  strength, 
large  deflections,  nonuniform  loading),  columns  (eccentric  loading,  com¬ 
bined  compression  and  bending),  and  buckling  of  bars,  columns,  flat  and 
curved  plates.  Reference  30  provides  a  discussion  of  the  design  of  members 
in  tension  (skin,  stringers,  spar  caps),  bending  (beams),  torsion  (shafts), 
compression  (columns,  flat  plates,  curved  sheets),  and  shear  (webs).  In¬ 
cluded  in  this  reference  are  design  equations,  nondimensional  buckling  de¬ 
sign  curves,  and  a  discussion  of  the  manner  in  which  local  crippling  fail¬ 
ures  can  be  analyzed. 


Prom  a  crashworthiness  design  efficiency  standpoint,  it  is  desirable  to 
have  structural  elements  which  exhibit  a  load- deflection  curve  that  is 
approximately  flat  (zero  slope)  beyond  the  elastic  limit,  provided  the 
maximum  transmittal  load  does  not  subject  the  occupant  to  an  intolerable 
acceleration  magnitude.  There  are  many  load-limiting  devices  which  can  be 
incorporated  for  crashworthiness  considerations.  However,  the  selection 
of  a  particular  device  requires  that  a  trade-off  be  made  between  perfonnance, 
weight,  space  and  cost.  Table  IV,  obtained  from  Reference  2,  compares 
several  "one-shot"  load-limiting  devices  for  1000-to-4000-pound  loads.  In¬ 
cluded  is  a  comparison  of  the  pertinent  design  factors  for  eight  different 
devices.  "Long-term  reliability"  refers  to  the  ability  of  the  device  to 
perform  its  function  without  benefit  of  maintenance  throughout  the  life  of 
the  aircraft,  while  "specific  energy"  indicates  the  amount  of  energy  that 
can  be  absorbed  per  pound  of  weight  of  the  device. 

References  1,  3,  *+,  17,  25,  26,  27,  32,  33,  3*+  and  35  also  present  data 
related  to  structural  element  energy-absorption  requirements  or  capability. 
Table  V  (Reference  33),  for  example,  compares  energy  absorption  efficiencies 
for  several  different  materials  and  methods. 

References  32  through  35  discuss  the  Lockheed  Dynasorb  frangible  tube  ener¬ 
gy  absorber  which  has  been  tested  using  a  3400-pound  simulated  vehicle 
structure  with  a  ground  impact  velocity  of  112  feet/second. 
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TABLE  V.  COMPARATIVE  SUMMARY  OF  ABSORPTION  ENERGIES  ( 


The  several  methods  of  analysis  that  are  available  for  determining  load- 
deflection  data  for  typical  structural  elements  include: 

•  Classical  static  plastic  analysis 

•  Dynamic  mode  approximations 

•  Finite  differences 

•  Finite  elements 

Classical  static  plastic  analysis  is  applicable  to  most  structures  except 
wherein  instability  and/or  buckling  occurs.  Beam  and  simple-shell  type 
structures  are  generally  analyzed  by  this  method.  The  dynamic  mode  approxi 
mation  technique  is  described  in  References  45  and  6l.  This  approach  is 
relatively  simple  to  apply;  however,  it  does  have  the  following  limitations 

•  A  constant  acceleration  is  assumed  throu^iout  the  dynamic  de¬ 
formation.  Therefore,  the  use  of  load-deflection  characteris¬ 
tics  obtained  by  this  approach  would  not  accurately  reflect 
the  actual  acceleration  response  of  the  structure. 

•  This  technique  is  not  applicable  to  structure  which  exhibits 
instability  and/or  buckling  characteristics. 

•  For  complex  structures  it  is  difficult  to  define  a  unique 
mode  shape  which  is  kinematically  and  dynamically  admissible. 

Finite  difference  and  finite-element  techniques  are  discussed  in  References 
17,  38*  52  and  63.  These  approaches  give  reasonable  load-deformation  in¬ 
formation  for  simple  structures.  However,  for  complex  structures, these 
numerical  methods  are  not  yet  considered  as  reliable  or  as  economical  as 
experimental  techniques.  Analytical  techniques  and  numerical!  methods 
Tor  stability  analysis  and  post-failure  analysis  presented  in  References 
15  and  64  through  71  are  primarily  concerned  with  the  technique  itself 
rather  than  the  application  of  the  technique  to  practical  problems.  The 
analytical  approaches  presented  in  the  literature  are  devoid  of  methods 
by  which  designers  can  make  direct  use  of  the  analytical  procedures. 

•Strain  Rate  Effects 

If  materials  exhibit  strain-rate  sensitivity,  then  the  effect  on  the  load- 
stroke  characteristics  for  the  structure  will  be  a  change  in  the  yield  and 
ultimate  loads.  For  cases  in  which  structural  materials  exhibit  strain 
rate  sensitivity,  an  approximate  accounting  for  this  effect  can  be  made. 
Several  references,  including  5,  19,  31,  36,  37,  38,  40,  43,  44  and  46, 
discuss  strain  rate  effects.  There  is  no  universally  validated  and 
accepted  strain  rate  law.  One  approximation  described  in  Reference  38  and 
used  in  several  publications  is  shown  on  the  next  page. 


static  value  of  yield 
point 

yield  point 


constants  of  the  struc¬ 
tural  material 


strain 
strain  rate 


on  average 


The  approximation  regards  the  effect  of  strain  rate  as  raising  the  yield 
point  (fry)  above  the  static  yield  value  (fr  )  with  the  associated  strain 
hardening  portion  of  the  curve  kept  parallel  to  the  static  strain  harden¬ 
ing  curve. 


In  Reference  18  is  presented  a  simplified  method  for  solving  impulsively 
loaded  structure  for  rate-sensitive  materials.  The  results  of  the  study 
indicate  that  good  approximations  to  the  exact  solution  may  be  found  by 
utilizing  a  rate-insensitive  material  with  constant  yield  stress  equal  to 
the  initial  dynamic  yield  stress. 


Inertia  Effects 


If  inertia  forces  for  a  particular  type  of  structure  are  important,  then 
the  load-deflection  characteristics  for  the  structure  under  dynamic  load¬ 
ing  conditions  can  differ  from  the  load-deflection  characteristics  under 
static  loading  conditions.  References  1,  8,  24,  29,  38,  40  and  43  present 
data  related  to  inertia  effects.  In  Reference  1  the  load-deflection 
curves  are  obtained  under  dynamic  loading  and  compared  to  a  previously 
obtained  static  load-deflection  curve  for  the  same  specimen  type.  The  load- 
deflection  curve  for  the  specimen  is  also  obtained  analytically  using  an 
existing  finite  element  program. 

If  the  higher  modes  of  the  frequency  components  dominate  the  response 
initially  such  as  to  cause  failure  before  the  fundamental  mode  can  respond, 
then  the  effects  of  inertia  should  be  considered  in  the  load-deflection 
curves. 

Reference  38  shows  the  effects  of  impact  speed  on  the  dynamic  force  deflec¬ 
tion  relationship  for  steel  and  aluminum  curved  box-beam  columns. 
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LOAD- SENSITIVITY  STUDY 


GENERAL 

The  results  of  the  study  described  in  Reference  1  show  that  an  analytical 
procedure  is  capable  of  predicting  the  responses  of  airframe  structure 
during  crash  conditions.  The  results  of  the  study  also  indicate  that  the 
trade-off  between  cost  and  weight  increase  versus  improved  crashworthiness 
capability  is  an  integral  part  of  the  iterative  process,  for 'determining  a 
crashworthy  design.  Although  the  computer  program  developed  in  Reference 
1  provides  an  acceptable  analytical  tool  for  determining  crashworthiness 
capability  during  preliminary  design,  refinements  can  be  made  which  will 
make  the  program  more  designer  oriented.  The  purpose  of  the  load  sensiti¬ 
vity  studies  is  to  ascertain  how  accurate  the  modeling  representation  of 
the  structure  has  to  be  in  order  to  obtain  sufficiently  accurate  results 
fo"  preliminary  design  usage.  It  is  important  to  bear  in  mind  that  during 
a  preliminary  design, the  structural  data  that  is  generally  available  is 
related  to  concepts,  sizing,  and  location,  whereas  very  little  in  the  way 
of  detail  construction  is  known. 

The  load  sensitivity  studies  make  use  of  the  31  lumped  mass  model  of  the 
UH-1H  helicopter  described  in  Reference  1.  The  general  arrangement  of  the 
UH-1H  is  shown  in  Figure  6.  The  basic  math  model  internal  springs,  external 
springs  and  lumped  masses  are  noted.  The  detailed  lumped  mass  model  and  the 
mass  representation  geometry  are  shown  in  Figure  7  and  Table  VI,  respectively. 
For  detailed  information  regarding  the  model,  see  pages  8l-93»  Volume  I  of 
Reference  1.  For  information  regarding  the  related  theory  including 
coordinate  systems,  see  pages  25-86,  Volume  II  of  Reference  1. 

LOAD-DEFLECTION  VARIATIONS 


Engine 

The  changes  in  the  engine-mount  stiffness  factors  and  load-deflection  curves 
are  shown  in  Table  VII  and  Figure  8,  respectively.  The  load-deflection 
curve  is  obtained  from  the  expression 


K 


where 


X 


(KRj)  •  (DX^)  versus  X 


stiffness,  lb/in. 
incremental  distance,  in. 

stiffness  reduction  factors;  i.e.,  for  the  first 
row  in  Table  VII,  the  KR  vs  DX  curve  appears  as 
shown  in  the  sketch  accompanying  Table  VII. 

distance,  in. 
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General  Arrangement  of  the  UH-1H  Helicopter. 


Figure  7.  31  Lumped  Mass  Helicopter  Model  Representation  (Reference 


•PAX  denotes  passenger. 


The  KR  versus  X  curve  is  depicted  below : 


There  are  six  KR  values  for  each  element.  The  first  KR  value  corresponds 
to  a  zero  deflection.  The  incremental  deflections  (DX)  are  equally  spaced. 
The  KR  values  shown  in  the  sketch  above  are  for  the  base  axial  case  (31-52) 
shown  in  Table  VII. 

For  a  linear  spring,  KR  =1 

For  a  constant  force,  KR  =0 
For  a  decreasing  force,  KR  <0 


TABLE  VII.  ENGINE-MOUNT  AXIAL  AND  BENDING  KR  CHANGES 


Case 

Identification 

DX 

(in.) 

KR  Values 

at  DX.. 

© 

© 

© 

© 

© 

© 

AXIAL 

Base 

.1 

1.0 

1.0 

-.7 

-.7 

0 

0 

El 

.1 

1.0 

.8 

0 

-.5 

-.25 

0 

E2 

.1 

1.0 

.8 

0 

0 

0 

0 

E3 

.1 

1.0 

.6 

.5 

.4 

.3 

.2 

BENDING 

Base 

.0667 

1. 

1. 

1. 

0 

0 

EB1 

,1 

1. 

1. 

-  1 

0 

-1. 

-.8 

EB2 

.1 

1. 

l. 

. 

0 

1. 

.8 

.1 

1. 

1. 

1 

0 

-.5 

-.4  • 

M L- 

.1 

1. 

1. 

1 

0 

_ 

.4 

Kaxial  =  5*64  X  ^  lb/in* 
^b  ending  ~  ^  x  1°>  lb/ir. 
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The  engine  mount  is  a  bipod  arrangement;  thus,  the  vertical  force  is  in¬ 
fluenced  primarily  by  the  member  axial  and  vertical  forces  and  the  angle 
of  the  mount.  A  left  side,  looking  inboard,  view  of  the  engine  mount  arrange¬ 
ment  is  shown  in  Figure  9» 


Figure  9-  heft  Side,  Looking  Inboard, View  of  Engine  Mount. 


Transmission 

The  transmission  mount  is  a  near-vertical  member;  therefore,  the  vertical 
response  is  influenced  primarily  by  the  member  axial  (vertical)  stiffness. 
The  changes  in  the  transmission-mount  stiffness  reduction  factors  and, 
consequently,  the  load-deflection  curve  are  shown  in  Table  VIII.  The 
associated  load-deflection  curves  (K./*  (KR)  (DX)  vs  X)  for  the  transmission 
mount,  members  8  -  9>  are  shown  in  Figure  10. 

Landing  Gear 

The  changes  in  the  landing  skid  stiffness  reduction  factors  are  shown  in 
Table  IX.  The  associated  load-deflection  curves  for  the  forward  landing 
skid  (members  l6  -  17  and  l6  -  18)  and  the  aft  landing  skid  (members  10  - 
it  and  10  -  15)  are  shown  in  Figure  11. 


TABLE  VIII.  TRANSMISSION-MOUNT  AXIAL  AND  BENDING  KR  CHANGES 


Case 

Identification 


AXIAL 


Base 


DX 

(in.) 


KR  Value  at  D. 


©  © 


85  .00785 
85  .00785 


BENDING 


Base 

TB1 

TB2 


K  axial 

=  2  x  105  lb /in. 

K  bending 

=  1.9  x  10k  lb/in. 

TABLE  IX.  LANDING-GEAR  BENDING  AND  ROTATIONAL  KR  CHANGES 


Case 

Identification 

DX 

(in.) 

BENDING 

Base 

4.0 

SBi 

4.0 

SB2 

4.0 

SB  3 

4.0 

SB4 

4.0 

ROTATION 

Base 

.02 

SRI 

.005 

SR2 

.005 

KR  Value  at  DX. 


K  bending 
K  rotation 


2.7  x  10 J  lb/in. 

2.17  x  106  ,in<-lb/rad 


w 


The  effective  bending  stiffness  of  the  landing  skids  can  be  approximated 
from  the  expression 

K  (V2 

*55  33  *55 

where  =  bending  stiffness  (load  due  to  deflection) 

Kjjtj  =  rotational  stiffness  (torque  due  to  rotation) 

=  load  due  to  rotation 

In  program  KRASH,  stiffness  reduction  factors  are  input  separately  for  all 
the  diagonal  terms,  i.e.,  Kn,  K22>  K33,  K44,  K55,  The  stiffness  re¬ 

duction  factor  for  each  diagonal  tem  applies  for  all  the  K  values  in  the 
respective  column  (the  stiffness  reduction  factor  for  K33  =  the  stiffness 
reduction  factor  for  K33,  K^,  Kj^»  K—,  Kg-). 

Fuselage 

The  changes  in  the  fuselage  external  spring  load-deflection  r *rves  are 
shown  in  Figures  12  and  13.  The  parameters  that  are  varied  include  the 
load  level,  stroke  and  the  rate  of  loading.  Figures  12  and  13  represent 
the  load-deflection  curves  for  masses  10  and  16  (forward);  the  load  level 
for  the  external  springs  from  masses  11  (mid),  12,  13,  22,  and  23  (sides) 
are  equal  to  half  the  load  level  of  masses  10  and  16. 

Crushing  load- stroke  characteristics  are  given  in  Table  X. 

RESULTS 

Load  sensitivity  studies  using  computer  program  KRASH  and  the  existing  31 
lumped  mass  model  are  shown  in  Figures  l4  through  17  for  the  load  and  energy 
absorbed  versus  deflection  for  the  engine  (axial),  engine  (bending),  trans¬ 
mission  (axial),  and  landing  gear  (bending),  respectively.  The  results, 
shown  in  Tables  XI  through  XV,  compare  the  peak  load  and  peak  responses  for 
the  base  case  versus  the  various  load-deflection  changes  for  the  structural 
elements . 

Engine 

The  results  of  the  engine  load- sensitivity  runs  are  shown  in  Figures  l4  and 
15  and  Tables  XI  and  XII.  For  large  changes  in  axial  load- stroke,  the 
response  is  fairly  insensitive.  For  example,  case  El  versus  the  base  case 
shows  a  potential  inaccuracy  of  12%  if  the  shape  of  the  buckled  load- stroke 
curve  is  not  represented  exactly.  If  the  load-stroke  curve  is  considered 
flat  (E2),  then  a  20%  error  would  be  introduced. 


DEFLECTION,  IN. 

Figure  12.  Fuselage  External  Crushing  Load-Deflection  Curves. 


TABLE  X.  FUSELAGE  EXTERNAL  CRUSHING  LOAD- STROKE  CHARACTERISTICS 


Case 

Load  Level  * 

(M>) 

Crushing  Distance* 
(In.) 

Bottoming  Stiffness* 
(Lb/in.) 

BASIC 

1.65  x  10* 

3.5 

li 

3.3  x  10* 

FI 

I.65  x  10^ 

7 

3.3  x  10* 

F2 

k 

3.3  x  10 

7 

3.3  x  10* 

F3 

1.65  to  3.3  x  101* 

7 

3.3  x  10* 

Fk 

3.3  x  102* 

3.5 

3.3  x  104 

F5 

I.65  x  101* 

3.5 

6.6  x  10^ 

F6 

:..64  x  iok 

3.5 

1.65  x  101* 

F7 

k 

3.3  X  10* 

7 

6.6  x  10^ 

f8 

k 

3.3  X  10 

7 

1.65  x  10 

*  Values  shown  are  for  masses  10  (aft)  and  1 6  (fwd). 

All  other  external  springs  are  changed  in  the  same  proportion. 
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DEFLECTION,  IN. 

(b)  ENERGY  ABSORBED  VERSUS  DEFLECTION 

Figure  15.  Engine-Mount  Bending  Load  and  Energy- 
Absorbed  Versus  Deflection  Curves. 
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TABLE  XI.  ENGINE  MOUNT  AXIAL  LOAD  SENSITIVITY  STUDY  RESULTS 


Condition 

Peak  Load  (lb) 
g  Time  A 

max  • 

Peak  Response  (G’s) 
g  Time  A 

max  • 

Percent  Change 

Load 

Response 

Base 

500  * 

49 

- 

- 

El 

2,500  * 

55 

400 

12 

F-2 

7,200 

6o 

1400 

20 

E3 

12,000 

72 

2400 

44 

*  Peak  Load  =  7200  lb  (Figure  14)  at  A  =  .15"  to  .2".  However  peak 
response  occurs  later  in  run  when  axial  load  drops  off  (  A  =  .6")  and 
is  influenced  by  bending  forces. 


TABLE  XII.  ENGINE  MOUNT  BENDING  LOAD  SENSITIVITY  STUDY  RESULTS 

Peak  Load  (lb) 

Peak  Response  (G’ s) 

Percent  Change 

Condition 

Time  @  A  =  .38" 

Time  @  A=.38" 

Load 

Response 

Base 

26,000 

47 

- 

- 

EB1 

24,000 

46 

-7.6 

-2.1 

EB3 

25,000 

47 

-3.8 

0 

EB4 

30,000 

49 

+15. b 

+4.2 

L'B2 

32,000 

50 

+23. 

+6.4 

- 

Condition 

Peak  Load  ( lb ) 

Peak  Response  (G's' 

Percent  Change 

Time  @  A  =  .5" 

Time  @A=„5" 

Load 

Response 

Base 

26,000 

50 

- 

- 

EB1 

12,000 

30 

-54 

-4o 

EB3 

20,000 

40 

-23 

-20 

EB4 

36,000 

60 

+38 

+20 

EB2 

Wttz  \  = ======== 

45,000 

72 

+73 

+44 

miMiinr.wBm 


Transmission 

The  results  of  the  transmission-mount  axial  load- sensitivity  runs  are 
shown  in  Figure  l6  and  Table  XIII.  A  comparison  of  the  percentage  of 
change  in  load  and  corresponding  percentage  of  change  in  response  indicates 
that  the  latter  parameter  is  far  less  sensitive  than  the  former  parameter. 
Conditions  T2,  T3  and  T4  represent  situations  in  which  the  load-carrying 
capability  will  decrease.  While  they  vary  considerably  in  load  vs  deflec¬ 
tion  beyond  the  failure  point,  the  effect  on  the  peak  response  varies  from 
25$  to  36.5$,  as  compared  to  the  basic.  Conditions  T  and  T1  represent 
conditions  ift  which  the  rate  of  load  increases  after  failure  occurs.  Under 
these  conditions  the  response  varies  frcm  3.6$  to  17.5$  as  compared  to  the 
base  case,  which  has  the  steepest  load-def lection  curve. 

Landing  Gear 

The  results  of  the  landing  gear  bending  load- sensitivity  studies  are  shown 
in  Figure  17  and  Table  XIV.  The  engine  and  transmission  mounts  are  insen¬ 
sitive  to  these  changes  because  failure  of  the  struts  occurs  approximately 
30  milliseconds  after  ground  contact  while  the  engine  and  transmission 
loads  occur  between  60  and  70  milliseconds  after  impact.  The  DRI  peak 
values  which  occur  between  80  and  90  milliseconds  after  impact  show  a 
slight  sensitivity  (~ 10$)  to  changes  in  landing  gear  load-stroke  values 
within  the  constraints  of  this  type  of  gear. 

Fuselage 

The  results  of  the  fuselage  load- sensitivity  studies  are  shown  in  Table  XV. 
Increasing  the  load  stroke  frcm  3.5  inches  to  7  inches  affects  the  response 
of  the  transmission  and  DRI's  by  approximately  10$.  (The  DRI  model  is 
described  in  Reference  1.  Basically  the  DRI  is  a  quantitative  measurement 
of  the  occupant  response  which  ic  related  to  the  probability  of  a  spinal 
compression  injury,  A  DRI  value  of  19,  for  example,  would  relate  to  a  10$ 
probability  of  spinal  compression  injury.)  Increasing  the  load  level  by  a 
factor  of  two  during  the  flat  portion  of  the  stroke  has  the  further  effect 
of  reducing  peak  responses  of  the  transmission  (*5$)  and  DRI  ( *  5 $).  The 
most  significant  factor  affecting  the  response  in  the  UH-1H  model  is  the 
bottoming  spring..  The  engine  response  is  increased  23$,  and  the  transmis¬ 
sion  response  is  changed  27$  when  the  bottoming  spring  is  doubled  (case  F5). 
When  the  bottoming  spring  is  halved  (case  F6),  the  engine  response  is 
reduced  by  20$  and  the  transmission  response  by  30$.  The  DRI's  responses 
in  the  UH-1H  model  are  more  sensitive  to  changes  in  load  level  and  stroke. 
The  DRI  peak  responses  occur  later  after  impact  compared  to  the  engine  and 
transmission  responses  and,  consequently,  are  affected  by  the  amount  of 
energy  that  is  absorbed  by  the  crushing  of  the  fuselage. 

• 

Conversely,  the  engine  and  transmission  responses  are  affected  by  the  rate 
of  load  change  due  to  the  stiffness  of  the  fuselage.  The  external  fuselage 
springs  bottom  (reach  peak  deflection)  at  between  60  and  70  milliseconds 
after  impact,  which  coincides  with  the  time  at  which  peak  engine  and  trans¬ 
mission  responses  occur. 
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TABLE  HU.  TRANSMISSION-MOUNT  AXIAL  LOAD  SENSITIVITY  STUDY  RESULTS 


Condition 

Peak  Load  (lb) 

@  Time  <  A 

max. 

Peak  Response  (G*s) 

@  Time  <  A 

max 

Base 

150,000 

86 

T1 

120,000 

83 

T 

95,000 

71 

T2 

85,000 

64.5 

T3 

75,000 

60 

t4 

70,000 

54.6 

TABLE  XIV.  LANDING  GEAR  LOAD  SENSITIVITY  STUDY  RESULTS 


Condition 


Relative  Energy- 
Absorption  to 
Failure 


33,000  (13) 

1+0,000  (15) 
35,000  (10) 
1+0,000  (12) 


Max. Response  (G’s) 

Max 

DRI 

Engine 

Transm. 

Aft 

Fwd 

53 

86 

63 

60 

53 

87 

63 

60 

52 

81* 

59 

55 

51 

86 

58 

54 

NOTE: 


*  Deflection  in  Inches 


DEFLECTION,  IN. 

(a)  LOAD  VERSUS  DEFLECTION 


0  k  8  12  16 


DEFLECTION,  IN. 

(b)  ENERGY  ABSORBED  VERSUS  DEFLECTION 

Figure  17.  Landing  Gear  Load  and  Energy  Absorbed 
Versus  Deflection  Curves . 
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TABLE  XV.  FUSEIAGE  LOAD  SENSITIVITY  STUDY  RESULTS 


ime  in  Milliseconds 


SVBSMtemB  AMLKSIS 


OBJECTIVE 


The  purposes  of  the  analysis  are  to: 

•  Develop  simplified  techniques  to  obtain  approximate  load-deflection 
characteristics  of  failure  and  post-failure  regions  of  typical  air¬ 
craft  substructure. 

•  Develop  ways  of  incorporating  empirical  and  analytical  procedures 
to  correlate  with  experimental  data. 

•  Extend  the  procedures  to  apply  the  analytical  techniques  to  struc¬ 
tures  which  have  a  different  geometrical  configuration. 

•  Predict  the  crashworthiness  capability  of  substructure  segments  to 
be  used  in  vehicle  crashworthiness  design  by  the  use  of  the  techni¬ 
ques  and  design  procedures  to  be  developed  in  this  study. 

METHODS  OF  ANALYSIS 


The  model  specimen  used  (Figure  18) .which  represents  the  lower  fuselage,  is 
a  composite  built-up  substructure  made  of  stiffened  panels  with  longitu¬ 
dinal  reinforcement  by  main  beams  along  the  sides.  The  specimen  dimensions 
are  46  x  l8  x  6.125  inches.  The  drawing  for  specimens  1  and  2  is  presented 
in  Figure  19-  Complete  details  are  given  in  the  SUBSTRUCTURE  TEST  PROGRAM 
section  which  follows.  Drawings  for  all  12  specimens, which  are  analyzed  and 
tested  during  the  study  described  herein, are  presented  in  Volume  II.  The 
specimen  is  approximately  half  size  except  that  the  skin,  web  and  angle 
stiffeners  are  actual  thicknesses.  Under  vertical  impact  conditions  the 
loads  are  transmitted  from  the  fuselage  segment  to  the  transmission  housing 
via  the  supports  at  the  four  edges  noted  in  Figure  20. 

In  this  study,  identifying  dominant  failure  mechanisms  and  post-failure 
behavior  is  of  primary  concern.  Emphasis  is  placed  on  examining  the  appli¬ 
cability  of  various  existing  analytical  and  semierapirical  approaches. 

The  structural  elements  selected  in  this  program  are  stiffened  panels,  which 
are  subject  to  crushing  deformation  modes  (or  end  shortening),  starting  with 
the  initial  buckling  of  the  webbed  structure, followed  by  post -buckling 
behavior,  which  is  characterized  by  geometrical  nonlinearity  as  well  as 
material  nonlinearity  due  to  plastic  deformation.  As  the  loads  are  further 
increased,  an  ultimate  point  is  reached,  and  beyond  this  point,  as  the  dis¬ 
placement  is  further  increased,  the  load  reduces  drastically. 

For  most  plates  and  stiffened  panels  under  in-plane  compressive  loads, 
failure  is  characterized  by  one  of  the  following: 

•  For  long  plates  and  long  panels,  the  structure  fails  mainly  due  to 
column  instability  behavior. 
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Figure  18.  Model  Specimen  of  Lower  Fuselage  Substructure. 


•  For  short  panels,  crippling  is  the  dominant  failure  mechanism. 

•  For  intermediate  length,  failure  is  a  result  of  combined  crippling 
and  column  behavior  with  possible  torsional  effects  added. 

On  the  other  hand,  the  longitudinal  beams  (top  and  bottom  skins)  are  subject 
to  bending  and  possibly  stretching  if  the  ends  are  constrained.  These  beams 
and  skins,  with  stiffened  panels  and  bulkheads  placed  between  them,  are 
expected  to  give  rise  to  very  complicated  interactive  response. 

In  conjunction  with  the  analysis  described  in  this  section, available  methods 
were  explored  and  evaluated.  The  applicability  as  well  as  complexities  of 
these  methods  is  discussed.  In  addition,  simplified  methods  for  the  pre¬ 
sent  study  are  proposed  and  applied  to  the  model  specimen  structures. 

Accurate  large -deflect ion  failure  and  post-failure  analyses  of  a  general  stif¬ 
fened  panel  require  complex  mathematical  formulations  including  the  incorpor- 
ation  of  plasticity  theory.  The  behavior  of  the  flat  plates  within  the  initial 
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Figure  19 .  Drawing  for  Specimens  1  and  2. 
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post-failure  range  is  generally  available.  However,  the  methods  of  esti¬ 
mating  maximum  strength  from  the  post-buckling  analysis,  References  87 
through  90,  are  limited  in  that  there  appears  to  be  a  lack  of  adequate  cor¬ 
relation  between  analytical  procedures  and  experimental  work.  In  addition, 
the  influence  of  parameters  such  as  material  and  fabrication-induced  imper¬ 
fections,  boundary  condition  influence,  and  geometry  changes  has  not  be«n 
thoroughly  examined.  Compared  to  analysis  and  testing  to  determine  failure 
points  of  structure  designs,  failure  beyond  the  ultimate  strength  point 
appears  to  have  received  very  little  attention.  Yet,  it  is  this  region, 
from  the  crashworthy  design  point  of  view,  that  is  vital  in  terms  of  energy¬ 
absorbing  capability  and  in  determining  whether  this  portion  of  substructure 
would  maintain  structural  integrity  throughout  the  crash  sequence. 

Except  for  the  cases  of  simple  geometry  and  boundary  conditions,  available 
methods  generally  require  a  fairly  sophisticated  level  of  computer  coding. 
Unless  the  use  of  sophisticated  programming  leads  to  reasonably  simple 
design  charts,  the  effort  would  not  meet  the  primary  objective  of  the  study, 
namely,  designer-oriented  tools.  Therefore,  the  development  of  such  com¬ 
puter  programs,  which  not  only  account  for  the  elastic  stability  behavior 
but  also  include  plasticity  phenomena  for  carrying  out  failure  and  post¬ 
failure  analysis  is  potentially  desirable.  However,  it  will  undoubtedly 
take  same  time  before  such  techniques  are  verified  as  reliable.  Therefore, 
it  is  important  to  develop  simplified  methods  which  fully  utilize  both 
empirical  and  analytical  approaches  for  predicting  large  deflection  behavior 
in  an  approximate  sense.  They  would  provide  not  only  some  insight  as  to  how 
one  can  get  load -deflect ion  behavior  for  similar  structures  without  going 
through  costly  test  and  computer  runs,  but  also  some  perspective  as  to  which 
way  to  pursue  this  relatively  unexplored  area. 

The  use  of  test  data  combined  with  analysis  bus  resulted  in  the  development 
of  semiempirical  approaches,  particularly  for  structure  subject  to  buckling. 
Buckling,  post -buckling  and  compression  strength  of  flat  plates,  composite 
elements  and  stiffened  panels  received  great  attention  during  the  last  four 
decades.  Although  substantial  theoretical  development  has  been  achieved  in 
recent  years  via  the  general  stability  theory  as  described  in  References  64 
and  65,  the  literature  does  not  indicate  that  there  has  been  any  additional 
extensive  experimental  investigation  since  the  comprehensive  undertaking  by 
NACA,  described  in  Reference  12,  in  the  late  1950' s. 

Although  the  literature  contains  comprehensive  information  regarding  buck¬ 
ling  of  plates  (via  design  charts)  obtained  by  the  classical  analysis,  and 
an  estimate  of  ultimate  strength  based  on  the  empirical  correlations  with 
experimental  test  data,  this  data  excludes  the  use  of  end -shortening  char¬ 
acteristics  and  post -failure  behavior.  Furthermore,  structural  imperfec¬ 
tions  dealing  with  rivet  strength  and  associated  failure  modes  (rivet  pitch, 
offset,  rivet  spacing)  as  well  as  geometrical  imperfections  were  rarely 
considered  in  modern  analytical  approaches  to  nonlinear  stru 'tural  analysis. 

References  87  and  89  describe  the  attempts  that  have  been  made  to  predict 
collapse  loads  by  their  analyses  for  flat  plates  under  in-plane  compression.- 
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Recently,  in  Reference  90  it  was  suggested  that  by  extending  Koiter's 
theory  (Reference  64)  and  combining  the  elastic-perfectly  plastic  theory, 
a  method  to  estimate  the  collapse  loads  for  the  lateral  buckling  of  a  rec¬ 
tangular  beam  can  be  achieved.  Although  some  of  these  approaches  must 
prove  to  be  feasible,  the  literature  indicates  that  the  current  state  of 
the  art  is  such  that  no  post-failure  analysis  of  composite  structures,  made 
of  stiffened  panels  and  beams,  exists  for  typical  fuselage  substructure. 

Therefore,  in  the  following  sections  the  emphasis  has  been  placed  on  using 
semiempirical/analytical  methods  to  develop  simplified  procedures  for  pre¬ 
dicting  crashworthy  (or  post-failure)  load-deflection  characteristics  of 
typical  built-up  structure  used  in  fuselage  design. 

PROCEDURE 


The  general  procedure  for  determining  the  total  load-deflection  character¬ 
istics,  including  failure  and  post -failure  behavior,  for  the  composite 
structure  selected  for  the  test  program  is  described  in  the  following  sec¬ 
tions.  Included  in  the  procedure  is  a  step-by-step  process  with  the  fol¬ 
lowing  sequence: 

•  Prediction  of  failure  loads  for  stiffened  panels 

•  Post-failure  analysis  of  stiffened  panel 

•  Main  beam  and  bottom  skin  analysis 

•  Total  load -deflect ion  curve 
Prediction  of  Failure  Loads  for  Stiffened  Panels 

The  various  failure  modes  of  general  stiffened  (riveted)  panels  are  shown 
in  Figure  21  (Reference  12).  Depending  upon  the  rivet  pitch,  spacing,  and 
strength,  a  ^iven  structure  may  fail  either  by  interrivet  failure,  wrink¬ 
ling  failure,  or  monolithic  behavior. 

Using  the  semiempirical/analytical  methods  described  in  References  12,  76, 
and  83  for  stiffened  short  panel  elements,  the  failure  loads  are  estimated. 
In  addition,  the  analysis  outlined  is  based  on  the  following  assumptions: 

•  At  the  threshold  of  failure  loads,  full  plastic  hinges  are  developed 
at  constrained  supports  and  mid  cross  section  of  the  panels. 

•  The  free  warping  energy  of  the  flange  of  the  stiffening  angles  is 
neglected. 

9  The  effect  of  strain  hardening  is  neglected. 

•  The  influence  of  the  axial  force  on  the  plastic  hinge  mechanisms  of 
the  stiffened  panels  is  neglected. 
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RIVET  SPACING 

Figure  21.  Various  Failure  Modes  of  Short  Riveted  Panels  (Reference  12). 


•  The  local  failure  pattern  changes  during  the  post-failure  stages  are 
neglected. 

•  The  effect  of  geometrical  imperfection  sensitivity  is  neglected. 

Since  the  present  stiffeners  are  of  one-corner  angles,  the  suggested  fail¬ 
ure  stress  formula  by  Gerard  in  Reference  12  is: 


t  t  /  — p  \  l/2 

m 

Of  =  0f 

w  s  MjEA 

^cy 

A  \  aCy  / 

Correlating  test  data  from  Reference  92  and  casting  them  into  the  Schuette's 
method  (Reference  93), the  constants  and  m  for  equal  angles  are  given  in 
Reference  12  to  be: 


sf 


m 


0.665 

O.85 
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(2) 

(3) 


Wrinkling  failure  is  expressed  in  Reference  9^  as  shown  below  in  equa¬ 
tion  (It): 


K  .  it 
w 


T)  •  T|  •  E 


aw  = 


(*)' 


(4) 


12  (1  -  v  ) 

where  and  t]  can  be  estimated,  respectively,  from  References  93  and  94. 

The  empirical  criterion  for  at  least  wrinkling  failure  rather  than  inter¬ 
rivet  failure  given  in  Reference  12  is  shown  to  be 


~<  1.27  K  1,2  (5) 

b  w 

s 

Instead  of  using  Equation  (4),  which  requires  Kw  to  determine  CTw»  a  simpli¬ 
fied  approximation  is  given  in  Reference  92  as 

5w-.«H|iil(ij  ’(£)(£)'*  (6) 

Another  approach  to  determine  the  effects  of  rivet  fastening  on  the  failure 
stress  based  on  the  monolithic  structural  configuration  is  shown  empiri¬ 
cally  in  References  96  and  97*  The  empirical  data  is  reproduced  from  Refer 
ence  12  in  Figures  22  and  23. 

When  the  sheet  material  is  different  from  the  stiffener  material,  the  effec 
tive  yield  stress  can  be  approximated  as 


a 

cy 


+  a 

cyw 


-  1) 


s 


(7) 


where 


=  £ 


(8) 


The  method  of  determining  the  number  of  flanges  to  be  used  in  Equation  (l) 
is  shown  in  Figure  24. 

Equations  (l)  through  (6)  provide  sufficient  information  for  the  present 
study  to  determine  the  failure  stress  for  a  given  stiffened  panel  segment 
under  inplane  compressive  loads. 


Figure  22.  Effect  of  Rivet  Tensile  Strength,  Pitch,  and  Diameter  Upon 
Compressive  Strength  of  Short,  Riveted,  Aluminum -Alloy 
Z-Panels,  2117-TU  Rivets  (Reference  12). 


Figure  23.  Effect  of  Rivet  Tensile  Strength,  Pitch,  and  Diameter 
Upon  Compressive  Strength  of  Short,  Riveted,  Aluminum- 
Alloy  Z-Panel,  Various  Rivet  Materials  (Reference  12). 
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(a)  z  STIFFENED  PANEL 


f  =  6  =  NUMBER  OF  FLANGES 
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(b)  HAT  STIFFENED  PANEL 


(c)  ONE-ANGLE  STIFFENED  PANEL 


Figure  2h .  Method  of  Determining  Number  of  Flanges. 
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In  order  to  determine  the  corresponding  failure  load  from  the  failure 
stress,  the  so-called  "effective  width  concept"  is  introduced.  Based  on 
the  work  described  in  References  98  and  99  >  an  extensive  test  made  by  C.I.T. 
during  the  early  1930's  suggested  the  following  formula  which  is  obtained 
from  Reference  83: 


where  n  =  0.37 


(£■) 


for  the  sheet  effective  width. 


(9) 

(10) 


The  method  of  determining  the  effective  width  for  various  configurations  is 
shown  in  Figure  25,  adopted  from  Reference  83. 


Once  acr  from  the  classical  buckling  theory  is  determined  as  described  in 

Reference  100,  the  effective  area  of  the  sheet  can  be  calculated.  Then  the 
failure  load  can  be  obtained  by 


b  ■  <s  +  E  V  °f 


(11) 


Post-Failure  Analysis  of  Stiffened  Panels 

A  representative  schematic  sketch  of  a  typical  segment  of  stiffened  panels 
during  the  post-failure  deformation  is  shown  in  Figure  26. 


From  the  energy  consideration, 

J  P  .  6u  =  2 


[/ ve+/ 
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From  the  geometry, 


tan  ^  Z 


(12) 


(13) 


where 


J  iu  -  u2A 


69  = 


ae 

az 


az 

au 


6u 


(14) 

(15) 


Substituting  Equation  (15)  into  Equation  (12)  and  applying  the  variational 
method,  Equation  (12)  can  be  shown  as 


/I 
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dZ 

"Su 


55 


6u  =  0 


(16) 


(a)  HINGE  MECHANISM  OF  POST-FAILURE  PATTERN  OF  A  STIFFENED  PANEL. 


VIEW  A-A 

(b)  INTERMEDIATE  WARPING  OF  FLANGE  DUE  TO  FORCED  BENDING  OF  WEB. 


(c)  LIMITING  CASE  OF  FLANGE  FOLDING  DUE  TO  FREE  WARPING. 


(d)  LIMITED  CONSERVATIVE  CROSS  SECTION  FOR  OBTAINING  PLASTIC 
BENDING  MOMENT. 

Figure  26.  Plastic  Hinge  Mechanisms  and  Effective  Cross  Section  To 
Determine  Plastic  Moment. 
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From  Equation  (l6),  after  sane  mathematical  manipulation 


_5 e _ 

\/fu  -  u2/l» 

where 


Mp  =  Mpx  +  ^ 


(17) 


(18) 


Calculation  of  Plastic  Hinge  Moment 


It  has  been  shown  in  Reference  75  that,  after  buckling,  the  bending  stress 
in  the  plate  increases  in  proportion  to  |  .  (a  -  aQT)  j  H  ,  whereas  the 

direct  compressive  stress  decreases  in  proportion  to  (<y  -  a^T  ).  The  pro- 

irtionality  factor  depends  on  the  boundary  conditions,  aspect  ratio  of  the 
plate,  and  buckled  wavelength. 


In  particular,  for  an  aspect  ratio  of  a./b  =  1,  and  with  the  plate  under 
uniaxial  compression  and  free  to  expand  laterally,  the  bending  stress  can 
be  shown  as 


(19) 


Therefore  if 


Vv  ■ 6 

then  the  bending  stress  at  the  outermost  fiber  of  the  sheet  skin  becomes 

<rb  =  t  0.88  CTf  (19a) 


and  if 


—  =  0.25 
°cy 


then 


01  =  t  0.22  O 

b  cy 


(20) 


Under  this  circumstance,  assuming  linear  bending  distribution  across  the 
cross  section,  the  yielding  will  set  in  on  the  stiffener  at  the  location 
approximately  2.3  times  the  sheet  skir,  thickness  from  the  neutral  axis  of 
the  sheet  skin. 


Although  it  is  premature  to  assume  that  at  collapse  loads  the  plastic  hinges 
have  formed  as  shown  in  Figure  26,  from  the  foregoing  assessment  it  can  be 
argued  that  this  assumption  is  not  too  unreasonable. 

The  effect  of  axial  force  on  the  plastic  hinge  mechanism  can  be  obtained  by 
applying  the  interaction  mechanism  described  in  Reference  101  and  shown  in 
Figure  27* 

Then 


It  can  be  shown  that  at  the  failure  load 


Figure  27*  Approximate  Yield  Function  for  Combined  Bending 
and  Axial  Force  (Reference  101) . 


Equation  (23)  implies  thp*  the  neglect  of  axial  force  would  give  about  a 
6-percenb  overestimate  of  the  plastic  moment  at  the  highest  direct  com¬ 
pressive  stress  state. 

As  will  be  discussed  below,  the  neglect  of  axial  force  and  the  assumption 
of  full  plastic  hinge  formation  at  the  threshold  of  maximum  load  are  criti¬ 
cal  when  compared  with  how  the  M  value  is  obtained  due  to  the  free  warping 
flange  portions  of  the  stiffened*paneis. 

An  assessment  of  the  effect  of  frge  flange  warping  instability  is  made  qua¬ 
litatively  in  Reference  102  using  Vlasov's  theory.  It  is  concluded  from  the 
assessment  that  the  effect  of  warping  on  the  simplified  plastic  failure 
mechanism  can  be  neglected  without  degrading  the  accuracy  of  the  post- 
failure  analysis. 

Therefore,  the  conservative  full  plastic  moment  of  the  stiffened  panel  may 
be  when  the  flange  portion  is  assumed  to  have  folded  completely  against  the 
web,  so  the  effective  conservative  cross  section  thickness  of  the  stiffener 
would  be  twice  the  web  thickness  as  shown  in  Figure  26.  However,  an  accu¬ 
rate  assessment  of  the  upper  bound  plastic  moment  is  difficult  due  to  the 
complex  behavior  of  the  free  warping  flange  during  the  post -failure  stage. 

Main  Beam  and  Bottom  Skin  Analysis 

The  subelement  of  bottom  main  beams  and  skin  is  shown  in  Figure  28(a).  In 
order  to  simplify  the  subsequent  analysis,  the  following  assumptions  are 
introduced: 

•  The  subelement  is  modeled  as  a  single  beam  with  simple  supports  at 
bulkhead  stations  A,  B,  E,  and  F,  respectively,  as  shown  in  Figure 
28(b). 

•  In  view  of  the  solid  support  system  between  stations  C  and  D,  it 
is  assumed  that  the  bending  rigidity  for  this  portion  of  beam  has 
an  infinite  value,  i.e.,  I^/lg  =  0. 

•  Although  top  beams  at  station  A  and  F  are  expected  to  undergo  a  cer¬ 
tain  bending,  and  some  shear  flow  for  the  corresponding  side  stif¬ 
feners,  it  is  assumed  that  this  bending  is  not  too  influential. 
Therefore,  the  surface  of  the  top  skin  and  beams  throughout  the 
entire  deformation  is  assumed  flat. 

If  the  above  assumptions  are  introduced,  the  beam  deflection  of  Figure 
28(b)  at  C  and  D  will  be  the  same  as  the  shortening  distance  of  the  panels 
between  stations  C  and  D. 

Following  the  procedures  described  in  References  103  and  104,  it  is  found 
that  the  plastic  hinge  forms  at  stations  C  and  D  simultaneously,  and  tie 
corresponding  applied  force  and  deflection  are 
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1  30E  Ix  P 

As  load  is  further  increased  the  bending  moments  at  B  and  E  reach  the 
plastic  moment  and  the  corresponding  force  and  deflection  are 

4  .  M 


Mp  .  r 

2  =  2EI^  (27) 

Beyond  this  stage,  the  load  does  not  increase,  and  the  deflection  increases 
until  the  hinge  rotation  reaches  complete  failure  at  the  location. 

H  e  foregoing  analysis  has  some  theoretical  drawbacks  when  the  displacement 
becomes  large,  due  to  the  introduction  of  the  simply  supported  assumption. 

For  example,  when  the  crushing  distance  becomes  about  3.0  in.,  assuming  that 
the  axial  strain  is  uniformly  spread  throughout  the  beam,  the  average  axial 
strain  is  about  0.03  in. /in.  This  is  certainly  more  than  failure  strain. 
However,  if  roller  guide  supports  rather  than  simple  supports  are  remodeled 
for  the  main  beam,  then  the  axial  influence  can  be  eliminated.  The  assump¬ 
tion  would  compensate  for  bending  of  the  cantilevered  portion  of  the  speci¬ 
men.  This  arbitrary  drawback  seems  to  be  contradictory  in  a  strict  theore¬ 
tical  sense ;  however,  the  approximate  load -deflect ion  curve  obtained  by 
this  assumption  is  believed  to  represent  the  dominating  deformation  mode. 

Total  Load  Deflection 


In  the  present  model  substructure,  the  stiffened  panels  are  assembled  to 
form  rectangular  cross  sections.  At  each  edge,  at  .least  one  or  more  stif¬ 
feners  join  the  panels  by  rivets,  together  with  sheet  skins. 

In  order  to  obtain  a  load  deflection  curve  of  this  typical  cellular  cross 
section,  two  factors  must  be  considered: 

•  The  effect  of  the  jointed  edge  constraint  on  the  initial  local 
buckling  of  the  composite  cellular  cross  sections. 

•  The  effect  of  the  constraint  of  the  joining  stiffening  members  at 
the  corners  during  the  failure  and  post-failure  stages. 
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Neglect  of  the  first  effect  can  be  alleviated  by  the  fact  that  the  so-called 
"four  moment  equations",  for  buckling  of  rectangular  tube  sections,  would 
render  no  more  reliable  results  than  the  separate  buckling  estimate  of  the 
individual  panels  when  the  tube  itself  has  stiffening  members  as  in  the  pre¬ 
sent  case. 

The  effect  of  corner  constraint  on  the  failure  load  presents  greater  com¬ 
plications,  and  this  will  require  indeed  a  fairly  involved  analysis.  If 
a  conservative  failure  load  is  emphasized,  it  is  assumed  that  the  connecting 
corner  stiffener  belongs  either  to  the  side  stiffener  or  the  bulkhead; 
and  neglecting  the  interaction  between ,two  neighboring  stiffened  panels , 
the  approximate  failure  load  can  be  estimated  from  the  individual  failure 
load  of  stiffened  panels. 

Another  complication  arises  for  the  nonuniformly  compressed  portion  of  the 
side  skin  between  stations  B  and  C,  and  D  and  E.  These  plates  may  be 
approximated  as  subjected  to  a  linearly  varying  edge  compressive  load  as 
described  in  Reference  88. 

While  the  involved  analysis  to  determine  the  failure  load  for  this  portion 
is  desirable,  in  view  of  the  effective  width  and  average  maximum  stress  con¬ 
cept  shown  in  Figure  25  plus  the  assumption  that  the  bulkhead  stationed  at 
B  and  E  would  remain  within  the  elastic  range,  another  approximation  is 
introduced  which  assumes  that  only  one-half  of  the  skin  width  participates 
in  failure  load  and  post-failure  deformation.  It  is  noted,  however,  that 
no  quantitative  assessment  was  sought  to  evaluate  this  approximation. 
Nevertheless,  preliminary  analysis  vs.  experimental  observation  tends  to 
justify  this  approximation. 

Column  Effect  on  the  Failure  Strength  of  Stiffened  Panels 

In  the  preceding  section,  the  crushing  stress  (or  failure  stress)  formulas 
are  based  on  the  value  of  L'/o  =  20.  For  a  sufficiently  larger  value  of 
L'/P  the  panels  behave  as  columns,  and  the  ultimate  stress  may  oe  calcu¬ 
lated  by  the  well-known  Euler's  formula 

<r  =  tt2  E/(L'/P)2  (28) 

However,  for  the  transition  range  as  shown  in  Figure  29,  the  parabolic 
approximation  based  on  the  value  of  cr  at  L'/p  =  20  can  be  obtained  from 
the  following  equation  derived  and  presented  in  Reference  12. 
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Figure  29.  Critical  Stress  Versus  L'/P  Ratio. 

For  aluminum  material,  the  stress  at  L'/P  =  20  is 

a20  =  ^  EA2°)2  =  ^/in.2  (30) 

To  determine  (L'/P)»  first  must  be  obtained  from  Equation  (29).  In  cal¬ 
culating  ,  first  P  and  L'  must  be  determined.  These  values  are  obtained 
following  procedures  similar  to  those  described  in  Reference  83. 

Figure  30  shows  the  scatter  values  of  end  fixity  obtained  from  testing  flat- 
ended  panels.  It  is  observed  in  Reference  83,  from  the  test  data,  that 
about  75  percent  of  the  test  points  lie  within  t  5  percent  of  the  curve  for 
C  =  3*0.  Therefore,  it  is  suggested  that  a  value  of  C  =  3*0  be  used  as  a 
good  average  value  of  the  panel  end  fixity.  The  same  value  will  be  adopted 
in  this  study  since  there  is  limited  data  available.  Using  C  =  3.0,  the 
following  expression  for  L'  is  obtained: 

L'  =  L/  n/c"  =  .  57735L  (31) 

In  order  to  determine  P,  the  necesrary  steps  are: 

1)  Calculate  b  as  shown  in  Figure  31  based  on  the  failure  stress 

e 

af. 

2)  Calculate  p  for  stiffeners  alone. 

w 

3)  Calculate  the  distance,  S,  from  the  mid  cross  section  of  the  skin 
to  neutral  axis  of  the  sbiffener  as  shown  in  Figure  31* 

4)  Use  Figure  32  or  obtain  from  Equation  (32): 
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Figure  30.  Failure  Stress  Versus  Slenderness  Ratio  of  Stiffened  Panels. 


This  completes  all  the  necessary  informati.on  to  obtain  a  new  cr  .  The  next 
steps  to  obtain  the  ultimate  load  are: 

•  Calculate  new  we/b  based  on  a  new  o  (see  Equation  9  and  Figure 

25).  S  C0 

•  Calculate  £A  based  on  the  new  value  of  b  . 

se  e 

Then  the  ultimate  load  is  calculated  as  previously  by  Equation  (ll). 


General  Application  of  the  Procedure 


The  procedure  presented  in  the  previous  sections  to  determine  the  failure 
stress  for  the  angle-type  stiffened  panel  is  also  applicable  to  a  variety 
of  other  types  of  panels,  if  the  parameters  such  as  P,  n  and  m  (in  Equa¬ 
tion  34)  are  adjusted  to  the  experimentally  evaluated  data. 

In  general,  according  to  the  data  in  Reference  12,  the  failure  strength  can 
be  generalized  by  the  following  form 


—  =  e 

a  g 
cy 


g  .  t  .  t 
B  w  s 


(ill 


(34) 


where 


g  -  the  number  of  flanges  +  the  number  of  cuts 


m 

S  (-2) 
c  V 

for  multi-corner  elements 

(35) 

f  m 

o  (£) 
f  V 

for  one-corner  elements 

(36) 

8  (!) 
g 

for  two-corner  elements  or  plates 

(37) 

In  Figure  33,  representative  types  of  stiffeners  with  from  one  to  six  cor¬ 
ners  are  shown.  In  Figure  34  and  Table  XVI,  values  of  B  ,  g,  n  and  m, 
obtained  from  Reference  12,  for  different  configurationssare  listed.  A 
detailed  illustrative  sample  problem  using  the  analytical  procedure  is  pre¬ 
sented  in  the  Design  Procedure  section.  The  following  subsection  presents 
analytically  predicted  subelement  and  total  load-deflection  curves  and  dis¬ 
cusses  the  manner  in  which  the  curves  are  obtained. 


Predicted  Subelement  and  Total  Load-Deflection  Curves 


Following  the  procedure  outlined  in  the  preceding  portion  of  the  Substruc¬ 
ture  Analysis  section,  the  failure  load  for  the  stiffened  panels  (and  stif¬ 
fener  and  bulkhead)  aro  determln*'!.  Next,  the  post-failuve  load-deflection 
curve  shape  is  obtained  for  each  of  the  stiffened  panels.  The  failure  de¬ 
flection  is  then  derived  from  the  intersection  of  the  ultimate  failure  load 
(failure  analysis)  and  a  coinciding  load  from  the  post- failure  analysis. 

The  assumptions  regarding  the  plastic  hinge  formation  at  the  threshold  of 
failure  implies  that  the  yield  strain  («y)  has  been  reached  which  can  be 
considered  as  the  upper  bound  strain  at  the  failure  load.  Thus,  the 
derivation  of  the  deflection  from  the  failure  point,  although  approximate, 
is  in  effect  based  on  the  upperbound  value  of  the  yield  strain.  The 
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(a)  ONE-CORNER  ELEMENT 


(b)  TWO-CORNER  ELEMENT 


(c)  THREE-CORNER  ELEMENT 


□  T.  Jl_ 


(d)  FOUR-CORNER  ELEMENT 


(f)  SIX-CORNER  ELEMENT 


Figure  33.  Stiffeners  With  Various  Corners. 


ANCLE 


PLATE 


TUBE 


(a)  EXTRUDED  ANGLE  TYPE 


T-SECTION 


CZZZZZZTZZZZZ3 


H-SECTION 
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(c)  FORMED  ANGLE  STIFFENERS 


(d)  FORMED  HAT  STIFFENED  PANEL  (e)  FORMED  Y  STIFFENED  PANEL 
Figure  34.  Various  Stiffener  Configurations. 
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SUBSTRUCTURE  TEST  PROGRAM 


TEST  PROGRAM  OBJECTIVES 

Bie  objectives  of  the  test  program  are  to: 

(1)  Develop  testing  techniques  for  determining  load-deflection 
characteristics  of  structural  elements  subjected  to  loads 
involving  crash-oriented  large  def orraations . 

(2)  Establish  the  validity  of  simplified  analytical  techniques 
for  use  in  predicting  load-deflection  characteristics  for 
crashworthiness  analysis. 

(3)  Verify  the  application  of  simplified  analytical  techniques 
for  different  types  of  structural  elements  used  in  aircraft 
design. 

(4)  Obtain  information  which  can  be  used  in  future  studies  for 
improving  crashworthiness  analytical  capabilities. 

(5)  Determine  limitations  associated  with  simplified  analytical 
techniques  in  the  prediction  of  crashworthiness  capability. 

GENERAL 

It  is  important,  in  the  establishment  of  the  test  program,  that 
consideration  be  given  to: 

•  representing  typical  structured  designs 

•  ascertaining  the  effects  of  detail  design  differences  on  load- 
deflection  characteristics 

•  representing  a  crash  loading  condition 

•  relating  the  test  data  to  procedures  for  predicting  load- 
deflection  characteristics 

•  relating  the  test  data  to  analytical  modeling  data  require¬ 
ments. 

Twelve  tests  were  performed  to  aid  in  evaluating  whether  or  not  simplified 
analytical  techniques  and  experimental  tests  provide  data  of  sufficient 
accuracy,  when  applied  in  a  crashworthy  analysis  of  helicopters,  to  yield 
reliable  results. 

In  examining  the  type  of  structure  used  in  helicopter  design,  the  most 
critical  area  from  the  standpoint  of  energy  absorption  and  habitable 
volume  is  the  fuselage.  In  addition,  this  area,  because  of  its  built-up 
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construction,  is  most  difficult  to  analyze  for  large  deformation.  Struc¬ 
tural  elements  representative  of  other  regions  of  the  aircraft,  i.e., 
landing  gear,  engine  mount,  and  tail  boom ,are  of  lesser  complexity  than 
the  fuselage  structure,  and  in  many  instances  they  are  not  primary  energy 
absorbers.  In  addition,  there  is  more  test  and  analytical  data  available 
for  these  types  of  structural  elements  than  for  multi -veb  construction. 

Testing  of  several  different  types  of  representative  fuselage  substructure 
provides  the  opportunity  to  obtain  data  which  produces  meaningful  load- 
deflection  characteristics  such  as  failure  modes  due  to  the  design  changes 
associated  with  the  various  elements  within  the  substructure. 

The  tests  are  to:  (l)  ascertain  if  the  analytical  techniques  to  be 
developed  in  this  study  for  determining  approximation  of  large  deformation 
response  characteristics  are  valid,  and  (2)  obtain  experimental  load- 
deflection  data. 

Inasmuch  as  an  infinite  number  of  combinations  of  loading  conditions  and 
edge  restraints  are  available,  care  is  exercised  to  select  only  those  com¬ 
binations  that  will  provide  the  necessary  data  to  fulfill  these  require¬ 
ments.  Once  analytical  techniques  are  compared  with  experiments  and  veri¬ 
fied  regarding  their  limitations,  they  can  be  applied  to  variations  in  de¬ 
sign,  geometry,  and  constraints  with  a  greater  degree  of  confidence. 

TEST  SPECIMENS 


A  representative  pylon  support  structure  is  selected  as  a  basic  test 
specimen  from  which  normal  types  of  detail  design  variations  are  made 
to  alter  its  load -deflect ion  characteristics  under  large  deflection 
loading  conditions.  The  pylon  support  structure  consists  of  an  inner 
and  an  outer  skin,  representative  of  the  basic  floor  and  the  helicopter 
lower  fuselage  skin  respectively,  tied  together  by  a  series  of  beams 
running  fore  and  aft  and  bulkheads.  This  structure,  during  helicopter 
accidents,  is  the  means  by  which  the  vertical  ground  impact  is  transmitted 
to  the  transmission.  Accordingly,  its  modes  of  failure,  the  amount  of 
energy  absorbed,  and  the  pulse  shape  produced  are  of  interest  in  assessing 
the  crashworthiness  capabilities  of  a  helicopter. 

Figures  43  and  44  show  specimens  of  the  two  different  depths  that  are 
represented.  The  specimen  shown  in  Figure  43  is  6.125  inches  deep  and  is 
applicable  to  specimens  1  through  9,  Table  XVII.  The  specimen  shown  in 
Figure  44  is  12.125  inches  deep  and  represents  specimens  10  through  12. 

The  specimen  configurations  and  test  conditions  include: 

static  load 

dynamic  impact  (l4  ft/sec  <  velocity  <  30  ft/sec) 
four-edge  support 
two-line  support 

skin  web  thickness  from  .025  in.  to  .064  in. 
angle  stiffener  thickness  from  .032  in.  to  .072  in. 
specimen  weight  from  13.6  lb  to  31*25  lb 
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TABLE  XVII.  SPECIMEN  DESIGN  CONFIGURATION,  DIMENSIONS 
AND  TYPES  OF  TESTS  PERFORMED* 
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3m  Volune  II  Figures  2  through  6  xor  detail  drawings  of  Speclatana 
Inpect  hand  and  carriage  weight  •  8bO  lb  (Teat  5),  •  1069  lb  (Teat  8) 
bfact  head  and  carriage  weight  •  3550  lb  (Teats  W,  6),  -  3585  lb  (Teat  7) 
Six  load  cells  were  installed  below  the  specimen  support 
Una  supported,  all  other  teats  were  with  four-edge  supports 


number  of  angle  stiffeners  from  12  to  32 
specimen  depth  from  6.125  in. to  12.125  in. 

Each  test  specimen  is  46  inches  long,  including  five  bays,  and  is  18  inches 
wide  with  a  side  panel  on  each  side  to  represent  the  main  bulkhead.  The 
test  specimen  basic  overall  dimensions  and  bulkhead  and  stiffener  spacing 
are  approximately  half  size. 

The  typical  lower  fuselage  bulkhead  and  stiffener  arrangement  shown  pre¬ 
viously  in  the  Substructure  Analysis  section  (Figure  20)  is  presented 
again  in  Figure  45.  Illustrated  in  this  figure  is  the  location  of  the 
test  specimens  relative  to  the  entire  lower  fuselage,  the  landing  gears, 
the  transmission  pylon  and  the  pylon  support.  The  drawing  for  Specimens 
1  and  2  is  presented  in  Figure  46.  Drawings  for  all  12  specimens  fabri¬ 
cated  and  tested  during  the  program  are  provided  in  Volume  II.  Table  XVII 
shows  the  dimensions,  thickness,  number  of  angles,  number  of  stiffeners, 
nominal  weight,  design  configuration,  and  type  of  test  performed  for  each 
of  the  12  specimens. 

INSTALLATION 


The  planned  test  program  required  chat  both  static  and  dynamic  tests  be 
performed.  The  static  tests  are  performed  using  a  Baldwin  Universal 
testing  machine  (Figure  47).  The  static  setup  with  a  specimen  installed 
is  shown  in  Figure  48.  The  impact  head  is  used  throughout  all  the  tests 
— static  and  dynamic.  A  quasi-static  loading  procedure  is  followed  in 
which  a  controlled  constant  deflection  rate  is  applied.  Initially,  the 
rate  is  0.05  inch/minute.  After  the  initial  failure  load  is  reached,  the 
loading  rate  is  increased  to  at  least  2  inches/minute.  The  dynamic  tests 
are  performed  using  the  frame  structure  test  tower  and  carriage  shown  in 
Figure  49.  The  minimum  carriage  weight  and  aluminum  impact  head  is  840 
pounds.  The  drop  test  facility  has  sufficient  clearance  to  perform  30- 
ft/sec  impact  tests;  however,  the  weight  cf  the  carriage  and  impact  head 
and  the  energy-absorbing  capacity  of  the  structure  are  the  determining 
factors  in  the  final  impact  velocity  selected  for  a  particular  test.  The 
initial  drop  tests  (test  specimens  4,  5  and  6)  employed  a  four-edge  wooden 
support  between  the  specimen  and  the  ground.  The  load  was  obtained  from 
accelerometers  mounted  on  the  impact  head.  The  second  set  of  dynamic 
tests  (specimens  7  and  8)  was  performed  with  a  set  of  six  load  cells  in¬ 
stalled  between  the  specimen  and  the  ground.  This  setup  is  shown  in 
Figure  50  for  the  l4-ft/sec  and  27-ft/sec  impact  test  conditions.  In  the 
lower  velocity  test,  mass  was  added  to  the  impact  head  to  provide  compar¬ 
able  energy  levels  as  would  be  experienced  at  higher  impact  velocity  level. 
The  five  500-pound  weights  U3ed  can  be  seen  on  top  of  the  impact  head  in 
Figure  50a, 
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Figure  1*6.  Drawing  for  Specimens  1  and  2  (All  Dimensions  in  Inches 
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Figure  47.  Baldwin  Universal  Static  Testing  Machine. 


Figure  48.  Static  Test  Setup. 


preceding  page  blank 


85 


Best  Available  Copy 


w .  • 


it  .  t  -  •  ’  *<■»<•; 


INSTRUMENTATION 

For  the  initial  series  of  tests,  seven  strain  gages  were  installed  to 
measure  compression  and/or  tending.  The  locations  of  the  strain  gages 
are  described  in  Figure  51  and  Table  XVIII.  For  the  static  tests,  the 
load  is  measured  directly  with  a  force  transducer,  and  the  deflection  is 
obtained  using  a  linear  variable  differential  transducer  (LVDT) .  All  data 
is  recorded  directly  on  the  Lockheed  Central  Data  System  (CDS),  from  which 
response  time  history  and  force  acceleration  versus  deflection  plots  are 
obtained . 


In  the  dynamic  tests,  in  addition  to  the  strain  gage  and  displacement 
measurements,  three  vertically  oriented  accelerometers  (  >  100  g  range 
and  with  frequency  response  «  800  Hz)  are  used. 

For  the  first  set  of  dynamic  tests  (4,  5  and  6),  accelerometers  are  placed 
on  the  Impact  head  and  acceleration  data  (mass  x  g's)  is  used  to  determine 
the  load.  In  the  second  set  of  dynamic  tests  (7  and  8),  in  addition  to  the 
accelerometers  mounted  on  the  impact  head,  load  cells  are  installed  between 
the  test  specimen  support  and  the  ground. 

TEST  SEQUENCE 

The  purpose  of  the  initial  tests  is  to  verify  the  significance  of  inertia 
effects  and  design  configurations  (i.e.,  lightening  holes,  modified 
stiffener  arrangement)  on  the  load-deflection  characteristics  of  typical 
fuselage  structure.  Prior  to  the  actual  tests,  preliminary  tests  were 
performed  to  check  out  the  test  setup,  motion  of  the  carriage,  and  in¬ 
strumentation.  The  preliminary  tests  were  run  with  a  rubber  pad  in 
place  of  the  test  specimen  and  performed  for  a  drop  height  of  2.5  feet. 

The  results  of  the  preliminary  tests  indicated  that  the  drop  tower  and 
data  reduction  systems  were  satisfactory  fo<  performing  the  dynamic  tests. 

The  testing  sequence  used  is  as  follows: 


Static  test 
Static  test 
stiffeners) 
Static  test 
Dynamic  test 
Dynamic  test 
city) 

Dynamic  test 
Dynamic  test 
Dynamic  test 
specimen  7) 
Static  test 
Static  test 
Static  test 
specimen  10) 


(specimen  with  reduced  number  of  angle  stiffeners) 
(specimen  with  representative  number  of  angle 

(same  as  specimen  2  except  for  lightening  holes) 
(same  as  specimen  3>  impact  velocity  17  ft/sec) 
(specimen  2  stiffened  for  30  -  ft/sec  impact  velo- 

at  14  ft/sec(same  as  specimen  5) 

at  14  ft/sec  (load  cells  installed) 

at  27  ft/sec  (load  cells  installed  -  same  as 

(same  as  specimens  7  and  8) 

(l2.125- inch  depth) 

(l2.125-inch  depth;  skin  and  angles  thicker  than 
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DENOTES  STRAIN 
GAGE  LOCATION 


TABLE  XVIII.  STRAIN  GAGE  LOCATIONS 


Strain 
Gage  No. 


Location 


Web  of  center  bay,  end 
beam  inside 

Web  of  center  bay,  ?nd 
bean  outside 

I  Web  of  bulkhead  V 

Web  end  beam  at  junction 
of  bulkhead  "c" 

Web  of  end  beam  at 
junction  of  bulkhead  "d" 

Bottom  of  beam  flange 
2-1/2  inches  from  junc¬ 
tion  of  bulkhead  "d" 

Bottom  of  flange  2 
inches  from  junction  of 
bulkhead  "e" 


1 

Direction  Of 
Measurement 

Strain  Gage 
Mounted  On 
Test  Specimens 

Compression  and/or 
bending  of  end  beam 

1  thru  12  * 

Compression  and/or 
bending  of  end  beam 

1  thru  1? 

1  thru  0 

Axial  compression  of 
end  beam 

1  thru  3 

1 

Axial  compression  of 
end  beam 

1  thru  3 

Axial  strain  for  bend¬ 
ing  of  beam 

1  thru  3 

Axial  strain  for  bend¬ 
ing  of  beam 

1  thru  0 

On  specimen  9,  strain  gage  1A  is  placed  or-  .ne  free  flange 
of  the  outside  angle  stiffener  in  the  center  of  the 
middle  bay,  since  this  specimen  is  closed  on  the  bottom. 

**  Strain  gage  locations  are  shown  in  Figure  51. 


[*] 


12.  Static  test  {12. 125-inch  depth;  same  as  specimen  10  except  for 
line  supports) 

PHOTOGRAPHY 

During  the  static  tests,  time  lapse  (2  frames/sec  and  6  frames/sec)  photo¬ 
graphic  coverage  and  35  mm  slides  were  taken  during  the  application  of  load. 
In  addition,  still  photographs  of  test  setup,  test  specimen  and  post-test 
failure  modes  were  obtained. 

During  the  dynamic  tests,  high-speed  motion-picture  (500  frames/sec  and 
1000  frames/sec)  coverage  was  employed  in  synchronization  with  the  tape- 
recorded  data.  The  location  of  the  cameras  was  such  that  two  different 
angles  were  covered  on  each  test. 

RESULTS 

The  results  of  the  tests  are  shown  in  Table  XIX.  The  load-deflection 
curves  for  each  of  the  test  specimens  are  shown  in  Figures  52  through  63. 
Because  the  initial  dynamic  tests  were  performed  with  accelerometers  ■  count¬ 
ed  on  the  impact  head,  the  acceleration  time  histories  have  oscillati^rs 
from  which  it  is  difficult  to  accurately  interpret  the  proper  loads. 

Since  the  accelerometer  response  is  influenced  by  the  Impact  tower  vi  na¬ 
tion,  impact  head  and  carriage,  the  load-deflection  curves  for  specimens 
4,  5  and  6  (Figures  55,  56  and  57,  respectively)  have  oscillations  whjch 
are  believed  to  fluctuate  about  a  mean  value. 

Because  of  the  difficulty  in  interpreting  accelerometer  data,  tests  7,  8 
and  9  were  designed  to  test  the  same  configurations  (with  r,he  exception 

of  a  lower  skin  on  specimen  9)  under  the  following  three  equal  energy 

input  conditions: 

•  27  ft/sec  impact  velocity 

•  14  ft/sec  impact  velocity 

•  static  load 

The  results  of  these  three  tests  are  shown  in  Figures  58,  59  and  60, 
Specimens  10,  11  and  ]2  were  designed  to  be  approximately  twice  as  deep 
(12.25  inches)  as  the  other  specimens  (6.125  inches).  Specimens  10  and  12 
were  identical.  Specimen  11  was  designed  with  skin  and  angle  gage  thick¬ 
ness  twice  that  of  specimens  10  and  12. 

Specimens  10  and  11  were  cested  using  the  same  four-edge  metal  support 
chat  was  used  for  specimens  7,  8 -and  9.  Previously  a  four-edge  wooden 
block  support  was  used,  but  when  load  cells  were  added  to  the  dynamic  test 
setup,  the  stiffer  support  was  considered  to  be  more  desirable.  The  load- 
deflection  curves  are  shown  in  Figures  6l,  62  and  63.  Specimen  12  was 
tested  using  line  supports  beneath  the  side  panels  (Figure  64).  Specimen 
11  (138,000  lb)  was  considerably  stiffer  than  specimen  10  (39,000  lb),  and 
this  is  easily  observed  In  the  failure  loads  for  each  specimen. 
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Figure  54.  Lo&d  Versus  Deflection,  Test  3. 


Figure  55.  Load  Versus  Deflection,  Test  4. 
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Figure  64.  Setup  for  Specimen  12. 


As  anticipated, the  failure  load  of  Specimen  12  was  higher  them  that  experi¬ 
enced  by  specimen  10  (68,000  lb  vs.  39,000  lb). 

Figures  65  through  70  show  the  post-test  condition  for  specimen  1  through 
6,  respectively.  Figures  71  and  72  shows  a  side-by-side  comparison  of 
specimens  5  and  6  after  equal  energy  dynamic  impact  tests  at  30  ft/sec  and 
1.4  ft/sec,  respectively.  The  post-test  condition  is  very  similar  for  the 
two  specimens.  Figures  73  through  77  show  the  post-test  condition  of 
specimens  7;  8  and  9  individually  and  side-by-side.  The  test  of  these 
specimens  provides  data  regarding  the  comparative  effects  of  dynamic  versus 
static  testing  on  load-deflection  characteristics.  The  test  results  indi¬ 
cate  that  the  lower  the  impact  velocity,  the  more  severely  the  specimen  is 
deflected,  the  most  extreme  deflection  occurring  under  static  load.  In 
addition,  the  lower  the  impact  velocity,  the  less  springback*  is  evident 
immediately  after  the  load  is  released  from  the  specimen.  From  Table  XIX, 
the  springback  effect  is  shown  to  be  as  high  as  40$  of  the  maximum  deflec¬ 
tion  for  a  27-ft/sec  impact  velocity  (specimen  8)  but  is  reduced  to  2 % 
at  a  14-ft/sec  impact  velocity  (specimen  7)*  For  specimens  5  (30  ft/sec) 
and  6  (l4  ft/sec),  the  springback  effect  is  31$  and  12$,  respectively. 


.  ,  ,  .  ,  .  ,  /Maximum  deflection-permanent  deflection) 

*  Springback  is  defineo  hare  as  ( - wirimi-  - ) 
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Figure  71.  Comparison  of  Post-Test  Damage  "Or 
Specimens  5  and  6  (End  View). 


Figure  72.  Comparison  of  Post-Tes'  Damage  for 
Specimens  5  and  6  (Top  View). 
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Figure  .  Post-Test  Condition  of  Specimen  3. 
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Figures  78  through  82  show  the  post-test  condition  for  specimens  10,  11  and 
12  (12. 12 5- inch-deep  specimens).  Unlike  all  the  other  tests,  specimen  No. 
12  was  tested  with  two  line  supports.  The  post-test  condition  differs 
considerably  from  specimen  10  with  which  specimen  12  was  identical  except 
for  the  load  restraint  system  (four-edge  versus  two-line  support).  There 
is  more  of  a  beam  bending  mode  of  failure  present  in  specimen  12.  Speci¬ 
men  11  is  similar,  in  general  configuration,  to  specimen  10  except  with 
thicker  angles  and  skins  and  as  anticipated  withstands  a  considerably 
higher  failure  load  (  69,200  lb  versus  39.000  lb)  and  absorbs  more  energy 
(  73,800  in. -lb  versus  51>750  in. -lb).  Figure  8l,  in  which  a  comparison 
of  specimens  10  and  11  is  shown,  indicates  a  similar  type  of  failure  for 
both  specimens. 

A  review  of  post-test  damage  photographs  (Figures  65  through  80)  indicates 
that  the  assumptions  regarding  the  formation  of  a  plastic  hinge  at  the 
threshold  of  failure,  the  cross-section  of  the  structure  participating  in 
forming  a  plastic  hinge  (flange  warping),  and  the  symmetry  associated  with 
the  mode  of  failure  for  the  stiffened  panels  are  generally  valid.  How¬ 
ever,  the  failure  patterns  for  the  stiffened  panels  are  irregular  ^rom  one 
test  to  another  test  and  for  each  stiffener  location.  While  plastic 
hinges  do  form  at  many  locations  they  do  not  always  form  at  every  con¬ 
straint  support  or  necessarily  at  each  midpoint  of  all  the  stiffened 
panels.  There  appears  to  be  evidence  of  hinge  formation  in  varying  de¬ 
grees  for  specimens  2  through  11.  In  nearly  all  the  cases  the  angle 
stiffener  flanges  warp.  An  exception  to  this  is  specimen  1  which  had  the 
least  number  of  angle  stiffeners  and  failure  was  initiated  by  rivet  pull¬ 
out.  However,  in  about  half  the  test  results  there  is  noticeable  warping 
of  the  main  beam  flange.  For  several  of  the  specimens  (Figures  69-7*0 
warping  of  the  main  beam  flange  is  quite  noticeable.  The  failure  pattern 
of  the  stiffened  panels  indicates  symmetrical  failure  modes  for  nearly 
half  the  specimens  (l,  3,  7,  8,  .11)  and  an  equal  number  clearly 
unsymmetrical. 

The  effect  of  each  of  the  above-mentioned  analytical  assumptions  on  the 
comparison  between  the  analysis  and  test  results  is  discussed,  in  detail, 
in  the  following  section  entitled  Correlation.  Where  applicable,  poten¬ 
tial  modifications  to  the  analyses,  which  utilize  empirical  data,  that 
will  enhance  prediction  capability,  are  discussed. 

The  most  significant  implication  of  the  above  observations  is  that  the 
overall  failure  modes  of  the  specimens  are  predictable  although  the  exact 
failure  mode  of  each  individual  element  is  difficult  to  predict. 

Although  the  individual  stiffeners  ana  panels  deform  under  different  modes 
(due  to  differences  in  design  and  imperfections)  the  overall  deformation 
modes  indicate  a  consistent  trend.  For  example,  although  specimens  2  and 
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Figure  80,  Post-Test  Condition  of  Specimen  12 
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Figure  8l.  Comparison  of  Post-Test  Damage  for  Specimens 
10  and  11  (Top  View). 


Figure  82.  Comparison  of  Post-Test  Damage  for  Specimens 
10  and  12  (Side  View). 


3  are  supposedly  identical  in  design  and  energy-absorption  capability,  the 
post-test  configurations  show  entirely  different  deformation  inodes  for 
individual  elements,  but  the  load-deflection  curves  are  nearly  identical 
(Figures  53  and  5*0.  This  aspect  of  the  validity  of  the  analytical  assump¬ 
tions  relative  to  the  overall  results  is  discussed  in  more  detail  in  the 
Correlation  section  under  the  Remarks  on  the  Initial  Analytical  Assumptions 
subsection. 

All  the  strain  gage  data  that  was  obtained  during  the  test  program  is 
presented  in  Volume  II  under  the  section  entitled  Test  Data. 
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GENERAL 


A  comparison  of  test  and  analytical  data  for  the  12  test  specimens  is 
presented  in  Table  XX.  The  predicted  versus  test  load-deflection  curves 
are  presented  in  Figures  83  through  9^*  The  comparison  of  test  and  analy¬ 
tical  data  shovs  reasonable  agreement  for  the  crushing  characteristics  of 
the  model  specimens  selected  in  this  program  as  representative  of  a  seg¬ 
ment  of  the  UH-1H  fuselage. 

Test  specimens  1  and  2  afford  the  opportunity  to  assess  changes  in  design 
configuration  (difference  in  the  number  of  angle  stiffeners)  under  an 
identical  loading  condition.  Test  specimens  2  and  3  provide  an  opportunity 
to  compare  analysis  and  test  for  specimens  with  and  without  bulkhead 
lightening  holes,  under  identidaL  loading  conditions. 

A  comparison  of  load-deflection  characteristics  under  static  and  dynamic 
loading  conditions  is  obtained  from  the  analysis  and  test  of  specimens  3 
and  b  and  for  specimens  7>  8  and  9*  The  effect  of  the  impact  mass  on  the 
response  of  the  structure  can  be  obtained  from  the  analytical  and  experi¬ 
mental  results  between  specimens  5  and  6,  and  between  specimens  7  and  8. 

• 

Finally,  specimens  10,  11  and  12  provide  the  opportunity  to  exercise  the 
analytical  app.oach  for  deeper  (12.125  inches  versus  6.125  inches)  speci¬ 
mens  and  compare  the  effect  of  line  supports  in  lieu  of  four-edged 
supports . 

COMPARISON  OF  TEST  AND  ANALYSIS  LOAD-DEFLECTION  CURVES 

Results  of  the  analysis  and  tests  for  specimens  1  and  2  are  shown  in  Fig¬ 
ures  83  and  81*,  respectively.  The  analysis  predicts  the  peak  failure  load 
within  15 %.  The  deflection  at  peak  failure,  from  the  analysis,  is  a  lower 
value  (0.10  to  0.20  inch)  than  that  obtained  experimentally  (0.25  to 
0.1*0  inch).  The  deflection  at  the  analytically  predicted  peak  failure 
value  is  the  point  at  which  the  post-failure  load-deflection  curve  inter¬ 
sects  the  computed  peak- failure  load  value.  This  matching  process  is 
described,  in  detail,  in  the  sections  entitled  Substructure  Analysis  and 
Design  Procedures. 

The  lower  initial  slope  in  the  test  load-deflection  curves  (Figures  83 
and  81*)  is  attributed  to  the  indentation  of  the  bottom  skin  rivet  heads 
into  the  wood  support  used  to  -distribute  the  loads.  The  wooden  support 
was  used  for  tests  1  through  6.  Thereafter,  an  aluminum  support  was  used 
when  load  cells  were  installed  for  the  second  series  of  dynamic  tests. 

The  wooden  support  is  not  considered  to  have  significantly  influenced  the 
overall  results. 
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SUMMARY  OF  COMPARISON  BETWEEN  TEST  RESULTS  AND  PREDICTED  VALUES 
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Figure  83.  Predicted  Versus  Test  Load-Deflection  Curves 
for  Specimen  1. 
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Figure  8U.  Predicted  Versus  Test  Load-Deflection  Curves  for  Specimen  2 
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For  specimen  1  the  deformation  mode  started  with  the  buckling  of  the 
central  bay  (Bay  (3) #  Figure  47)#  side  skin#  and  the  two  transverse  bulk¬ 
heads  (C  and  D,  Figure  47).  The  center  bay  spacing  is  18  inches  and  the 
distance  between  the  end  stiffeners  is  also  18  inches.  As  the  deformation 
increased#  the  skin  tore  out  from  the  stiffeners  at  the  rivet  locations. 

At  this  point#  the  load  suddenly  increased  by  approximately  2000  pounds 
(Region  A,  Figure  83).  After  this  initial  tearing  incident#  the  load  in¬ 
creased  until  additional  skin  tore  out  from  other  rivet  locations.  The 
load  decreased  again,  as  is  noted  by  Region  B  in  Figure  83.  The  third 
peak  load  point  in  Figure  83  is  believed  to  have  occurred  after  tearing 
had  been  experienced  at  several  rivet  locations.  The  maximum  crushing 
strength  appears  to  have  been  encountered  in  the  second  or  third  peak  load 
value.  After  crushing#  the  load  gradually  decreases  to  approximately  5000 
pounds.  Beyond  a  deflection  of  3  inches,  which  corresponds  to  approximate¬ 
ly  half  of  the  specimen  height  (6.125  inches),  the  load  started  to  increase 
rapidly.  This  increase  is  attributed  to  increased  stiffness  developed  from 
the  trapped  material  as  the  specimen  is  crushed. 

The  predicted  value  shown  in  Figure  83  is  based  on  procedures  similar  to 
the  one  given  in  the  example  of  specimen  2.  First,  the  monolithic  failure 
stress  is  calculated;  then,  by  taking  into  account  rivet  strength/geometry 
effect,  a  modified  failure  stress  can  be  obtained.  From  the  latter  stress, 
the  failure  (or  peak)  load  is  estimated.  It  is  felt  that  the  failure  load 
by  this  approach  falls  somewhere  between  the  monolithic  failure  load  and 
the  complete  wrinkling  failure  load. 

The  discrepancies  between  the  test  result  and  the  prediction  are:  (l) 
for  the  deflection  in  the  range  between  0.2  to  1.0  inch,  the  successive 
(three  times)  abrupt  failure  mode  changes  from  monolithic  to  rivet  tear¬ 
ing  failure,  observed  during  the  test#  are  believed  to  be  responsible  for 
the  experimental  load  being  higher  than  that  obtained  from  the  analyses; 
and  (2)  after  a  deflection  of  1.5  inches  complete  tear-out  of  the  angle 
stiffener  occurred#  drastically  reducing  the  post-failure  plastic  bending 
moment,  and  in  turn  reducing  the  experimental  load  substantially  below  the 
predicted  load. 

For  specimen  2,  additional  angles  were  added  over  that  shown  in  the  basic 
drawing  (specimen  l);  1  inch  x  1  inch  x  0.032  inch  angles  were  added  out¬ 
side  of  bulkhead  stations  C  and  D.  In  the  middle  of  Bay  (3)  (see  Figure 
51),  the  skin  is  stiffened  by  two  angles,  one  on  each  side  of  the  web  (in¬ 
side  and  outside).  Eight  angles  were  added  to  the  specimen  along  the  side 
skins.  The  comparison  of  predicted  and  test  load- deflection  curves  is 
given  in  Figure  84. 

Because  the  side  skins  are  reinforced  by  angle  webs  at  both  sides,  the 
rivet  tearing  did  not  occur  prior  to  failure  load.  As  the  crushing  defor¬ 
mation  increased,  however,  the  specimen  failed  at  various  locations  (i.e., 
skin  was  torn  from  the  stiffener  web).  Free  warping  of  the  stiffener 
flanges,  assumed  in  the  post-failure  analysis  and  noted  in  Figure  66, 
justifies  the  conservative  estimate  of  rigid-plastic  hinge  moment  used  in 
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the  analysis .  In  Figure  84,  it  is  seen  that,  in  general,  close  agreement 
is  achieved  between  test  and  analysis  except  for  the  difference  between 
the  time  of  occurrence  of  the  peak  loads.  This  discrepancy  in  analysis  is 
believed  to  be  caused  by  two  factors:  (l)  the  neglect  of  stiffness  changes 
during  the  post-buckling  stage  up  to  the  failure  point,  and  (2)  the  assump¬ 
tion  of  rigid-plastic  hinge  mechanisms  while  ignoring  interactions  of  the 
complicated  stress  state. 

Specimen  3  is  identical  with  specimen  2  except  for  the  bulkhead  skins  which 
are  solid  (no  lightening  holes).  Even  though  the  initial  buckling  load 
might  be  different,  it  was  anticipated  that  the  ultimate  load  would  not 
be  noticeably  different  from  the  one  with  lightening  holes ,  since  the 
effective  width  of  the  skin  varies  very  little  at  the  failure  stress  level. 
The  comparison  of  the  experimental  load-deflection  curves  of  specimen  2 
versus  Specimen  3  (Figures  84  and  85)  shows  that  this  assumption  is  accur¬ 
ate.  Figures  66  and  67,  which  show  the  post-test  condition  of  specimens  2 
and  3,  also  indicate  similar  failure  modes  for  the  two  specimens. 

Specimen  4,  which  is  identical  to  specimen  3,  was  tested  under  dynamic 
loading  conditions  using  the  same  four-edge  boundary  condition  as  the 
first  three  static  tests.  An  840  pound  mass  was  dropped  with  a  maximum 
impact  velocity  of  17  ft/sec.  The  calibrated  load-deflection  curve,  based 
on  the  accelerometer  level  located  near  the  center  of  the  impact  head,  is 
given  in  Figure  86  together  with  the  predicted  curve.  For  specimen  4, 
the  maximum  deflection  is  about  3.75  inches  at  a  time  of  40  msec.  The 
maximum  g- level  is  estimated  to  be  approximately  31  g's,  which  is  equiva¬ 
lent  to  a  peak  dynamic  load  of  about  26,880  pounds.  The  final  deflec¬ 
tion  was  recorded  to  be  2.54  inches.  The  data  sampling  rate  used  during 
this  test  was  750  samples/second.  It  is  slightly  lower  than  desired  and, 
consequently,  may  have  resulted  in  a  lower  peak  load.  For  ail  the  sub¬ 
sequent  dynamic  tests,  a  data  sampling  rate  of  1500  samples/second  was 
used  to  ensure  obtaining  peak  values  within  the  desired  frequency  range 
(  <  300  Hz) .  Figure  86  shows  that  the  predicted  peak  failure  load  is 
approximately  6.8$  higher  than  the  test  value  using  faired  accelerometer 
data.  The  failure  deflection  point,  as  in  previous  comparisons,  is  lower 
for  the  analysis  as  compared  to  the  test. 

Specimens  5  and  6  are  the  same  as  specimen  4.  except  for  the  added  thicker 
gage  angles  between  the  bulkhead  stations  (C)  and  (d),  Figure  47,  and  for 
the  replacement  fasteners  that  attach  the  angle  to  the  main  beam  (stiffener) 
web  which  are  0.2 5- inch-diameter  screws.  The  screws  are  installed  to  pre¬ 
vent  a  separation  of  the  angles  from  the  side  stiffeners  due  to  tearing  or 
tension  failure  at  the  higher  impact  velocity  of  30  ft/sec  (specimen  5). 
Specimen  6  was  tested  at  14  ft/sec,  so  the  impact  head  weight  had  to  be 
increased  to  3550  pounds  to  ensure  that  both  specimens,  5  and  6,  were 
tested  at  the  same  energy  level.  Figures  87  and  88  show  a  comparison  of 
the  load-deflection  curves  for  specimens  5  and  6,  respectively.  As  in  the 
test  of  specimen  4,  the  load  values  were  obtained  from  faired  accelerom¬ 
eter  data.  As  can  be  noted  in  Figures  86,  87  and  88,  the  oscillations 
in  the  accelerometer  readings  make  it  difficult  to  define  the  continuous 
load-stroke  curve.  At  best,  it  can  be  ascertained  that  the  accelerometer 
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Figure  85.  Predicted  Versus  Test  Load-Deflection  Curves  for  Specimen  3 
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Figure  86.  Predicted  Versus  Test  Load-Deflection  Curves  for  Specimen  4 


response  level  is  faired  about  a  mean  curve.  All  subsequent  dynamic  tests 
(7  and  8)  were  performed  with  load  cells  installed  between  the  support  and 
the  ground.  This  setup  was  described  earlier  in  the  section  entitled  Sub¬ 
structure  Test  Program.  The  mode  of  failure  for  specimens  5  and  6  is  very 
similar,  as  can  be  noted  in  Figure  Jl.  Using  faired  accelerometer  data, 
the  difference  between  experimental  and  predicted  peak  values  for  speci¬ 
mens  5  and  6  is  between  22.5#  and  24#. 

Tests  7,  8  and  9  were  designed  to  test  the  same  configuration  of  specimens 
(except  for  lower  skin  for  specimen  9)  under  three  conditions:  27  ft/sec 
impact  velocity,  l4  ft/sec  impact  velocity,  and  statically.  In  these 
tests,  load  cells  were  installed  beneath  the  four-edge  support  and  the 
ground  surface.  The  results  of  these  tests  are  shown  in  Figures  89,  90, 
and  91*  The  load  cells  appear  to  provide  a  better  definition  of  the  load 
than  is  obtain'd  from  accelerometer  readings,  which  are  influenced  by  the 
frequencies  of  the  total  system,  including  impact  head,  carriage,  tower, 
and  specimen.  The  results  of  these  three  tests  also  show  that  for  this 
type  of  structure,  static  tests  are  sufficiently  accurate  to  determine 
load-deflection  characteristic  shape,  peak  load  and  post-failure  load. 
However,  in  dynamic  tests  the  springback  effect  is  immediate,  and  the  struc¬ 
ture  will  experience  slightly  lower  deflection  values  than  when  tested 
statically.  The  predicted  load-stroke  curve  for- the  static  load  condition 
was  used  for  all  three  cases,  and  the  peak  failure  load  varied  from  -3-9# 
to  4.5#  from  the  peak  load  obtained  from  the  load  cell  data.  The  maximum 
load  obtained  during  the  27- ft/sec  dynamic  impact  test  is  144,000  pounds, 
and  the  failure  load  experienced  in  the  static  test  is  132,000  pounds. 

The  difference  between  peak  static  and  dynamic  (27  ft/sec)  failure  load 
for  this  type  of  specimen  is  approximately  9#» 

The  comparisons  of  predicted  and  test  load-deflection  curves  for  speci¬ 
mens  10,  11  and  12  are  shown  in  Figures  92,  93  and  94.  The  discrepancy 
between  predicted  and  test  failure  loads  obtained  for  specimens  10,  11 
and  12  is  between  +1.3#  and  -8#.  A  slight  slippage  of  one  of  the  supports, 
which  was  noted  at  the  conclusion  of  the  test,  may  account  for  the 
"diagonal"  like  buckling  failure  mode  Indicated  in  Figures  80  and  82. 

COMPARISON  OF  TEST  AND  ANALYSIS  ENERGY  ABSORPTION 

Table  XX  summarizes  the  comparative  values  of  energy  absorption  and  the 
percentage  of  difference  between  test  and  analysis.  The  analysis  ranges 
from  -27.1#  (specimen  12)  to  26.3#  (specimen  5)  of  the  measured  data.  For 
specimens  1  through  4  the  analysis  is  between  -1#  and  -8.6#.  For  speci¬ 
mens  6,  7,  8,  9  and  10,  the  predicted  energy  absorbed  is  +10#,  -6#,  -lS#, 
-1.5#  and  -13#,  respectively.  The  shape  of  the  post- failure  load-deflec¬ 
tion  curve  will  greatly  influence  the  amount  of  energy  that  is  absorbed. 

The  analysis,  with  its  assumed  plastic  hinge  moment,  presents  a  lower  bound 
post-failure  curve  and  thus,  in  general,  is  expected  to  underestimate  the 
structural  energy  absorption  capability.  This  occurs  in  10  of  the  12 
specimens  tested. 
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TEST  LOAD-DEFLECTION  CURVE 


PREDICTED  LOAD-DEFLECTION  CURVE 


DYNAMIC  TEST  CONDITIONS 


STRIKING  MASS  3535  LB 
IMPACT  VELOCITY  lfc  FT/SEC 


DEFLECTION,  IN. 

Figure  89.  Predicted  Versus  Test  Load-Deflection  Curves  for  Specimen  7 
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TEST  LOAD-DEFLECTION  CURVE 


PREDICTED  LOAD-DEFLECTION  CURVE 


DYNAMIC  TEST  CONDITIONS 


STRIKING  MASS  1385  LB 
IMPACT  VELOCITY  127  FT/ SEC 


DEFLECTION,  IN. 

Figure  90.  Predicted  Versus  Test  Load-Deflection  Curves  for  Specimen  8. 
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Figure  91.  Predicted  Versus  Test  Load-Deflection 
Curves  for  Specimen  9. 
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Figure  92.  Predicted  Versus  Test  Load-Deflection 
Curves  for  Specimen  10. 


Figure  93.  Predicted  Versus  Test  Load-Deflection  Curve  for  Specimen  11 
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Figure  Predicted  Versus  Test  Load-Deflection  Curves  for  Specimen 


The  following  sub-sections  discuss  the  effect  of  the  initial  assumptions  on 
the  comparison  between  analysis  and  test  results  previously  described.  In 
addition,  empirically  based  correction  factors,  which  are  designed  to  improve 
predictions  of  failure  and  the  post-failure  load-deflection  curve,  are  pre¬ 
sented,  where  applicable. 

REMARKS  OH  THE  INITIAL  ANALYTICAL  ASSUMPTIONS 


Plastic  Hinge  Formation 

The  post-failure  analysis  is  based  on  the.  plastic  hinge  forming  in  a  sym¬ 
metrical  manner  as  is  shown  in  Figure  26a.  While  the  analytical  prediction 
based  on  such  a  mode  yields  reasonable  results,  it  has  been  observed  from 
the  tests  performed,  during  the  study  described  herein,  that  an  asymmetric 
plastic  hinge  mode  occurred  in  some  of  the  panels.  An  asymmetric  plastic 
hinge  mode,  as  shown  in  the  sketch  below,  may  result  from  structural 
imperfections, 


colunn  length 
*  plastic  hinge 
moments 


The  energy  expression  for  the  above  asymmetric  mode  system  is  as  showr 
below: 
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By  adding  Equations  (4l)  and  (42)  and  using  trigonometric  identities. 
Equation  (43)  is  obtained. 


C0S  (01  +  02>  =  (2£  “  u)2  "  +  i2)  (43) 

Equation  (43)  can  be  rearranged  as  shown,: 

$  6  (ax  +  e2)  =  -  dz/d-z)1/2  (44) 

where 

2=  |(2l  -  u)2  -  (f^  +  «/)]  /2<x£2  (45) 


Therefore 
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Equation  (47)  is  obtained  by  substituting  Equation  (46)  into  Equation  (38). 
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From  Equation (47) 
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Equation  (49)  has  the  identical  form  as  Equation  (17)  except  for  the 
weighting  term  Q(u).  It  is  concluded  that  for  asymmetric  cases  the  post¬ 
failure  load  can  be  estimated  by  multiplying  the  expression  shown  in  Equa¬ 
tion  (17)  by  the  sensitivity  factor  Q(u)  described  by  Equation  (48). 

Figure  95  shows  the  sensitivity  factor  versus  deflection  (crushing  defor¬ 
mation)  for  panel  depths  of  6  inches  and  12  Inches. 
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Figure  95.  Sensitivity  Factor  for  Asymmetric  Hinge  Formation  Versus 
Deflection  as  a  Function  of  Column  Length. 
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It  can  be  seen  that  as  either  the  panel  depth  increases,  or  the  ratio  of 
l'*2  decreases ,  the  sensitivity  factor  converges  to  the  symmetric  case 
dl  *  I2).  In  general,  the  asymmetric  modes  cause  higher  post-failure 
loads  from  10^  to  35%,  for  column  iengths  of  12  inches  and  6  inches, 
respectively,  at  about  three  inches  of  deflection.  If  for  example,  the 
sensitivity  factor  shown  in  Figure  95  is  applied  to  the  analysis  of 
specimen  10  (Figure  92)  in  which  asymmetric  failure  modes  are  noted,  the 
post-failure  curve  load  level  at  three  inches  of  deflection  would  show 
improved  agreement  with  the  test  data  with  regard  to  energy  absorption. 
However,  the  energy  absorption  correlation,  without  the  application  of  the 
sensitivity  factor,  is  within  13%  of  agreement.  The  sensitivity  factor 
would  be  applicable,  as  shown,  only  if  all  the  observed  failures  are 
asymmetrical  modes.  Since  the  test  results  indicate  the  presence  of  both 
symmetric  and  asymmetric  failure  modes,  the  basic  assumption  of  a  symmetrical 
failure  mode,  while  slightly  underestimating  energy  absorption,  is  adequate. 

Empirical  Corrections  to  the  Post-Failure  Load-Deflection  Curve 

The  post-failure  load-deflection  curve  derived  in  Equation  (17),  subsection 
entitled  Post-Failure  Analysis  of  Stiffened  Panels,  has  the  general  shape 
shown  in  Figure  96  as  a  function  of  panel  depth  and  the  value  of  the  plastic 
hinge  moment  that  is  formed. 


Figure  96.  Post-Failure  Load -Deflection  Curves  . 

In  the  Substructure  Analysis  section  the  plastic  hinge  moment  (Mp)  given 
in  Equation  (17)  is  based  on  the  most  conservative  geometry  shape,  or 
minimum  value,  which  corresponds  to  the  condition  in  which  the  stiffener 
flange  folds  completely  over  due  to  warping.  This  condition  results  in  a 
much  lower  bending  stiffness  than  would  be  obtained  if  the  flange  is  ex¬ 
tended.  Although  the  correlation  of  analysis  and  test  shows  good  agreement 
with  regard  to  load-deflection  characteristics,  the  analysis  yields  a  con¬ 
sistently  lower  deflection  value  at  failure  when  compared  to  test  data.  In 
part  this  disci epancy  may  be  due  ^o  the  formation  of  asymmetric  hinge  modes 

125 


V 


* 

I 

I 

I 

1 

1: 


r 


as  opposed  to  the  formation  of  symmetric  hinge  modes,  as  initially  assumed. 
It  is  estimated  that  the  formation  of  an  asymmetric  failure  mode  in  lieu 
of  a  symmetric  failure  mode  will  influence  the  failure  load  predicted  by 
analysis  with  a  shift  in  the  failure  deflection  point. 


Thus,  the  primary  difference  between  the  failure  deflection  obtained  by  anal¬ 
ysis  and  test  is  attributed  to  the  fact  that  the  original  analysis  (Equa¬ 
tion  17)  assumes  an  overly  conservative  hinge  moment.  This,  for  example, 
would  be  reflected  in  post-failure  curve  A  of  Figure  96.  The  less  conserva¬ 
tive  the  hinge  moment,  the  more  the  post-failure  curve  shifts  to  the  right 
(B  and  C).  The  test  results  indicate  that  the  plastic  hinge  moment  could 
be  higher  for  the  transitional  range  (at  and  immediately  after  failure). 

The  same  phenomenon  is  noted  in  Reference  4l,  which  based  the  post-failure 
analysis  on  pure  bending.  The  following  simple  idealization  of  such 
behavior  in  view  of  the  experimented  data  is  suggested  in  Reference  1+1. 


M  =  M  ,  +  (M  -  M  .  )  e"K9  (50) 

p  min  '  max  min7  '  ' 

where  K  is  constant  and 


9  is  a  plastic  hinge  rotation  in  radians 

Mmin  13  the  most  Pi33**3  "ir*.  moment 

M  is  the  momfent  for  the  undeformed  cross-section 

H1SLX 


However,  no  attempt  is  made  to  determine  K  in  Reference  1+1,  In  the  fol¬ 
lowing  discussion,  the  parameters  which  influence  the  value  of  K  are 
assessed.  In  the  rigid  perfectly  plastic  hinge  moment,  the  only  variables 
which  govern  the  value  of  K  in  Equation  ( 50)  are  the  section  modulus  and 
the  column  length.  Thus,  the  exponential  function  of  e"K®  should  have  a 
decreased  section  modulus  value.  For  the  typical  angle-type  stiffener 
arrangement,  as  shown  in  Figures  2l+(c)  and  24(d),  the  following  simple 
form  of  constant  K  is  proposed. 
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Correlation  with  the  test  data  obtained  in  this  study  suggests  K  as 


Although  the  expression  for  S  given  in  Equation  (51)  has  been  limited  to 
application  to  angle  type  stiffened  panels,  it  may  be  applicable  to  other 
types  of  stiffeners.  Until  additional  crushing  test  data  is  evaluated, 
however,  the  expression  for  K  noted  above  should  be  limited  to  angle 
stiffeners. 

The  use  of  Equation  (51)  is  extended  to  zee  stringer  stiffened  panels  by 
using  a  value  of  a  =  10  in  Equation  (51)  based  on  post-failure  load- 
deflection  data  deduced  from  test  data  presente'’,  in  Reference  4.  Fig¬ 
ure  97  shows  the  experimental  load-deflection  curve  of  a  stiffened  panel 
with  zee  stringers  versus  the  predicted  curve  using  the  analytical  method, 
modified  by  Equation  (51)*  The  comparison  between  test  and  analysis 
matches  reasonably  well  up  to  a  vertical  crushing  deflection  of  0.2  in¬ 
ches.  However,  beyond  this  deflection,  especially  above  1.0  inch,  the 
predicted  load  increases  to  approximately  twice  the  load  obtained  from 
testing.  This  discrepancy  in  load  is  attributed  to  the  relaxation  of  the 
end  constraints  of  the  test  setup  in  the  large  post-failure  range,  thereby 
approaching  a  simply  supported  end  condition.  Photographs  of  the  test 
specimens  in  the  post-test  condition,  presented  in  Reference  4,  show  evi¬ 
dence  of  end  restraint  relaxation. 

If  the  end  restraints  remained  fixed  it  is  anticipated  that  the  test  load 
at  a  deflection  >  1  inch  would  be  approximately  twice  the  load  level  shown 
in  Figure  97  due  to  additional  hinge  formations  at  the  supports.  There¬ 
fore,  the  analytically  predicted  values  shown  in  Figure  97  are  expected 
to  show  good  correlation  with  Reference  4  test  data,  assuming  the  end 
restraints  are  held  fixed. 


Failure  Deflection 

Determination  of  the  Deflection  Value  at  Failure 

Although  extensive  U3e  of  semi -empirical  approaches  has  been  made  in  this 
study  to  predict  the  buckling  and  failure  load,  the  literature  evaluation 
indicates  that  information  regarding  failure  deflection  and  pre-failure 
stiffness  is  very  limited  with  the  exception  of  simple  configurations  such 
as  flat  plates.  Figure  98,  obtained  from  Reference  107,  shows  a  typical 
curve  of  changes  in  stiffness  versus  applied  load  for  a  square  plate.  The 
relationship  between  the  structural  imperfection  parameter,  a,  and  plate 
thickness  (ts),  initial  lateral  deflection  (W0),  critical  stress  (ocr) ,  and 
yield  stress  (ocy)  is  given  in  Reference  107  by  the  following  expression: 
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Figure  98.  Variation  of  Stiffness  with  Load  (Reference  107). 


For  example,  from  Figure  98  it  can  be  seen  that  as  the  initial  structural 
imperfection  become  greater  (increasing  or)  the  prebuckling  stiffness 
( P/ pcrit  -i)  is  Bore  strongly  affected  than  the  post-buckled  stiffness 
(P/pcrit  -l)*  Figure  98  shows  changes  in  stiffness  versus  P/pcrit  ^or  a 
flat  square  plate  in  which  crcr/<rcr  =  0.25  and  a  =  0.3,  and  is  simply  sup¬ 
ported  with  no  edge  restraint.  As  a  increases  it  is  expected  that  in  the 
limit,  the  stiffness  in  the  pre-buckled  region  will  nearly  equal  the 
stiffness  in  the  post-buckled  region.  Thus,  the  greater  the  initial  im¬ 
perfection  the  more  it  would  be  expected  that  a  constant  stiffness  will 
be  maintained  up  to  the  failure  point.  However,  the  prediction  of  the 
stiffness  for  complex  structures  up  to  failure  is  difficult  to  ascertain 
accurately  by  simple  techniques  due  to  the  sensitivity  6f  the  stiffness 
of  these  structures  to  initial  imperfections  and  their  nonlinear  response. 

The  test  data  for  the  specimens  tested  during  this  study  indicate  that 
there  may  be  significant  imperfections  in  the  structure,  in  that  the  ini¬ 
tial  slope  of  the  load-deflection  curve  in  the  pre-buckled  region  is  not 
markedly  different  from  that  in  the  post-buckled  region  up  to  the  failure 
load. 

For  the  analytical  method  developed  during  this  study  and  described  in  the 
Substructure  Analysis  section,  the  prefailure  stiffness  is  obtained  from 
the  intersection  of  the  predicted  failure  load  and  the  post-failure  analysis 
curve.  This  point  determines  the  failure  point  deflection.  The  correla¬ 
tion  between  analysis  and  test  shows  that  the  analysis  consistently  under¬ 
estimates  the  failure  deflection  (Figures  83  through  94)  although  predicted 
energy  absorption  and  failure  load  values  show  reasonable  agreement  with 
test  data.  In  general  the  test  results  show  failure  deflections  of  .20 
to  .35  inches  whereas  the  analytically  predicted  failure  deflection  values 
vary  from  .12  to  .25  inch.  Table  XXI  is  a  comparison  of  the  deflection 
at  failure  obtained  by  test  and  analysis  for  each  set  of  test  specimens. 


*  Average  Value 
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Unless  an  extensive  analysis,  perhaps  requiring  involved  computer  coding, 
is  performed  it  is  expected  a  semi-empirical  approach  using  simplified 
analytical  techniques  would  not  ensure  significantly  better  agreement  in 
determining  failure  deflection  than  is  obtained  using  the  approach  adopted 
in  this  study.  Implicit  in  this  approach  is  a  presumption  that  yield 
strain  is  reached  at  the  point  of  failure.  The  failure  criterion  which 
most  investigators  adopt  for  collapse  analysis  is  in  fact  an  upperbound 
solution  for  failure.  This  can  be  shown  by  the  following  simplistic  ap¬ 
proach  to  an  analysis  of  a  flat  plate  using  von  Karman's  plate  equations. 


(55) 


However, 


the  upper  bound  value  of 


u 


can  be  given  as 
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upper  bound 


c  dx 
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(56) 


where 


u  =  total  end  shortening 
w  =  lateral  deflection 
ex  =  axial  strain 
i  =  length  of  web 

=  partial  derivative  of  w  with  respect  to  the  coordinate  x 


In  view  of  the  fact  that  the  failure  stress  can,  at  most,  be  expected  to 
reach  the  yield  stress  it  follows  that  the  upper  bound  value  of  ex  can  be 
considered  to  equal  the  yield  strain  (eTO)  and,  therefore,  the  maximum 
probable  value  of  ex  can  be  approximateaPas  equal  to  t  .£ 

yp 


The  analysis  described  in  this  study  is  based  on  the  assumption  that  full 
plastic  hinge  moments  are  developed  at  the  threshold  of  failure  which 
implicitly  suggests  that  the  yield  strain  («  )  has  been  reached.  In 

principle  this  assumption  is  equivalent  to  determining  the  upper  bound 
end- shortening  using  the  general  collapse  criteria  which  assumes  the  stress 
at  the  supports  has  reached  yield.  For  all  practical  purposes  of  obtain¬ 
ing  approximate  energy  absorption  capability  and  characteristic  load- 
deflec.tion  shapes  for  crashworthiness  analysis,  the  approach  adopted  in 
this  study  is  reasonable,  especially  when  it  Is  realized  that  the  pre¬ 
failure  deflection  region  accounts  for  only  a  very  small  fraction  of  the 
total  crushing  deformation. 
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The  present  method  adopts  a  modified  conservative  plastic  hinge  moment  in 
line  with  a  similar  observation  made  in  a  study  described  in  Reference  108 
in  which  it  is  stated  "if  the  stiffener  is  riveted  to  che  web,  -  -  -  -  the 
bending  stiffness  that  comes  into  play  is  more  nearly  that  of  the  attached 
leg  alone,  rather  than  that  of  the  entire  stiffener  -  -  Therefore, 
the  end- shortening  at  the  failure  load,  determined  by  the  present  method, 
will  in  fact  be  an  intermediate  value  sanewhere  between  the  maximum  and 
minimum  possible,  plastic  hinge  moments  that  can  be  developed.  If  the 
maximum  possible  plastic  hinge  momept  is  to  be  used  in  the  analysis  it 
would  result  in  obtaining  an  upperbound  failure  point  deflection  which 
would  show  better  agreement  with  failure  deflection  values  obtained  from 
specimen  testing.  However,  the  use  of  the  maximum  hinge  moment  wculd  then 
result  im  an  overestimation  of  the  post-failure  load  and  energy  absorption. 
Therefore,  the  present  approach, alt hough  simplifying  the  complex  structural 
behavior  involved  with  regard  to  the  prediction  of  the  failure  deflection, 
does  account  for  the  dominating  crushing  modes  observed  from  the  test  per¬ 
formed  during  this  study. 

Strain  Data  and  Failure  Modes 


The  strain  gages  are  installed  to  detect  the  modes  of  failure  and  to  verify 
the  assumptions  introduced  in  the  analysis.  Strain  gage  data  has  limita¬ 
tions  when  seeking  quantitative  measurements  of  large  deflection  failures. 

In  addition,  wrinkling  and  interrivet  tearout  make  it  difficult  to 
accurately  interpret  strain  values.  On  the  first  series  of  tests,  seven 
strain  gages  were  used.  The  number  of  strain  gages  was  reduced  to  two  for 
the  last  six  tests  after  a  review  of  the  initial  strain  data  showed  that 
the  failure  modes  assumed  in  the  analysis  are  valid.  For  example,  Figure  99 
shows  load  versus  strain  for  the  strain  gages  installed  on  specimen  1. 

Gages  1A  and  IB  mounted  on  the  outer  and  innfer  sides  of  the  side  skin  re¬ 
sponded  to  strain  in  the  same  direction  as  a  result  of  rivet-tearout .  In 
subsequent  tests  these  two  gages  responded  to  both  compression  and  bending. 

In  this  particular  test  the  strain  exceeds  0.0048  in. /in.  which  is  the  yield 
point  for  2Q24-T3  material.  The  curve  for  strain  gage  2  shows  that  up  to 
the  failure -load  (l4,000  lb)  there  is  an  increasing  compressive  stress  to 
approximately  30,000  psi  (0.003  in. /in.)  at  the  stiffener  flange  surface. 
After  the  failure  load  is  reached  there  is  an  apparent  change  in  failure 
pattern.  Strain  gage  3  was  located  at  a  point  where  the  side  skin  tore  out 
from  the  rivets. and  no  meaningful  interpretation  is  given  to  the  data  be¬ 
yond  the  failure  load.  The  readings  from  strain  gage  4  indicate  that  the 

strain  at  that  particular  location  remained  within  the  elastic  range  (c  = 

max 

0.0015  in. /in.)  during  the  deformation  of  the  stiffener.  This  response 
tends  to  support  the  contention  resulting  from  analyses  of  the  main  beam  and 
bottom  skin,  that  a  built-in  support  exists  at  the  location  of  this  strain 
gage.  Strain  gage  4,  installed  for  the  purpose  of  detecting  a  plastic 
hinge  formation  at  its  location,  recorded  low  strain  levels  (»0.0021  in. /in.) 
at  the  time  of  initial  failure .  As  bending  increased  the  plastic  hinge 
zone  increased, as  is  noted  bv  „he  increase  in  tensile  strain  as  the  load 
decreases.  Observations  of  the  specimen  during  the  test  and  the  post-test 
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STRAIN,  (jl IN. /IN. 

Figure  99.  Load  Versus  Strain  for  Specimen  No. 


inspection  indicated  that  a  plastic  hinge  formed  approximately  six  inches 
from  strain  gage  location  5.  Strain  gage  6  underwent  very  little  strain 
throughout  the  deformation  of  the  stiffener.  The  maximum  strain  recorded 
is  0.0004  in. /in.  ( *  40001b/in.^).  The  low  strain  reading  at  this  location 
supports  the  assumption  of  a  double  plastic  hinge  model  in  the  bottom  beam 
and  skin  analysis. 

The  following  discussion  provides  some  insight  into  the  potential  modes  of 
failure  and  the  difficulty  associated  with  accurately  interpreting  strain 
gage  data.  In  the  subsequent  specimen  tests  (2  through  11),  in  which 
additional  angle  stiffeners  are  added  to  the  side  panels,  strain  gages  1A 
and  IB  are  located  on  the  outer  and  inner  surfaces ,  respectively,  of  the 
angle  stiffener  legs.  Figure  100  shows  exaggerated  local  deformation  modes 
exhibited  by  the  side  panel,  for  small  and  large  post-failure  deflections. 
Deformation  Mode  A  is  shown  by  the  solid  lines  while  deformation  Mode  B 
is  shown  by  the  dashed  lines.  In  the  case  of  Mode  A  the  stiffeners  move 
together, which  places  strain  gage  1A  in  tension  and  strain  gage  IB  in 
compression.  On  the  other  hand.  Mode  B  is  formed  by  the  two  stiffeners 
moving  in  opposite  directions.  In  this  case  both  strain  gages  will  be 
placed  in  tension.  The  direction  in  which  the  stiffeners  move  under  the 
compression  load  is  a  function  of  the  residual  strains  incorporated  during 
manuf acture ;  for  example,  in  the  installation  of  the  fasteners  at  x  and  y 
in  Figure  100.  A  change  in  mode  shape  could  also  take,  place  depending  on 
the  influence  of  the  global  modes  on  the  initial  local  modes. 

Table  XXII  presents  a  summary  of  the  stiffener  response  of  some  of  the 
specimens  employed  in  the  tests.  The  data  used  to  develop  the  table  are 
contained  in  figures  presented  in  Volume  II  in  the  section  entitled  Test 
Data. 

Behavior  of  the  stiffeners  during  pre-failure  and  post-failure  applica¬ 
tion  of  the  load  is  random.  That  is,  the  fact  that  the  stiffeners  may 
move  together  prior  to  failure  does  not  mean  that  'they  will  move  together 
after  failure  or  vice-versa. 


TABLE  XXII.  SUMMARY  OF  TEST  SPECIMEN  STIFFENER  RESPONSE 
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PROGRAM  "KRASH"  REFINEMENT 


GENERAL 

Program  KRASH  was  developed  as  an  analytical  tool  to  be  used  during  pre¬ 
liminary  phases  of  design.  In  Reference  1  the  results  of  a  parametric 
study  are  integrated  into  an  iterative  design  procedure  by  which  a  trade¬ 
off  between  potential  incremental,  cost  and/or  weight  versus  incremental 
improvements  in  crashworthiness  capability  could  be  accomplished.  The  re¬ 
sults  of  that  study  and  a  subsequent  study  for  the  U.S.  Army  (Reference 
106)  showed  the  potential  benefits  that  could  be  achieved  for;  designers 
using  a  comprehensive  but  unsophisticated  analytical  approach.  As  is  the 
situation  with  newly  developed  advances  in  the  state-of-the-art,  there  are 
generally  areas  in  which  the  analytical  techniques  can  be  simplified  for 
the  user. 

In  an  effort  to  facilitate  a  designer’s  usage  of  program  KRASH  and  to  in¬ 
corporate  the  results  of  the  literature  survey,  load  sensitivity  study, 
substructure  analysis,  and  tests,  the  program  is  revised.  1“  particular, 
the  input  data  format  is  changed  for  ease  of  data  input  and  subsequent 
data  changes,  as  would  be  required  during  parameter  tradeoff  studies.  In 
addition,  flexibility  is  added  in  the  manner  by  which  the  stiffness  reduc¬ 
tion  factors  (KR's)  are  input.  The  revised  program  was  run  for  the  follow¬ 
ing  conditions: 

•  Correlation  case  31-52  from  Reference  l(  to  demonstrate  that  the 
revised  program  format  is  compatible  with  previous  results:  23 
ft/sec  vertical  velocity  combined  with  a  18.5  ft/sec  lateral 
velocity), 

•  Three-dimensional  velocity  ( ~  40  ft/sec  combined  velocity  Impact  - 
27.5  ft/sec  longitudinal,  23  ft/sec  vertical,  18.5  ft/sec  lateral). 

•  Upper  mass  penetration  into  a  specified  occuplable  volume-. 

•  Simplified  blade  contact. 

e  Utilization  of  load-deflection  data  obtained  from  the  12  substruc¬ 
ture  tests,  performed  during  the  study,  and  related  to  actual  fuse¬ 
lage  structure  size. 

* 

Volume  II  contains  a  revised  input-output  format  description  as  well  as 
|  two  sample  problems.  One  sample  problem' represents  the  updated  correla- 

j  .  tion  run.  Case  31-52,  and  the  second  sample  problem  corresponds  to  the 

[  three-dimensional  velocity  impact  condition.  Also  included  in  Volume  II 

l  is  a  write-up  of  the  energy  balance  equations  and  a  sample  of  the  format 

|  in  which  energy  results  are .provided.  The  energy  balance  data  is  current- 

<  ly  in  the  program  as  an  option  since  only  limited  checkout  of  the  equations 

lias  been  performed . 
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PROGRAM  "KRASH"  INPUT  FORMAT  REVISIONS 

The  input  format  changes  are  divided  into  three  categories: 


•  Reordered  data 

•  Standardization  of  certain  inputs 

•  Allowance  for  more  general  'KR  curves 

First,  the  input  data  is  rearranged  so  that  all  mass  associated  data  (sub¬ 
script  i)  is  in  one  block,  followed  by  all  external  spring  data  (subscript 
ik)  and  then  all  internal  beam  data  (subscript  ij).  Next,  standard  values 
are  assigned  to  certain  input  data  items  unless  otherwise  specified.  Thus 
much  repetitive  data  input  is  eliminated.  The  following  quantities  are 
given  standard  values  unless  otherwise  specified: 


QUANTITY 

SYMBOL 

STANDARD 

VALUE 

Angular  Momenta 

Hej. 

0 

Euler  Angles 

01,  8 », 

0 

Aerodynamic  x,ift 

lCi 

0 

Stiffness  Reduction 
Factors  . 

KRij 

1 

Failure  Inflection 

Vraaxijl 

100 

Thus,  the  angular  momenta  of  the  masses  (He^)  and  the  aerodynamic  lift 
(lei)  are  normally  zero ,  as  are  the  Euler  angles  9^,  relating  mass 

body-fixed  axes  to  airplane  c.g.  axes.  For  linear  internal  beam  elements, 
KR  =  3-  for  the  entire  run,  so  only  nonlinear  KR  data  need  be  input. 
Similarly,  Vmax’g  need  only  be  Input  for  those  elements  where  rupture  is 
expected  to  occur;  for  elements  not  specified,  a  failure  deflection  of  100 
inches  and  rotation  of  100  radians  are  assumed.  Thus,  the  nonspecified 
beams  will  not  rupture  during  the  run. 

The  input  has  also  been  revised  to  allow  more  general  KR  (stiffness  reduc¬ 
tion  factor)  curves.  Previously,  each  KR  curve  could  have  only  six  data 
points,  equally  spaced.  Now  each  curve  can  have  up  to  15  data  points, 
with  any  desired  spacing.  !Ms  greatly  facilitates  modeling  more  com¬ 
plex  load-deflection  relationships.  However,  the  load-sensitivity  study 
results  indicate  that  for  many  structural  elements  r  very  simple  KR  repre¬ 
sentation  is  adequate.  The  KR  representation  for  the  various  load- 
deflection  categories  is  discussed,  in  detail,  in  the  Design  Procedures 
section. 


As  a  result  of  the  program  refinement  and  rearrangement  of  coding,  the 
program  capacity  has  been  expanded  to  treat  at  least  80  lumped  masses 
and  100  Internal  beam  elements  simultaneously.  Previous  capacity  was  40 
and  80,  respectively. 

THREE-DIMENSIONAL  IMPACT  AND  MASS  PENETRATION 
Case  Description 

After  program  KRASH  was  refined, as  described  previously,  correlation  run 
31-52  from  the  preceding  study  (Reference  l)  was  performed,  and  the  re¬ 
sults  indicate  that  the  program  conversion  is  successful,  Case  31-52, 
run  with  the  refined  program,  is  presented  as  a  sample  problem  in  Volume  II. 

The  basic  31  mass  model  described  in  Reference  1  is  modified  by  adding  a 
mass  32  to  represent  a  lumped  mass  on  the  rotor  blade,  with  the  blade  in 
the  fore-aft  position,  as  shown  in  Figure  101.  '  In  order  to  obtain  a  mass 
penetration  of  mass  32  into  the  cabin  area  within  a  reasonable  run  time, 
the  rotor  is  lowered  and  tilted  so  that  mass  32  starts  out  just  above  the 
helicopter  cabin.  Figure  101  is  a  left-hand  side  view  of  the  model;  the 
dashed  lines  represent  the  mass  penetration  control  volume,  which  is  refer¬ 
enced  to  mass  9*  Mass  32  is  located  just  above  the  control  volume 
since  the  purpose  of  the  computer  run  is  to  demonstrate  mass  penetration 
capability  and  not  to  correlate  with  any  particular  accident  condition. 

The  top  surface  of  the  control  volume,  which  mass  32  penetrates,  is  pro- 
perly  located. 

The  initial  c.g.  velocity  is  chosen  as  40  ft/sec,  vectored  as  follows: 

Vx  =  27  ft/sec 
V  =  18.6  ft/sec 
V_  =  23  ft/sec 

|U 

Thus,  a  three-dimensional  impact  with  a  horizontal  surface  is  simulated. 

The  initial  vehicle  attitude  is  flat  (zero  roll,  pitch  and  yaw  angles), 
and  the  initial  angular  velocities  are  zero. 

Results 

Mass  32  penetrates  the  control  volume  at  ,086  second  after  initial  impact. 
The  aft  skids  fail  at  .056  and  the  forward  skids  at  .062,  The  fuselage 
first  contacts  the  ground  just  after  skid  failure,  the  aft  fuselage  con¬ 
tacting  just  slightly  before  the  forward  fuselage.  After  fuselage  con¬ 
tact,  the  engine  and  transmission  accelerations  build  up  to  their  peaks 
as  follows: 


' 

Vertical 

Lateral 

Item 

Accel,  (g's) 

Accel.  Cg*s) 

• 

Engine 

Transmission 

31.3  9  .097 

30.0  §..089 

7.3  @  .087 

3.5  e  .096 

Longitudinal 

Accel.(g's) 

5.2  @  .086 
2.3  ®  .088 


The  engine  mounts  deform  nonlinearly  vertically;  the  transmission  mount 
goes  nonlinear  longitudinally.  Lateral  deformations  remain  linear. 

ROTOR  BLADE  CONTACT 


Case  Description 

A  rotor  blade  is  added  to  the  basic  31  mass,  38  beam  model  shown  in  Refer¬ 
ence  1,  resulting  in  a  36  mass  model  with  43  internal  beams.  Figure  102 
shows  the  revised  configuration.  Masses  32  through  36  were  added  to  repre¬ 
sent  the  rotor  hub  and  blade,  with  the  blade  stationary  and  perpendicular 
to  the  plane  of  the  longitudinal  axis  of  the  helicopter.  Figure  102  shows 
a  rear  view  of  the  right  side  of  the  vehicle.  The  rotor  blade  is  repre¬ 
sented  with  masses  lumped  at  33  end  34.  External  springs  in  the  y  and  z 
direction  are  placed  at  masses  34  (rotor  tip),  14  (skid  end),  and  12  and 
22  (lower  right  corners  of  fuselage). 

The  vehicle  is  initially  rolled  35  degrees  to  the  right,  so  the  right 
rotor  tip  and  landing  skid  contact  the  ground  simultaneously.  The  initial 
c.g.  velocity  vector  is  30  ft/sec  perpendicular  to  the  ground. 

Results 

The  only  beam  elements  which  fail  are  the  right-hand  landing  skids,  ele¬ 
ments  10-14  and  16-17.  These  fail  at  .055  and  .076  second,  respectively. 
After  their  failure,  the  fuselage  contacts  the  ground  at  mass  12  (aft  end) 
at  .086.  Within  the  run  time  analyzed  (.1  second),  the  rotor  blade  did 
not  fail.  The  maximum  vertical  deflection  of  the  rotor  blade  is  about  38 
inches  at  the  tip.  The  stiffness  and  strength  of  the  rotor  blade  are  only 
modeled  approximately.  The  intent  of  the  case  is  to  illustrate  that  rotor 
blade  contact  can  be  modeled,  not  to  determine  precisely  what  happens  to 
the  vehicle  in  this  situation. 

TEST  DATA  LOAD -DEFLECTION  CURVE 

Case  Description 

A  computer  run  was  made  to  utilize  load -deflect ion  information  obtained 
from  the  specimen  testing  performed  during  the  earlier  stages  of  the 
program. 

A  modified  load-deflection  curve  is  shown  in  Figure  103.  The  load -deflec¬ 
tion  characteristics  indicate  a  rapid  load  buildup  until  failure  load  is 
reached. 
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Thereafter,  the  load  decreases  until  a  near-constant  level  is  reached. 
Finally,  bottoming  occurs  as  the  structure  crushes.  This  curve  compares 
favorably  with  the  analytical  load-deflection  curves  used  in  Reference  1 
except  for  the  sharp  peak  at  the  failure  load. 

The  test  data  is  related  to  actual  fuselage  structure  design  (size,  thick¬ 
ness,  spacing,  number  of  angles)  and  is  incorporated  into  the  UH-1H  lumped 
mass  model  (Figure  7)  for  fuselage  masses  10  (aft),  11  (midsection),  and 
16  (forward). 
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Figure  103.  Test  Specimen  Data  Load -Deflection  Curve. 


Results 

The  responses  for  the  engine,  transmission  and  occupants  (DRJ's)  using 
the  test  specimen  related  load-deflection  data  are  compared  to  the  corres¬ 
ponding  parameter  response  for  the  correlation  results  from  Reference  1 
und  test  data  from  Reference  1.  These  results  are  given  in  Table  XXIII. 
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Table  mil.  COMPARISON  OF  RESULTS  USING  TEST  SPECIMEN  DATA  WITH 
PREVIOUS  CORRELATION  ANALYSIS  AND  TEST  RESULTS 

■ — - . 

Test  Specimen 
Data  Results 

Reference  1 
Correlation 
Analysis  Results 

Reference  1 
Test  Data 
Results 

Accel.* 

Time** 

Accel.* 

Time** 

Accel.* 

Time** 

Engine  Vertical 

23.8 

138 

26 

129 

26 

120 

Transmission  Vertical 

26.0 

113 

30 

117 

27 

no 

DRI-Forward 

20.0  . 

140 

22 

142 

- 

- 

DRI-Aft 

24.5 

_ 

145 

27 

1^8 

30 

120 

*Acceleration  in  g's 

**Time  in  milliseconds  after 

initial  impact 

The  comparison  of  results  shows  that  compared  to  previously  presented 
analytical  results,  the  responses  change  less,  8.5$  to  13.6 $,  as  far  as 
response  accelerations  are  concerned  and  that  there  is  a  shift  of  2  to  9 
milliseconds  in  the  time  of  occurrence  of  events.  For  brevity,  only  sev¬ 
eral  prime  comparisons  are  made.  The  deflections,  deformations,  lateral 
forces,  etc.,  are  changed  less  than  those  parameters  noted.  The  largest 
change  (13.6$)  is  for  the  transmission,  and  this  actually  improves  the 


comparison  with  test  data.  The  engine  vertical  response  is  still  within 
8.5$  of  the  test  data. 


DESIGN  PROCEDURES 


SPECIFICATION  OF  CRASHWORTHINESS  CAPABILITY 

New  helicopter  designs  are  currently  required  to  provide  occupant  surviv¬ 
ability  and  structural  integrity  to  specified  crash  environment  levels  and 
peak  acceleration  forces.  However,  to  verify  a  design  capability  would 
literally  require  that  crashworthiness  capability  be  demonstrated  by  a 
number  of  representative  tests.  The  cost  of  testing  a  number  of  vehicles 
is  prohibitively  high  and,  at  best,  provides  one  data  point  per  test.  On 
the  other  hand,  compliance  by  showing  crashworthiness  capability  with  a  few 
simple  calculations  for  a  vertical  impact  condition  using  unverified  trans- 
missibility  factors  through  complex  structure  from  ground  to  occupant  or 
structure  is  also  not  a  satisfactory  solution.  Therefore,  there  is  a  need 
for  a  compromise  in  the  manner  in  which  compliance  is  shown. 

Crashworthiness  capability  can  be  specified  by  the  responsible  agencies  by 
defining  impact  velocity  envelopes  and/or  peak  accelerations  for  which  air¬ 
frame  manufacturers  would  then  be  required  to  design  their  vehicle.  Since 
these  levels  would  presumably  be  related  to  percentile  accident  levels  based 
on  existing  crash  data,  they  would  reflect  a  certain  acceptable  crash  con¬ 
dition.  ■  This  approach  has  the  advantage  of  at  least  providing  a  design 
goal.  However,  it  has  the  potential  to  create  a  situation  in  which  the 
manufacturer  will  overdesign,  underdesign,  or  both,  depending  on  the  location 
of  the  occupants  relative  to  major  mass  items.  The  data  from  which  the  de¬ 
sign  criteria  are  obtained  could  very  easily  be  adequate  for  one  type  of 
vehicle  design  while  not  practical  for  other  types.  For  example,  helicopter 
types  vary  from  single  powerplant  lightweight  utility  (6000  lb)  to  triple 
powerplant  heavy  lift  helicopters  ( 116,000  lb)  having  greatly  different 
occupant  accommodations  and  airframe  construction. 
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A  second  alternative  is  to  allow  the  aircraft  manufacturers  to  specify  the 
level  to  which  tneir  vehicles  are  capable  of  providing  crashworthiness  capa¬ 
bility.  In  so  doing,  the  manufacturer  would  then  have  to  describe,  in  de¬ 
tail,  the  procedures  used  to  show  that  the  capability  can  be  achieved. 
Employing  this  approach,  the  manufacturer  would  then  have  the  option  to  take 
exception  to  a  part  of  the  requirement,  particularly  if  he  could  show  that 
the  basic  objective  of  occupant  survivability  is  attained.  As  part  of 
the  verification  of  crashworthiness  capability,  the  manufacturer  would  be 
asked  to  provide  specimen  test  data  to  substantiate  structural  load- 
deformation  characteristics  in  critical  regions, 

A  third  approach  is  to  specify  that  certain  procedures  be  followed.  For 
example,  the  iterative  design  process  described  in  Reference  1  provides 
for  a  consistent  crashworthiness  design,  in  which  compatibility  between 
the  essential  structural  elements  is  maintained,  and  takes  into  considera¬ 
tion  the  crash  impact  environment,  human  tolerance  limits,  and  incremental 
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weight,  cost  and  configuration  penalties  while  treating  the  vehicle  and 
its  occupants  as  an  integrated  system.  A  procedure  of  this  nature  can  be 
applied  by  the  responsible  agency,  the  manufacturer,  or  an  impartial 
designated  party.  As  in  the  previous  approach,  limited  test  verification 
could  be  required.  Of  paramount  importance  in  this  approach  is  the  ability 
to  specify  the  trade-off  between  penalty  (weight,  cost,  configuration)  and 
crashworthiness  capability  early  in  the  design  process.  At  this  early 
stage,  design  changes  which  would  not  interfere  with  the  primary  operational 
requirement  could  be  more  economically  incorporated.  This  approach  offers 
the  advantage  of  treating  the  aircraft,  including  occupants,  as  *  total 
system,  which  would  have  the  effect  of  minimizing  the  tendency  to  over- 
design  in  order  to  meet  arbitrary  design  criteria. 

GENERAL  APPROACH 


The  determination  of  load-deflection  characteristics  of  aircraft  structures 
will  enhance  the  ability  of  designers  to  predict  structure  and  occupant 
response  during  severe  yet  potentially  survivable  accidents.  However,  it 
is  important  that  this  data  be  presented  in  a  manner  that  will  aid  the 
designer  in  developing  the  desired  level  of  crashworthiness  capability  in 
the  airframe  stucture.  In  order  to  ascertain  crashworthiness  capability, 
it  is  necessary  to  define  the  crash  environment,  to  understand  the  crush¬ 
ing  characteristics  of  airframe  structure  and  the  potential  structural 
failure  modes,  and  to  have  an  analytical  technique  which  is  applicable  to 
realistic  crash  situations.  The  design  procedures  presented  herein  provide 
the  basis  for  the  formulation  of  a  design  manual  by  which  designers  can 
readily  obtain  pertinent  data  necessary  to  perform  a  crashworthiness 
analysis.  The  information  presented  in  this  section  follows  from  earlier 
developments  described  in  this  study  and  the  preceding  study  (Reference  1) 
as  well  as  information  contained  in  the  current  Crash  Survival  Design 
Guide  (Reference  2).  The  development  of  design  procedures,  to  be  meaning¬ 
ful,  must  form  a  part  of  an  overall  process  by  which  crashworthiness 
capability  can  be  specified. 

The  development  of  a  mathematical  model  to  perform  crashworthiness  analysis 
for  a  particular  helicopter  configuration  during  the  initial  stages  of  the 
design  formulation  requires  that  consideration  be  given  to  the  following. 

«  Identification  of  the  overall  failure  modes. 

•  '  Identification  of  the  individual  element  failure  modes. 

•  Prediction  of  the  load-deflection  curve  for  each  subelement. 

•  Prediction  of  the  total  load-deflection  curve  for  each  sub¬ 
structure  . 

•  Simplified  representation  of  structural  element  load-deflection 
curves  while  providing  response  values  within  acceptable  accuracy. 

By  following  the  approach  outlined  above,  the  designer  will  have  an  oppor¬ 
tunity  to  perform  meaningful  trade-off  studies  between  cost  and/or  weight 
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versus  crashworthiness  capability.  Although  this  approach  is  not  neces¬ 
sarily  rigorous,  it  will  be  consistent  with  the  amount  and  quality  of  de¬ 
tailed  data  that  is  available  during  a  preliminary  design  study.  Further¬ 
more,  if  properly  evaluated,  the  analysis  of  preliminary  designs  will  serve 
such  useful  functions  as:  (l)  ascertaining  critical  design  regions  wherein 
alterations  to  the  structural  response  will  be  uost  beneficial,  (2)  the 
extent  to  which  additional  energy  absorption  is  needed,  and  (3)  the  struc¬ 
tural  element  load-deflection  characteristics  (and,  consequently,  structure 
design  and  size  requirements)  that,  are  needed  to  meet  a  specified  or  desired 
crashworthiness  capability. 

As  a  first  step  in  establishing  the  mathematical  model,  it  is  necessary  to 
determine  the  types  of  failure  modes  for  the  design  involved.  There  are 
several  potential  types  of  failures  including  tension,  compression,  bending 
and  shear,  as  well  as  combined  failures,  e.g.,  axial  and  bending.  Tables  II 
and  III  show  some  of  the  typical  structural  elements  used  in  aircraft  along 
with  their  failure  modes.  For  example,  the  lower  fuselage  of  the  UH-1H  is 
composed  of  a  series  of  riveted  interconnecting  bulkheads  and  panels  with 
top  (floor)  and  bottom  (outer  contour)  skins.  The  panels  throughout  the 
fuselage  have  a  variety  of  different  sized  angle  stiffeners  attached  to 
them.  The  depth  of  the  structure  and  the  spacing  of  the  stiffeners  vary 
depending  on  their  function  in  the  airframe  of  the  vehicle.  Although  it  is 
desirable  to  represent  the  entire  fuselage  with  finely  divided  segments  to 
accurately  describe  all  the  forces,  acting  on  each  element ,  it  is  doubtful 
that  sufficient  detail  to  do  so  would  be  available  during  the  preliminary 
design  stage.  Even  if  detailed  information  were  provided,  the  task  of  es¬ 
tablishing  such  a  rigorous  representation  of  the  entire  vehicle  would  be 
unwarranted  from  practical  and  economic  considerations ,  especially  since 
the  adequate  representation  of  the  crushing  behavior  of  the  fuselage  struc¬ 
ture  is  the  primary  objective. 

To  establish  a  mathematical  representation  of  an  aircraft  for  crash  analysis 
requires  an  awareness  of  the  structures'  behavior.  Up  to  the  point  of 
failure,  the  anticipated  behavior  of  the  structure  is  expected  to  be  defin¬ 
able  since  the  design  for  normal  operation  will  be  such  as  to  preclude 
reaching  failure.  To  model  for  a  dynamic  impact  such  as  a  crash  condition 
requires  that  the  mass  characteristics,  as  well  as  the  stiffness  charac¬ 
teristics,  of  the  structure  be  represented.  Regardless  of  the  analytical 
methods  employed  in  treating  a  complex  problem,  some  basic  assumptions  are 
always  required  which  compromise  the  actual  physical  system,  The  extent 
to  which  the  physical  system  is  compromised  depends  on  many  factors 
including:  available  detail  design  data,  accuracy  required,  available 
computer  systems,  economic  and  time  considerations.  The  above  noted  factors 
will  influence  the  degree  of  refinement  to  be  used  in  mathematically  repre¬ 
senting  the  structure. 

Eeam  model  representations  of  the  major  structures  with  lumped  masses 
adopted  in  program  KRASH  normally  result  in  a  more  coarse  modeling  of 
a  structure  than  would  be  obtained  using  a  more  detailed  finite  element 
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technique,  although  both  approaches  are  similar.  As  described  herein  and 
in  Reference  1  program  KRASH  emphasizes  crashworthiness  analysis  during 
preliminary  design.  Thus,  the  application  of  an  analytical  model  con¬ 
sisting  of  beam  elements  and  lumped  masses,  to  be  of  benefit  to  the  de¬ 
signer,  must  represent  the  gross  behavior  of  the  structure  with  reasonable 
accuracy.  The  representation  of  the  load-deflection  characteristics  depends 
on  the  type  of  structural  element,  its  location  in  the  structure,  and  the 
manner  in  which  it  interfaces  with  other  structural  elements.  For  example, 
in  the  region  of  the  UH-1H  lower  fuselage,  the  overall  crushing  deformation 
is  primarily  influenced  by  the  static  structural  modes  and,  therefore,  sub¬ 
structure  sections  may  be  represented  by  their  static  crushing  character¬ 
istics.  However,  as  is  shown  by  the  results  of  the  study,  the  total  load- 
deflection  characteristics  of  a  large  subsection  can  be  obtained  with  rea¬ 
sonable  accuracy  from  piecewise  superposition  of  the  individual  subelement 
load-deflection  curves. 

Those  mass  items  whose  response  is  influenced  by  the  energy  absorbing  capa¬ 
bility  of  the  lower  fuselage,  the  location  of  the  masses  in  the  fuselage, 
the  manner  in  which  the  mass  items  are  attached  to  the  fuselage  along  with 
the  type  and  size  of  the  attachment  members,  influence  the  load-deflection 
characteristics.  As  a  consequence,  the  load-deflection  curve  may  have  to 
be  obtained  in  a  different  manner  than  by  piecewise  superposition  of 
individual  subelement  stiffnesses,  For  the  mass  items  such  as  engine  and 
transmission  the  subelement  stiffness  is,  in  effect,  the  element  stiffness. 
Although  each  condition  has  to  be  reviewed  and  evaluated  on  its  own  merits, 
simplified  load-deflection  representations  sure  possible  without  knowing  the 
exact  post-failure  characteristics  of  the  entire  structure. 

The  results  of  the  load  sensitivity  studies  indicate  that  the  post-failure 
curve  can,  depending  on  the  type  and  location  of  the  particular  structured, 
element,  be  modeled  in  em  approximate  manner.  In  particular ,  several  load- 
deflection  categories  can  be  established  which  will  yield  response  values 
of  acceptable  accuracy  for  preliminary  design  purposes.  In  the  literature 
review  and  evaluation  it  was  noted  that  structural  elements  can  be  classi¬ 
fied  into  load-deflection  categories.  The  different  post-failure  load- 
deflection  categories  are  shown  nondimensionalized  in  Figure  104.  Included 
for  each  type  of  post-failure  curve  is  the  corresponding  stiffness  reduc¬ 
tion  factor  (KR)  and  distance  (X)  terms  which  can  be  used  to  model  the 
representative  beam  element  for  the  respective  structural  elements  in  pro¬ 
gram  KRASH. 

The  three  categories  shown  in  Figure  104  are  described  as  follows. 

I  -  The  load-carrying  capability  increases  as  the  deflection  increases. 
This  category  is  applicable  to  locations  where  secondary  structure 
is  available  to  carry  loads  after  the  initial  structure  fails; 
e.g. ,  a  soft  transmission  mount  on  the  pylon  support  structure. 
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(a)  Category  I 


(b)  Category  II 


Figure  104.  Simplified  Load-Deflection  Categories. 
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II 


-  The  load-carrying  capability  remains  essentially  constant  as 
the  deflection  increases  beyond  the  yield  load.  This  category 
is  applicable  to  energy  absorbers  like  Dynasorb  or  honeycomb 
in  compression;  e.g.,  seat  systems  in  which  load  limiters  are 
used. 

Ill  -  The  load-carrying  capability  decreases  beyond  the  failure  load. 
This  category  is  applicable  to  support  struts  for  some  major 
mass  elements 'such  as  an  engine  on  pin  ended  column  supports. 

Note  that  in  Figure  104  the  terms  represent  the  deflection  at  initial 
failure  or  yield  and  are  obtained  from  the  pre-failure  mode  which,  in  general, 
is  elastic.  A  KR  of  1  means  that  a  linear  stiffness  applies  to  the  struc¬ 
tural  element  until  deflection  is  reached.  The  program  is  coded  with 
sufficient  generality  to  handle  any  practical  type  of  load-deflection  curve; 
thus,  in  order  to  represent  a  transition  from  one  stiffness  to  another  stiff 
ness,  a  small  incremental  change  in  deflection  is  a  required  input.  This 
small  incremental  change  is  expressed  as  Xg  -  X-^  for  each  of  the  categories 
in  Figure  104.  A  KR  value"  of  0  provides  a  zero  stiffness  for  the  associated 
increment  of  deflection  whi< h  also  results  in  a  constant  force.  Negative 
KR’s  result  in  decreasing  loads  for  the  respective  increments  of  deflection. 
For  the  category  III  type  elements,  two  representations  are  shown.  The 
solid  line  represents  a  failure  which  results  in  a  rapid  loss  of  load¬ 
carrying  capability  in  the  post-failure  region.  The  dashed  line  is  a  more 
gradual  loss  of  load  capability.  From  the  UH-1H  engine  sensitivity  given 
in  Tables  XI  and  XII  for  a  42  ft/sec  vertical  crash  condition,  the  differenc 
in  load  between  the  dashed  and  solid  lines  would  be  approximately  12%. 

DETAIL  APFROACH 

Detail  Design  Procedure  and  Illustrative  Sample  Calculations 

The  following  procedure  is  a  description  of  the  step-by-step  process  that 
a  designer  would  follow  to  obtain  the  load-deflection  characteristics  of 
a  composite  built-up  substructure  made  of  stiffened  panels  with  main  beam 
longitudinal  reinforcements  along  the  sides.  The  procedure  is  based  on  the 
simplified  analytical  techniques  described  earlier  in  the  report  and  illus¬ 
trates  the  types  of  curves,  nomographs,  tables  and  equations  which  can  be 
used  to  formulate  a  designer-oriented  procedure.  At  present,  the  simpli¬ 
fied  analytical  techniques  have  been  verified  by  test  for  four-edge  and 
two-line  supported  configurations.  An  example  is  included.  The  analysis 
presented  in  this  section  is  based  on  the  representation  of  the  actual 
structure  in  the  region  from  Fuselage  Station  (FS)  102  to  180,  Water  line 
(WL)  10  to  22,  Butt  Line  (BL)  -l4  to  +  14.  The  angle  stiffener  and  bulk¬ 
head  arrangement  is  shown  in  Figure  105. 


V'igure  105.  Structure  U8ed  in  Example. 


The  pertinent  dimensions  for  the  structure  are: 

Bulkheads;  t  *  0.032  in. 
s 

t  =  0.032  in. 
w 

angle  stiffener;  0.75  in.  x  0.75  in.  x  0.032  in. 

Main  Beams;  1  in.  x  0.063  in.  x  78  in. 

Side  Webs;  t  =  0.032  in. 
s 

Side  Stiffeners;  t  =  0.032  in. 
v 

Fillers;  1  in.  x  0.063  in.  x  8.8  in, 

Material;  2C2U  -  T3,<r  =  *)5>000  psi 

cy 

A  sketch  of  the  main  beam,  stiffener  web,  side  stiffener  and  filler  is 
shown  on  the  next  page.  Dimensions  in  figures  and  sketches  in  this 
section  are  in  inches  unless  otherwise  noted. 


MAIN  BEAM 


SECTION  A-A  (from  Figure  105) 

The  assumed  deformation  mode,  as  described  in  the  analysis,  for  the  four- 
edge  support  condition  is  shown  in  Figure  106. 


8  =  VERTICAL  CRUSHING  DISTANCE 
O  ~  PLASTIC  HINGE  FORMATION 


Figure  106.  Assumed  Deformation  Mode. 


The  three  bulkheads  at  stations  C,  D  and  E  (Figure  106)  and  the  side 
stiffener  webs  in  this  region  are  subjected  to  the  in-plane  compression. 
The  side  webs  between  B  and  C,  and  between  E  and  F  are  under  a  varying 
compressive  load.  The  main  beam  deformation  mode  is  shown  in  Figure  106. 
This  is  the  case  where  both  beams  (upper  and  lower  beams  as  shown  in 
Figure  105) behave  as  individual  b'ams  restrained  by  the  web. 
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The  step  by  step  procedure  is  as  follows: 
A.  Bulkhead  Analysis 


The  bulkhead  is  represented  as  shown  below : 


*28- 


1..  Failure  Load  Analysis 

a.  Examine  the  column  instability  effect. 

For  stiffened  panel  determine  the  effective  radius  of  gyration 
(p)  as  follows  for  each  bulkhead  and  side  panel. 

(l)  determine  the  effective  neutral  axis: 


NEUTRAL 

r~~j 

_ 1 

AXIS 

n  A  X, 

Neutral  Axis  N/A  =  £  --j . —  (56) 

i=l  Ai 

where 

=  area  of  segment  i 

X.  =  distance  to  c.g.  of  segment  i  from  outside 
face  of  the  web 


(2)  determine  the  area  moment  of  inertia  (I) 

3 

.+  A. 


n 

£ 

i*l 


Vi. 


12 


W 

- 


(57) 


where 


=  width  of  segment  i 
=  thickness  or  height  of  segment  i 
d.^  =  distance  from  c.g,  to  N/A  for  segment  i 


(3)  sum  up  each  area  segnenr  to  obtain  total  area  (A) 

(4)  calculate  radius  of  gyration 

p  =  >Jl/h  (58) 

(5)  determine  ratio  of  length  to  radius  of  gyration 
(L'/p),  where  L'  is  the  effective  column  depth. 

If  L'/p <30  then  continue  with  section  A  through 
D  described  herein.  If  L'/p  >3°  then  modify 
calculations  for  long  column  effects  (subsection 
E  of  the  Detail  Approach). 

Following  the  above  procedure  a  determination  is  now  made  to  see  if  column 
effects  have  to  be  included ; 


t 

e  =  0.0620431  s  —  e  -  “■  =  0.0460431 


.752  +  .718  (.032) 

eo  =  1/2  .75  +  .718 

e  =  0.1994  in. 
o 

I  =1/3  [Be3  -b  (e  -t  )3  +  t  (H-e  )3] 

0  '  o  o  w  w  '  0  j 

-1/3  .0059461  -  .0033681  +  .0053414]  =  2.6398  X  10“3  in.1* 
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Therefore: 


It  ■  IQ  +  1/12  | 5.4  (.032)3J+  .1728  (.0U60U31)J 

+  .0^698  ( ,1994  +  .032  -  .0620U31)2 
It  =  4.3682  x  10‘3 


=  -14098  in* 


and  taking  end  fixity  as  C  =  3 


P'  ”(Vc)(P)"  V3  (.14098) 


=  1+9.14 


Therefore,  for  the  bulkhead,  correction  for  long  column  effects  is 
required. 

b.  Determination  of  critical  buckling  stress  (a  ) 

op  cr 

k  tt  e  ft  r 


(1)  a  =—777 
cr  12(1 


TT  E  ft 

" . — “K-  — 

lb  I 


t  =  skin  thickness 
s 

b  =  skin  width  between  angle  stiffeners 
8 


2  rt 

=  KJS.  "her.  1^  .  [^j 


For  aluminum,^  can  be  calculated  from  the  above  expression  or  obtained  from 
Figure  107  knowing  t 


Kc  can  be  obtained  from  Figure  108. 

For  the  bulkhead  in  the  illustrative  procedure  the  buckling  load  for  the 
sheet  skin  of  one  angle  stiffener  spacing  is  obtained  as  follows: 


a  =  12  in. 

b  =  5.4  in. 

s 


a/b  = 


=  2.2 


boundary  condition  A  (all  edges  clamped) 
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c  =  clamped 
ss  =  simply  supported 
f  =  free 


Figure  108.  Compressive -Buckling  Coefficient  for  Flat 
Rectangular  Plates  (Reference  12). 
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Therefore,  from  Figure  107  for  condition  A  and  a/b  =  2.2 


X  =  7.8 
c 


_  _  v  — -  a _ 

or  c  12  (l-v2)  b 


=  7.8  =  2476  psl 

\12  (1-.32)/  \  5.4/ 

or  obtain  from  Figure  107  and  multiply  by  Kc 
c.  Determination  of  type  of  failure 

(a)  Wrinkling 

(b)  Inter rivet 

(l)  Calculate  wrinkling  coefficient 

A.)* 

*2T\  T)  E  U8  / 

where  wrinkling  stress  can  be  obtained  as  follows: 

rt  1  rt  iV^rt 
-  .48  Ti  fi  s  f-  f 

w  J  1  W  J  l  S  J 


Therefore 


t  \  V3  /t  ,  1/8  /b  .  3/2 


■--ft)  ft)  ft) 


k  1/2  =  yir 

w  v  w 


(62) 

(63) 

(64) 

(65) 


i 

\ 


For  the  example  structure  K  is  obtained  either  from  Figure  109  or  calculated 
as  follows;  w 


Kw  =  17.577 
therefore 

K^1/2  =  4.1925 

Typical  values  are  as  follows: 

t  ;  from  .03  in.  to  .12  in. 
w 

t  ;  from  .02  in.  to  .09  in. 
s 

b  ;  from  3  in.  to  6  in. 
s 

b  ;  from  .5  in.  to  1  in. 
w 

f  ;  from  .25  in.  to  .75  in. 
w 
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The  next  step  requires  a  check  to  see  'whether  failure  is  due  to  wrinkling 
or  monolithic  behavior  and  not  inter rivet  failure.  The  following 
expression 


(66) 


must  be  satisfied  if  the  failure  mode  is  due  to  wrinkling  or  monolithic 
behavior. 

For  the  example  the  calculations  are  as  follows: 

p  =  1  in. (rivet  pitch) 

b  =  5.4  in.  (angle  spacing) 
s 

"b"  ”  fj =  -185185 
s 


hSzL-  =  ~'KL-  _  3Q2Q2 

|K  jt 72  4.1925  39 


P_  <  1.27 


K 


1/2 


w 


Therefore,  bulkhead  fails  by  wrinkling  or  monolithic  behavior  and  not  inter¬ 
rivet  failure. 

If  failure  is  of  the  inter  rivet  type  then  the  monolithic  failure  (o^.) 

determined  in  the  next  step  would  have  to  be  reduced  to  o  . 

cr 

d.  Determination  of  yield  stress 


' '  I 
**  x  1 


*2 


.032 


1 


_ Ic 


/ 


b2  =  .718 

=  .032 


( 


YZZZZZZ2ZZZZZZZZ2L 
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(67) 


°cy  ” 


°cys  As  +«cyw  Aw 

______ 


cy 


(~w  -  1) 
°-cys  ^cyw  \ 

*A 


n 

i=l  A 

where  t  -  - 

w  n 

Ebi 

i=l 


From  the  sketch  above 


-  _  £bj  ti  _  bi  \  *  b2  t2  +  b3  t3 

W  2bi  bl  +  b2  +  b3 

5- U( .032)  *  .7l8(.032)  +  .75(.032) 

57TT7?IF‘+TT5 


t  =  .032 

w 


For  the  same  material  nr  a 


(68) 


(69) 

(70) 


(71) 


e.  Determination  of  failure  stress  (a^) 


(1) 


For  angle 
follows: 


stiffeners  failure  stress  can  be  determined  as 


(72) 


3  =  .66$  for  angle  stiffeners  (obtained  from  Table  XVI 

and  Equation  (36)) 

(2)  According  to  the  data  in  Reference  12,  the  failure  strength 
can  be  generalized  by  the  following  form: 


(73) 


wher* 

g  =  the  number  of  flanges  +  the  number  of  cuts 
/c\m 

0=0  |-~)  for  multi-corner  elements 

g  c  \g/ 

m 

0  =  0.  (—)  for  one-comer  elements 
g  f  \g; 

m 

0^.  =  0  j  for  two-corner  elements  or  plates 


(74) 

(75) 

(76) 


In  Figure  33  representative  types  of  stiffeners  with  from 
one  to  six  comers  are  shown.  In  Figure  34  and  Table  XVI 
values  of  0^,  g,  n  and  m,  obtained  from  Reference  12,  for 

different  configurations  are  listed. 

A  is  the  cross-sectional  area  indicated  in  the  sketch  below. 


Determine  number  of  flanges,  f,  from  Figure  24. 

Use  cladding  factor  =  1  (fully  cladded);  Reference  12 
notes  other  cladding  factors. 

(3)  Correct  af  for  rivet  effect  (Figure  22  or  23) 

a  =  a  times  the  value  obtained  from  Figure  22 
or  Figure  23 

For  the  example  structure  the  monolithic  failure  stress  for  the 
bulkhead  is  determined  as  follows: 
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cy 


n^(-) 

A  cy 


T1Ev3/2]  0.85 


Fran  Table  XVI,  g  =  5,  P  =  0.55 
A  =  As  +  Aw  =  5#1*  (°-°32)  +  1.468  (0.032)  =  0.219776 


ts  =  0.032,  tw  =  0.032,  ay  *=  45,000  psi,  n  *  1 


af 

-  0.55  . 

'5  x  0.0322  - 

i'1000\  ’  d] 

cr 

cy 

.  0.219776  ' 

i  4.5/  J 

=  0.223843 


To  include  rivet  strength  effect, 
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Frcm  Figure  22 
Therefore: 


(medium  value) 


—■  =  .179  °fr=  8058  psi 

cy 

The  failure  stress  (<rfr)  is  corrected  for  the  long  column  effect  from 
Equation  (29)  using  the  following  values: 

<rc  =  59495  psi  (Equation  (28)) 

<r20  =  24674  psi  (Equation  30)) 

<rcr  =  2476  psi  (Equation  (6l)) 

<rfr  =  8058  psi  (Equation  (72)  and  Figure  22) 

From  Equation  (29) 

crco  =  0.979  <rfr  =  7888  psi 

°fr  ttan  be  U3e(i  as  a  conservative  upper  value  if  the  column  effect  is  ignored, 
f.  Determination  of  failure  loed  (pf)  from  failure  stress  (vfr) 


where  n  =  .3^ 
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a  is  obtained  in  Step  2 
cr 

a  is  obtained  in  Step  k 
cy 

o^r  is  obtained  in  Step  5 
(2)  From  Figure  25  obtain 

fe) 


— i  wei|—  -j  “i  (-".a-]  h*H*eih__i!r 


(3) 

(4) 


Knowing  b^  and  K  then  wg^  can  be  obtained  from  the 
expression  above  in  which  (w^/bjj*  K 

Sum  overall  effective  widths 

i~l 


(5)  Multiply  Zw  x  t  =  ZAo 
c  s  s 


(6)  Xaw  =  Awl 


+ 


(7)  Failure  load  =  (Za  + Za  )  o_  per  panel 

S  W  I 

e 

For  the  example  structure  the  effective  width  is 
given  by 


=  .25 


•37( -179) 


. 30212 


C,JMT  T 


N  =  Number  of  effective  widths 
e 

A  =  (N  W  +  2  We_)  t 
se  e  e  e2  s 

-  0°  (.30212)  5.4  +2  (.5)  .30212]  .032 

=  .53173 

=  Number  of  Angles 

Aw  =  NaAfl  =  4( -75  +  .718)  .032  =  .1879 

F.  =  0-  (Aee  +  A  )  =  5800  lb 
f  fr  se  w 

2.  Post -Failure  Analysis 


a.  Determination  of  plastic  hinge  moments 

(1)  obtain  N/A  =  — 


(2)  obtain  M 


-  2  A,  e, 


i  'll  cy 


(3)  calculate  M  as  a  function  of  deflection,  MminJ 

plastic  hinge  rotation  9,  and  constant  K. 

For  the  sample  structure  the  maximum  and  minimum  plastic  hinge 
moments  for  post-failure  analysis  are  obtained  as  fellows: 


M0,  min  =  M_,  max  =  l/k  bt  *  a  =  l/U  (28)(  .032)c(4?000) 
d  d  s  cy 

_ «  322.6  in. -lb _ 

*Mp  for  one  solid  rectangular  cross-section  =  jj-  bt?‘'ocy 


e 


The  neutral  axis  N/A  for  the  stiffener  and  skin  is: 
e  =  .062 

M,,  max  *[28(.032)(.062-.0l6)  +  4(  .75)  (.032)  (.062  -  .032  -.016) 

+  .7l8  ^3^.  +  .032+. 032 062 j  (.032)  4]ocy  =  .0757373(45000)  =  3408  in. -lb 

To  allow  for  the  effect  of  brake  forming  on  the  tensile  and  compres¬ 
sive  stress  as  noted  in  Reference  105  a  1 %  increase  in  the  plastic 
hinge  moment  is  introduced. 

•U*1-15  WJ 

=  4290  in. -lb 

For  we  assume  the  flange  of  the  angle  has  folded  as 


m  -  f?8( .032) (.02626-. 016)  +4  ( .75+. 718) (.032) (.032+. 032 

I,  Min  *- 

-  .02626)3  45000  =  733  in. -lb 
Therefore 

\in  “  1*15  (322  +  733 }  =  1213  in‘"lb 
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The  post-failure  plastic  hinge  moment  is  given  as: 

M  -  \in  *  *  MW  '  ■K9 

Therefore, 

M  *  1213  +  3077  e'K® 


where 

K  .  5  (f)2  (fy  4.2 

9  “  t“"1  [4^T/(l-t)] 

and  the  post-failure  load-deflection  cur/e  can  be  obtained  from  the 
following  expression: 


1213  +  3077  e'U*23 


where 

£  a  12 

u  =  vertical  crushing  distance 
B.  Side  Stiffener  Analysis 


The  side  stiffener  web  is  represented  as  shown  below: 


1.  Failure  Load  Analysis 


a.  Examination  of  column  instability  effect. 


Determine  if  column 
°'*fr 


J.  U  ouuuxu 


1=2^.  (.032)3  =  2.048x10 
s  12 


-5 


e  = 
0 


I  = 
w 


1  *  (1,T.!.Q^L^  =  .262 
2(2-.032T 


21  =  1.24532  x  10 
w 


-2 


I  *  I  +  21*  +  .06298  (.262  +  .048)2  (2)  +  21 

S  x  W 

+  2  (1)  (.032)  (.032)2  =2.465  x  l0"2in.4 


A  =  Ag  +  2Aw  +  2Af  =  7.5  (.032)  +  2  (I.968)  .032 
+  2  (1)  x  .032  =  .43  in.2 


.2394  in. 


12 

*f5  (.2394) 


28.9  <30.0 


From  the  flat  curve  range  shown  in  Figure  29  the  column  instability  effect 
is  not  considered  significant.  Therefore,  any  further  consideration  of 
column  effect  on  is  ignored  in  the  subsequent  analysis. 

b.  Determination  of  critical  buckling  stres? 


a/b  =  12/7.5  =1.6 

s 

k  =  8.5  from 
c  Figure  100 
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2  „ 

a  =  k  JL- B — _ 
°c r  c  12  (l-v^) 


t 

b 

L  s 


(86) 


=  8.5 


2  7 

TT  X  10 ' 


12  (1-0. 32) 


0.032 

L  7.5 

=  1398.5  psi 


c«  Determination  of  type  of  failure 

For  this  particular  configuration  the  side  stiffener  panel 
consists  of  five  double  angle  stiffeners  between  the  supports. 
This  is  considerably  stiffer  than  the  bulkheads  in  the  same 
region.  Therefore,  inter rivet  failure  is  ruled  out  in  this 
case.  However,  in  general,  the  following  expression 
must  oe  evaluated : 


t< 


1.27 


s  K 


1/2 


(87) 


w 


as  previously  performed  for  the  bulkhead, 
d.  Determine  Yield  Stress 
Due  to  same  material 


a  -a  =  a 
cy  cys  cyw 

e.  Determine  the  Failure  Stress 
The  monolithic  failure  stress  is: 


From  Table  XVI, 


(88) 


(89) 


Pg  =  .0583,  g  =  7 

tw  =  2  x  t  (due  to  double  stiffener) 
=  0.064 
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Av  =  1.968  (0.032)  =  0.062976 
A  =  7.5  (0.032)  =  0.24 

O 


A=A  +  2A  =  0.365952}  a  =  45000  psi 
s  w  y 


’ey 


«  0.583 


(0.064)  (0.032) 

0.365952 


=  0.3690681 

To  account  for  the  rivet  strength  effect 


»  67 


Therefore, 


-^=  0.33216  °fr«  14,947  psi 
cy 


f.  Determination  of  failure  load. 

In  calculating  p  previously  the  filler  is  included  in  both 
the  A  and  I  values.  This  is  necessary  because  the  fillers 
are  assumed  ?o  participate  in  the  bending  stiffness  property. 
However,  in  the  calculation  of  crippling  load  (or  failure 
load)  the  fillers  will  not  carry  the  direct  compressive  load. 
Therefore,  their  contribution  is  neglected  in  determining 
the  failure  load. 


2A  =  2(5)  A  -  10  (.06298)  =  .6298 
w  v 

Agfi  =  10  (7.5)  .032  =  .49039  f*e 

Pf  =  °fr  (  ~Aw  ♦Age)  =  i6»T43  lb* 


obtained  from 
Equation  (77) 


The  plastic  hinges  are  assumed  to  occur  at  Sections  AA  and  BB 


The  schematic  model  of  symmetric  plastic  hinges  is  shown 
above. 


Therefore: 

Mmax  ~  jMl’  max  +  M2’  max]  1.15 

=  2  (1A)  b  tg2  +  20  Ay  Sy  *  10  Af  Sf]  <ry 

=  [2  (1A)(1*2).0322  +  20  (.06298M.262  +.032  +  .016) 

+  10  ( .032)  ( .032)  ]  cr  (1.15) 

••  *  » 

=  A222  (1*5000)  (1.15)  -  21850  in. -lb 


(92) 


For  w  the  cross-section  has  the  following  geometry ; 


\  min  *  [V1*  11  's'  +  10  Av  %  +  10  Af  S'f]  K) 

M2,mln'[1/1*>>t32tl0AvS;)(-y) 

Mmin  =  |2  UA)(42)(G.0322)  +  20  (0.062976) (0.111) 
+  10  (0.063)  O.OU75]  (<ryy 
=  0.19123572  (  45000)=  8605.6  lb-in. 


for  modified  plastic  moment 


K  =  5 


/  20  \2  /0.032\  „  , 
\23.57/  \O.O32j 

M  =  8605.6  +  13036  *  e"3,6e 


where 


0  =  tail 
!■  =  10  in. 


The  length  is  based  on  the  distance  between  section  BB  rather  than  the 
total  length  (12  in.)  because  the  main-beam  flanges  are  assumed  not  to 
bend  laterally.  The  post-failure  load/deflection  curve  for  the  side 
web  is  given: 


/§  d  -  u/2) 


(93) 


C.  Main  Beam  and  Bottom  Shin  Analysis 

Perform  calculations  for  each  side.  The  appropriate  model  is 
shown  in  Figure  28. 


m 


1 .  Load-Def lection  Curve 

a.  Determination  of  neutral  axis  (n/A) 

n  . 

b.  Determination  of  M  =  £  (A.  e.|  <r  for  section  AA  (9*0 

P  i=1\  i  V  cy 


c.  Determination  of  I  for  Section  AA 

-  T  +  y  i  ^ 

‘  1  12  ^  Ad 

d.  Computation  of  P  ,  from  Figure  110 

6  ,  62  from  Figure  111 


or  from  expressions 

pi  =  #mp 


a  -  *3  M£‘ 

61  30  El,  P 


p2  =  7mp 

M  l 2 
62  =  2PEI 


(95) 


(96) 

(97) 


For  the  sample  structure  side  stiffener  web  the  following  sketches  and 
calculations  pertain. 


UPPER  MAIN  BEAM 
LOWER  MAIN  BEAM 

I, 


B  C 

~ T — ! — ! — 1 

P=30*0; 

( 

A 

r 

E  F  / 

G 

0 - s - 1 

2 

- A — “D 

■30' 


4-  12' 


■12 


The  lumped  beam  model  representation  of  the  lower  main  beam  is  given  above; 
where  due  to  rigid  support  the  center  section  of  the  oe am  is  assumed  to 
have  infinite  bending  rigidity,  or  Irt/l_  =  «= . 

The  first  stage  of  hinges  forms  at  C  and  E 

16m  16m 


I 


4 


ft 
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PLASTIC  HINGE  MOMENT  (M  KiIN. 


where  £=  12 


.  13  MlL  ,  oli  A-6  M 

61  =  30E,  p  -  6'24  x  10  -E- 


and  the  maximum  load  and  the  corresponding  displacements  are 

J2!p  i,  i 

p2=~  "15^ 

where  £  =  12. 


2  ~  2EI, 


=  7.2  x  10“6  -2- 


Calculate  M  and  I.,  for  the  beam. 
P  1 

Idealized  beam  cross-section  is: 


MAIN  BEAM 


BOTTOM  SKIN 


T 

a  »  0.126— ^  -f*— 

JL _ HU 


1.0  H 


•2.0* 

TI 

B 

H 

r-b/2 

1 

2  2 
aH  +  bd 

2 

aH  +  bd 

1 

Be  3  -  b(< 

3 

L  0 

aH 

+  bd 

b  s  B-a 


-0.126 , 


0.032  -t  t-o.063 


(99) 

(100) 


=  (0.126)(l.)  +  (2-0.126) (0.032)  =  0.l86l  in.‘ 
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£ 


i 


1  T  (0.126)12  +  (1.874)  +  (0.032)2] 

e  =  —  L - — J  =  0. 

'-/N  ^  A/*. 


■o  2 


0.1061 


343632  in. 


I  =—  [2(0. 343632) 3  -  1.874  (0.343632  -  0.032)- 

o  3  l 

+  0.126  (1  -  0.3l*3632)^j 
=  0.059707  in.1* 

Ag  =  28  x  0.032  =  O.896  in.2 

I  =  —  x  28  x  0.0323  =  7.61*6  x  10“5  in.1* 
s  12 


(28)  (0.032)  (0.01b)  +  (2)  (0.1861)  (0.31*3632  +  0.032) 
28(0.032)  +  (2)(0.l86l) 

=  0.10236  in. 


L.  =  I  +  21  +  O.896  (0.10236  -  0.016)2 

iso 

+  0.1*88124  (0.343632  +  0.032  -  0.10236)2 
=  0.036452  in.1* 

M  =  a  .  [0.896  (0.10236  -  0.016)  + 

p  y  l 

2(2)(o.o63)(o. 10236  -  0.0635)  + 


s 


I 


/0,0232\  0.9142 

- ’j  (2)  (0.126)  +— - - (2)(0.126; 

=  8656  in-lb. 

therefore 

P .  =  Jt  Mp  =  2308.3  lb 
15 

fi  8656  \ 

5l  =  6-2l*x10  (oT^64?)=  1-1*83  in' 

P0  =-Mp  =  2885.33  1b 
2  3 

B2  =  1.731  in. 
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D .  Total  Load-Deflection  Curve 


The  failure  load  (P_)  obtained  from  section  for  the  bulkhead 
and  side  prnels  is  depicted  below: 


FAILURE  LOAD 


INITIAL  FAILURE 


The  total  load-deflection  curve  is  obtained  for  the  sample 
structure  from  the  calculated  values  of  the  previous  sub¬ 
sections: 

A  -  Bulkhead 

B  -  Side  Stiffener  Web 

C  -  Main  Beam  and  Bottom  Skin  ' 

If  no  pre-failure  stiffness  properties  are  available  or  easily 
obtainable  the  following  procedure  to  superimpose  individual 
sub’element  load-deflection  curves  into  a  composite  load- 
deflection  curve  can  be  adopted. 

(1)  For  the  stiffened  panel  (this  pertains  to  bulkhead 
and  side  stiffener  but  excludes  main  beam)  plot 
the  post-failure  load-deflection  curve. 

(2)  Asymptotically  match  the  maximum  failure  load  as 
shown  in  Figure  112. 

(3)  The  approximate  total  load-deflection  curve  is 
obtained  by  superimposing  the  individual  bulkheads, 
side  stiffener  webs  and  main  beam  curves.  For  the 
sample  problem  presented  this  includes  three  bulk¬ 
heads  (C,  D,  E),  two  side  stiffener  webs  and  the 
main  beam. 


Figures  113  through  115  depict  the  load-deflection  curves  for  the  bulkheads, 
side  stiffener  webs,  beam  with  botxom  akin  individually.  Figure  116  depicts 
the  total  composite  load -deflection  curve  which  combines  the  results  shown 
in  Figures  113,  114,  and  115. 
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APPLIED  LOAD 


****-»  -psWigi*?  >  *• 


rf  h - A 


FAILURE  LOAD 


POST-FAILURE  CURVE 


/v—  OBTAINED 

AFTER  POST- FAILURE  - - - - 

CURVE  AND  FAILURE  LOAD  VALUE 
ARE  ASYMPTOTICALLY  MATCHED 


CRUSHING  DISTANCE 

Figure  112.  Load-Deflection  Curve  Obtained  by  Asymptotic  Matching. 

E.  Long  Column  Effect 

In  the  event  that  a  L’/p  value  >30  is  calculated  for  either 
the  bulkheads  or  side  stiffener  webs,  it  is  suggested  that 
the  column  effect  on  the  failure  strength  of  stiffened  panels 
be  considered.  This  is  done  in  the  manner  outlined  below: 

Determine  p  as  follows: 

(1)  Determine  af  (subsection  A-k  of  the  Detail  Approach) 

(2)  Determine  W  /b^  (Figure  25) 

(3)  Determine  b0 


Figure  115.  Load-Deflection  Curve  for  the  Main  Beam  and  Bottom  Skin,  Sample  Problem 


B&aSSS&wgEOM  - 


(4)  Detemine  f°r  stiffener*/-^  jgiSrj  <101> 


(5)  Calculate 


b  t 

e  s  3  s 

A  ’A  »  0 

stiff  stiff  w 


and  using  Figure  32  obtain  (p/p  ) 


w 


(6)  Multiply  value  for  (P/ )  by  pw  to  obtain  p 

(7)  Calculate  Euler  Stress  a 


Euler  Stress  =  o 


ttLe 


e  (L/pr 

(8)  Obtain  stress  for  L/p  =  20 

“20  =  ^0®  "  2l,6'7l‘0  lb/in'2 

(9)  Obtain  column  stress  a 

co 


(102) 


(103) 


O  =1 
co 


acr\  |  20 


1/2  .  -  1/2 


a  V2  -  a  !/2^ 

a20  cr 


(104) 


(10)  Obtain  ultimate  load  (P^,)  as  in  subsection  A  or  B-lc 


P=cr  (l  A  +  Z  A  ) 
f  co  \  se  w/ 


(105) 


(ll)  Continue  with  post-failure  analysis,  main  beam  analysis 
and  total  load-deflection  curve  analysis. 


Structural  Element  Modeling  Guidelines 

This  subsection  provides  some  guidelines  by  which  designers  can  simplify 
modeling  requirements  and  utilize  approximate  techniques  for  performing 
preliminary  design  crashworthiness  studies.  Although  the  information 
presented  is  directly  applicable  to  the  use  of  program  KRASH,  it  is  also 
general,  and  therefore  can  be  used  for  other  analytical  approaches. 
Included  in  this  subsection  are  discussions  regarding  the  determination 
of  stiffness  and  load-deflection  characteristics  for  d:  x’erent  types  of 
structural  elements.  Design  charts,  curves,  tables  and/or  equations  are 
presented  for: 
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•  columns 

•  honeycomb  plates 

•  frangible  tubes 

This  section  contains  information  applicable  to  a  crashworthiness  design 
manual  in  that  it  r.hov 3  the  type  of  information  and  the  format  in  which 
data  can  be  provided.  The  data  presented  in  the  following  subsection  are 
available  in  current  literature. 


A.  Category  Type  III  Structural  Elements 

Structural  elements  which  exhibit  greatly  reduced  load-carrying  capability 
after  failure  load  is  reached  a^e  characterized  mainly  by  buckling  failures. 
As  described  earlier,  in  both  the  load  sensitivity  study  and  general  design 
approach,  this  type  of  structure  can  be  classified  as  a  Type  III  load- 
deflection  category.  Typical  structures  which  fit  into  this  category  are 
shown  in  Tables  II  and  III  of  this  report.  To  illustrate  design  guide  pro¬ 
cedures  for  obtaining  approximate  load-deflection  curves  for  category  III 
type  structure,  Ihe  compressively  loaded  column  structure  is  used  as  a 
sample  element  in  uhis  subsection. 

1.  Solid  Cross-Section  Column 

Column  failures  can  be  characterized  as  primarily  flexural  buckling  or 
combined  torsional-flexural  buckling.  Thick  walled  and  solid  cross- 
section  columns  tend  to  fail  due  to  pure  f.  'ural  buckling  whereas  thin 
walled  columns  exhibit  failure  as  a  result  of  a  combined  torsioi.il- 
flexural  buckling  mode. 


Of  primary  concern  in  column  buckling  and  post-buckling  analysis  is 
the  intermediate  slenderness  ratio  range  (30  <L'/psl2C'.  Using 
Equations  (U3)  and  (1»8)  in  Reference  6l  for  this  range,  the  following 
expression  for  critical  stress  and  lateral  deflection  in  terms  of 
generalized  nondimens ional  quantities  is  obtained: 


g  (T,y)  =  y 


3/2  (y  "  T)  " 


o  (2T)' 


(106) 


where 


T  =  a/a 

y 

y  =  6/S 


186 


i,'  .*  s:'v"  ,  ■ 

mgc*'-" 

?.?<• 


assess® 


8 


P;' 


at,- 


im 


mV 


m 


a  = 


(L’/P)‘ 


TT2  E 

a  =  compressive  stress 
6  =  lateral  deflection 
S  =  core  radius  of  the*  cross-section 
K  =  eccentricity  factor 


L',  p,  E,  o  have  been  defined  earlier. 

y 


The  critical  stress  and  the  corresponding  deflection  are  obtained  from 
the  following  expressions,  obtained  by  setting  9T/9y  =  0: 


*„  -  *» -T1) 3 

5r»*fc-,«KK 


=  a 


(107) 

(108) 


where 


cr  critical  stress 


cr 


cy 


yield  stress 


Equation  (107)  is  solved  numerically  from  known  values  of  o  and  an 
assumed  value  of  K. 

From  equation  (106)  a  (T,y)-curve  is  obtained  by  linear  perturbation 
as  follows: 

(a)  differentiate  equation  (106)  with  respect  to  y. 

1/2  r  ,  rl/2 


_*g  _ 

dy 


1  (*  -  T)  • 


y  -  K 


|  (|  -  t)  -y  -  K 


(109) 


(b)  differentiate  equation  (106)  with  respect  to  T, 

-1/2 


Sg 

8T 


JL 

4 


(t2  +  i) 


rll(|--T)  -y-K  (T)  3/2(no) 


187 


•J 

4 


'Ur‘9  && 


vfrtertGiz*  - 


For  Sg/aT^O 


(ill) 


By  incremental  solution  the  approximate  post-buckling  (T,y)-curve  for 
a  general  column  can  now  be  obtained  as  follows: 


(a)  T£r  from  equation  (107) 

(b)  ycr  from  equation  (108) 


(c)  Assign  increment  of  Ay  and  obtain  AT  from  equation  (106), 
The  new  T  vs  y  becomes: 


ynew  =  y 


T  =  T 
new 


old 

+  Ay 

(112) 

old 

+  AT 

(113) 

(d)  Ur'ng  the  procedure  described  above  a  typical  (T,y)~curve 
is  obtained  for  a  column  with  the  following  values: 

(L'/p)  =  78.3 

o  =  105  lb/ in. 2 

y 

E  •-  3  x  10T  lb/in.2 


K  =0 

Figure  117  shows  the  stress  ratio  vs.  lateral  deflection,  (T,y)-curve, 
for  the  column  with  the  above  noted  parameter  values. 

To  obtain  axial  load  versus  end  shortening  (deflection)  the  following 
procedure  is  to  be  followed: 

(a)  Obtain  core  radius.  Table  XXIV  shows  the  core  radius  for 
typical  cross-sections. 

(b)  Obtain  actual  length  from  the  reduced  length  as  a  function 
of  end  fixity.  Table  XXV  shows  the  reduced  length  as  a 
function  of  end  fixity. 
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Figure  117 •  Stress  Ratio  Versus  Normalized  Lateral  Deflection.  Curve  for  Axially  Loaded 


Figure  118.  Load  versus  Deflection  Curve  for  Axially  Loaded  Coluam 


TABLE  XXIV.  CORE  RADIUS 


S  =  h/b 


s  =  rA 


TABLE  XXV.  REDUCED  COLUMN  LENGTH  AS  A  FUNCTION  OF  END  FIXITY 


Both  Ends 
Simply- 
Supported 


Fixed-Fixed 


Fixed- Simply- 
Supported 


L’=O.O99  0 

0=1/0.699 


Fixed-Free 


(c)  Obtain  load  (P)  and  end- shortening  deflection  (u)  as 
follows: 

F  =(A)(T)(oy) 

_2  2  2  2  /i  h 

u  =  C  tt  S  y  /k£ 


(11*0 

(115) 


Equation  (115)  is  an  approximate  end- shortening  based  on  the  assumption 
that  the  mode  shape  is  sinusoidal  and  only  the  fundamental  mode  is 
retained.  For  the  same  sample  parameter  used  to  obtain  the  (T,  y)-curve 
(Figure  117)  the  corresponding  load  versus  axial  deflection  (P,u)- 
curve  is  shown  in  Figure  ll8. 

For  Euler  column  range  (L’/p  >80)  another  formula  obtained  from 
Reference  111  is: 

g  (T,  y)  =  y  |3  -  i)  -  y  -  K  -Ut^-i)3  =  o  (116) 


Following  the  same  procedure  outlined  previously  in  this  subsection 


if  K  =  0  Equation  (117)  coincides  with  the  Euler  Buckling  stress 


and 

ycr  =Y~  -  1  -  K/3  (118) 

cr 

^  =  3  (7-  i)  -  y  -  k  3^  -  1)  -  3y-K  (119) 

1t  =  1*“(|-1)2(t+1)  ■-$-  [3(|-1)->r-K]  <120> 

To  obtain  load  versus  axial  deflection  the  procedures  described  pre¬ 
viously  are  to  be  followed. 


For  columns  with  thin-walled  cross-sections,  failure  generally  occurs 
as  a  result  of  a  combined  torsional-flexural  buckling.  References 
100  and  102  discuss  the  theoretical  aspects  of  this  problem.  Reference 
111  presents  design  charts  for  typical  cross-sections  of  thin-walled 
columns.  The  design  charts  provide  preliminary  design  oriented  data 
in  that  the  torsional-buckling  loads  are  presented  in  terms  of  Euler 
buckling  loads,  column  geometry  and  support  boundary  conditions. 

A  typical  curve  is  shown  in  Figure  119  for  column  cross-section  with 
unequal  angle  lengths.  Reference  111  contains  design  charts  for 
other  types  of  cross-sections  including  inwardly  and  outwardly 
lipped  channels,  T  sections  and  plain  channels.  It  is  assumed,  for 
simplicity,  that  the  (T,y)-curve  for  thin-walled  columns  follows 
that  of  solid  columns  by  the  factor  a  shown  in  Figure  119.  It  is 
suggested  that  Figure  119  be  used  in  the  following  manner: 

(a)  determine  Euler  stress  load  (Am  ).  For  solid  columns 
use  Equation  (28),  page  63.  fc 

(b)  determine  t  from  t  ^ —  (121) 

x 
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where  t  ~  angle  thickness 

bx  =  short  angle  length  (see  Figure  119) 


(c)  determine  r/ru  where  2  2m2l  (ape  Figure  119) 

(d)  determine  a  from  Figure  119  using  the  r/m  and  by/bx  values 

(e)  determine  critical  load  (P  )  from  P  =  (P_  ,  )  (o') 

v  '  '  or  cr  '  Euler' 

(f)  determine  axial 'load  versus  deflection  curve  (P  vs.  u) 
following  procedures  outlined  for  solid  cross-section 
columns  and  reduce  the  load  curve  proportional  to  the 
value  of  a. 


Figure  119.  Critical  Load  Factor  for  Thin  Walled 

Column  with  Cross-Section  of  Unequal  Angles. 

B.  Category  Type  II 

The  load  sensitivity  study  and  literature  evaluation  presented 
previously  in  the  report  describe  structural  elements  which  exhibit  a 
constant  load  carrying  capability  after  the  failure  load  is  reached, 
as  a  category  Type  II  element.  Included  in  this  category  are  structure 
types  which  are  used  as  load  limiters  or  energy  absorbers  such  as 
frangible  tubes,  honeycomb  and  other  multiple  cell  structure,  stain¬ 
less  steel  straps,  inversion  and  tel  :-scoping  tubes.  Type  II  elements 
are  prcminantly  used  in  design  regions  wherein  control  (limit)  of  the 
load  level  is  desired.  In  particular,  this  type  of  structural  element 
is  used  where  controlled  deformation  is  required;  e.g.,  seat  systems 
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DEFORMATION 

(a)  STATIC  LOADING  (.15  IN./MIN} 


DEFORMATION 

(b)  DYNAMIC  LOADING  (27.6  FT/SEC) 


Figure  120.  Load  Deformation  Characteristics  for  Compressively 
Loaded  Honeycomb  Panel  (Reference  112). 
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and  landing  gears.  Of  primary  concert!  for  load  and  shock  attenuation 
systems  sire  characteristics  such  as; 


•  Energy  absorption  capacity 

•  Elastic  deformation  per  unit  load 

•  Peak  loads  to  initiate  crushing 

•  Load  on-set  rates 

•  Rebound 


Another  consideration  in  the  selection  of  a  load  limiting  system  is 
the  efficiency  of  the  design  with  regard  to  absorbed  energy  per  pound 
of  material  required.  Table  V  (page  18)  in  this  report  presents  a 
comparative  summary  of  absorption  efficiencies  for  different  kinds  of 
configurations . 


To  illustrate,  design  guide  procedures  for  obtaining  load-deflection 
characteristics  for  Type  II  structural  element  honeycomb  and  frangible 
tube  design  data,  from  the  literature,  are  presented  in  the  following 
subsections : 


1.  Honeycomb 


Extensive  testing  has  been  performed  to  determine  crushing  load 
and  characteristics  for  different  sizes  and  configurations  of 
honeycomb  panels.  Reference  23  provides  test  data  regarding  the 
effects  which  changes  in  specimen  configuration,  cell  geometry, 
test  conditions  and  environment  have  on  the  characteristics  of 
crushable  energy  absorbing  materials.  Since  the  load-deflection 
characteristics  of  a  honeycomb  structure  are  sensitive  to  the 
tolerances  for  foil  thickness,  adhesive  cell  size  and  orientation, 
cross-sectional  area  and  outer  skin  dimensions,  test  data  have  been 
generated  for  a  wide  variety  of  honeycomb  configurations.  Refer¬ 
ence  112  is  a  series  of  technical  service  bulletins  provided  by 
a  major  manufacturer  of  honeycomb  structure.  The  reference 
reports  include  comprehensive  design  data,  design  charts  and 
procedures.  The  reference  material  also  includes  an  excellent 
discussion  on  the  general  crushing  characteristics  of  honeycomb 
under  static  and  dynamic  loads  and  the  mode  of  honeycomb  crushing. 
The  following  figures,  tables  and  discussions  obtained  from 
information  in  Reference  112  are  presented  to  illustrate  the 
type  of  data  that  is  contained  in  the  reference. 

Characteristically,  honeycomb  load-deflection  curves  will  look 
like  those  shown  in  Figure  120(a)  and  120(b)  for  static  and 
dynamic  loading.  The  average  level  of  the  load  in  the  relatively 
uniform  area  is  called  the  crushing  stress  and  generally  varies 
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between  40^  and  6of>  of  the  compressive  strength  of  the  honeycomb, 
depending  on  type.  The  difference  can  be  seen  in  the  sharp  peak 
load  at  yield  for  the  static  load.  For  design  purposes,  the 
crushing  stress  for  the  ccmpressively  leaded  honeycomb  panel  can, 
as  an  approximation,  be  considered  to  decrease  as  a  function  of  the 
loading  angle  as  shown  in  Figure  121.  Figure  122  illustrates 
other  types  of  available  data  provided  in  Reference  112  which 
is  useful  to  the  designer  for  preliminary  crashworthiness  studies. 
Some  basic  formulas  which  will  also  aid  designers  are: 


Stopping  distance  (s) 
.  V2 


V  =  Impact  velocity 
G  =  acceleration  level 
g  =  32.2  ft/sec2 


(122) 


minimum  core  thickness  (t  ) 


c  ~  .7(2gG) 


(123) 


crush  strength  (fcr) 

f  _  WG 
cr  A 


W  =  weight 
A  =  area 


(124) 


kinetic  energy  (KE) 


KE  =  f  AS 
cr 


(125) 


Nomographs  are  provided  in  Reference  112  for  all  the  expressions 
presented  above. 


Detailed  design  and  test  data  is  available  in  References  32-35 • 
Application  areas  of  such  energy  absorbing  devices  include  crash- 
worthy  passenger  seat  systems,  landing  gears,  and  other  load 
limiters. 

The  average  load  suggested  in  Reference  34  is 


P  =  K .  A  +  4  it  f  A  a 
t  c  c  y 


(126) 


in  which  K.  is  the  axial  tube  stress  caused  by  the  resistance  to 

ribbon  bending  and  splitting,  A  is  the  cross-sectional  tube  area, 

f  is  the  friction  coefficient,  A  is  the  cross-sectional  area  of 
c  7  c 
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one  side  of  the  control  ring  and  Oy  is  the  tensile  yield  stress 
of  the  control  ring.  Figure  123  snows  a  sketch  of  frangible 
tube  energy  absorber  showing  the  tube  control  ring  and  splitter 
ring.  For  the  sizes  tested  in  References  32-35 >  K^.  and  f,  are 
given  in  Table  XXVI. 

TABLT:  XXVI.  TYPICAL  VALUES  OF  Kj  A!®  f  (REFERENCE  34) 


Material 

1  f 

c 

2024-T3 

4000 

II 

0.067 

Brass 

7000 

0.156 

Az  -  31B 

1750 

0.110 

Steel  1015 

11000 

0.119 

Steel  4130 

25000 

-i 

o.io4 

A  typical  load-deflection  curve  for  a  tube  made  of  mild  steel 
1015  is  shijwn  in  Figure  124.  From  Equation  (126),  the  average 
load  is 

P  =  (11000)  (tt)  (2)  (0.049)  +  4n. (0.119)  (0.0764)  (55000) 

=  9670.3  (lbs) 

This  value  is  about  a  7.5$  higher  load  than  the  test  obtained 
load  and  overestimates  the  energy  absorption  for  6.2  inch  travel 
by  approximately  %. 

C.  Category  Type  I  Structural  Elements 

As  discussed  previously,  Type  I  elements  are  characterized  by  load- 
deflection  curves  which  show  an  increase  in  load-carrying  capability 
after  initial  failure  is  experienced.  Generally,  structure  of  this 
nature  is  representative  of  combined  responses  from  two  or  more 
elements  with  different  modes  of  failure.  For  example,  the  fuselage 
bumper  structure  tested  and  analyzed  during  the  study  described  in 
Reference  1  illustrates  Type  I  structure.  The  fuselage  bumper  load- 
deflection.  curve  is  shown  in  Figure  1.  The  crushing  of  the  spherical 
cap  occurs  within  the  first  two  inches  of  deflection.  However,  as  the 
cap  collapses,  the  load  is  transmitted  to  two  bulkheads  which  exhibit 
a  second  region  of  increased  stiffness  (2  to  3  inches  deflection). 

As  the  bulkhead  plastically  deforms,  the  stiffness  reduces  and  a 
relatively  constant  load  is  realized.  This  is  shown  in  the  region  in 
which  the  structure  deflects  3  to  9  inches.  Another  example  of  a 
Type  I  loading  situation  is  a  transversely  loaded  beam  on  bulkhead 
supports  at  the  ends.  The  beam  will  yield  locally  in  bending  and  pull 
the  bulkhead  into  tension.  The  load-deflection  curve  for  this  condition 
is  characterized  in  the  following  sketch: 
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Deflection 


The  general  discussion  regarding  this  particular  type  of  structural 
/failure  indicates  that  the  load-deflection  curve,  in  reality,  consists 
of  load-deflection  curves  of  particular  elements,  each  one  of  which 
may  be  of  eicher  the  Type  II  or  Type  III  category.  One  approach  to 
developing  a  load-deflection  curve  for  Type  I  elements  is  to  super¬ 
impose  the  load-deflection  characteristic  of  the  individual  elements 
in  the  sequence  in  which  they  occur.  For  example,  in  the  UH-1H,  the 
transmission  mount  can  be  idealized  by  a  soft  stiffness  region 
representative  of  its  rubber  mount  followed  by  a  muc  higher  stiff¬ 
ness  region  which  represents  the  pylon  support  structure,  the  stiffness 
of  which  is  important  after  the  rubber  mount  has  bottomed. 
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RESULTS  OF  THE  PROGRAM 


The  results  of  the  program  are  summarized  as  follows: 

LITERATURE  SURVEY  ARP  CRITIQUE 

A  literature  survey  and  critique  of  technical  reports  and  publications  are 
completed*  Sixty  reports  are  reviewed.  The  literature  is  discussed  with 
regard  to  eight  areas  of  applicability  to  the  development  of  designer- 
oriented  crashworthiness  design  procedures'.  A  literature  Index  is  provided 
in  Volume  II  which  further  divides  each  report  into  specific  content  and 
area  of  applicability.  There  are  29  such  divisions.  The  basic  contents  of 
the  literature  sure  reviewed  and  commented  upon.  Some  notable  comments  re¬ 
garding  the  literature  include: 

•  The  analytical  procedures  and  empirical  data  which  are  available  in 
the  literature  are  not  oriented  toward  simplified  design  procedures 
because 

1)  they  are  related  to  unique  specimens  or  loading  conditions. 

2)  they  require  the  use  of  complex  computer  programs. 

3)  the  empirical  data  is  not  related  to  analytical  procedures. 

•  Analytical  techniques  and  numerical  methods  of  stability  analysis 
and  post -failure  analysis  are  primarily  concerned  with  the  tech¬ 
nique  rather  than  the  application  of  the  techniques  to  practical 
design  problems. 

LOAD - SEU SI TI VI TY  STUDY 

Program  XRASH  and  an  existing  31  mass  model  of  the  UH-1H,  developed  under  a 
previous  study  contract  with  the  U.  S.  Army  Air  Mobility  and  Research  Lab¬ 
oratory  (U3AAMRDL),  Fort  Eustis,  Virginia  (Reference  l),  is  used  to  deter¬ 
mine  the  sensitivity  of  the  response,  for  several  major  structural  elements, 
to  changes  in  load-deflection  characteri sties  in  the  post-failure  region, 
Variations  were  made  in  the  load-deflection  characteristics  for  the  engine, 
transmission,  landing  gear  and  fuselage.  A  total  of  27  parameter  variation 
cases  are  analyzed. 

The  results  of  the  load- sensitivity  study  using  the  UH-1H  math  model  indi¬ 
cate  that  during  a  vertical  crash  impact  with  a  95th  percentile  potentially 
survivable  accident  (42  ft/sec  vertical  velocity): 

•  The  energy-absorption  capability  of  the  structural  elements  in  ques¬ 
tion  and  their  particular  location  within  the  airframe  will  influence 
the  degree  of  refinement  with  which  they  are  to  be  modeled. 
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•  For  "noncritical  energy-absorbing  structure"  in  the  post-buckling 
region, the  load -deflection  characteristics  do  not  require  a  soph¬ 
isticated  modeling  representation. 

•  The  "noncritical  ene 'gy-absorbing  structure"  for  the  UH-1D/H  air¬ 
craft  modeled  (i.e.,  .Landing  gear,  engine  mount,  transmission  mount 
in  the  UH-1H)  can  be  represented  with  three  or  less  nonlinear 
load-deflection  categories. 

•  The  fuselage  structure  load-deflection  characteristics  including 
the  bottoming  spring  stiffness  exert  a  significant  influence  on 
the  response  of  the  major  upper  masses  (engine  and  transmission) 

as  well  as  the  occupants  due  to  the  fact  that  the  fuselage  crushing 
is  the  primary  absorber  of  energy  in  the  airframe. 

SUBSTRUCTURE  ANALYSIS 

A  general  procedure  for  determining  the  total  load-deflection  characteris¬ 
tics,  including  failure  and  post-failure  behavior,  for  the  structure  se¬ 
lected  for  the  test  program  is  presented.  The  structural  elements  analyzed 
and  tested  in  this  study  are  representative  of  the  multiweb  construction 
of  helicopter  lower  fuselage,  which  consists  of  stiffened  panels  and  angle 
stiffeners  in  the  region  where  a  transmission  pylon  is  supported.  The  gen¬ 
eral  procedure  for  determining  the  load -deflection  characteristics,  in¬ 
cluding  failure  and  pest -failure  behavior,  for  the  structure  includes  a 
step-by-step  process  with  the  following  sequence: 

•  Prediction  of  failure  loads  for  stiffened  panels 

•  Post-fail’  re  analysis  of  stiffened  panels 

•  Main  beam  and  bottom  skin  analysis 

•  Formulation  of  total  substructure  load-deflection  curve 

The  procedure  takes  into  account  monolithic,  winkling  and  interrivet 
failure  modes.  The  basic  equations  are  developed  for  a  nondimensional 
L/p  ratio  of  20.  For  L/p  ratios  greater  than  20,  column  effects  are  also 
included. 

Although  developed  for  angle  type  stiffened  panels,  the  procedure  is  ap¬ 
plicable  to  a  variety  of  panel  types,  including: 

•  T-type  stiffeners 

•  Formed  angle  stiffeners 

e  Extruded  angle  stiffeners 

•  Formed  multicorner  sections 
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•  Formed  hat  stiffened  panels 

•  Extruded  Y  stiffened  panels 

Individual  load-deflection  curves  are  presented  for  specimens  1  through  6 
depicting: 

•  Lateral  bulkhead  crushing 

•  Side  stiffener  web  crushing  (representing  major  bulkheads  at  BL+14 
and  BL-l4) 

•  Side  stiffener  beam  bending 

The  individual  load-deflection  curves  are  superimposed  upon  each  other  to 
obtain  a  total  load-deflection  curve  for  the  total  structure.  Predicted 
total  load-deflection  curves  for  all  the  different  specimens  are  presented. 

The  simplified  analytical  technique  developed  in  this  study  provides  a 
reasonable  representation  of  the  load-deflection  characteristics  of  multi- 
webbed  angle  stiffened  structure  typical  of  helicopter  structure.  While 
the  method  cannot  accurately  predict  the  failure  behavior  of  each  individ¬ 
ual  stiffener  element  it  does  satisfactorily  reproduce  the  overall  crushing 
characteristics  and  general  deformation  behavior  of  the  structural  segment. 
The  assumption  of  a  symmetric  failure  mode  appears  justified  in  a  gross 
sense  although  some  failure  modes  will  be  asymmetrical  in  nature.  While 
asymmetrical  failure  modes  will  introduce  a  small  error  in  the  predicted 
load  compared  to  loads  based  on  the  symmetric  mode  failure  assumption,  the 
representation  of  the  peak  failure  load  will  not  be  impaired  and  the  energy 
absorption  capability  of  the  structure  will  be  sufficiently  accurate  for 
preliminary  design  purposes.  The  primary  interest  in  the  analysis  pre¬ 
sented  herein  is  the  development  of  analytical  tools  for  preliminary  design 
analysis  and,  as  such,  the  emphasis  is  placed  on  adequate  prediction  of 
failure  loads  and  energy  absorption.  The  simplified  analytical  procedure 
does  not  account  for  structural  imperfections.  However,  these  imperfec¬ 
tions  influence  the  pre-failure  region  (pre-  and  post-buckled  behavior) 
more  than  the  post-failure  region  where  the  imperfections  become  small 
relative  to  the  actual  deflection  of  the  structure. 

The  analytical  technique  described  herein  tends  to  predict  a  lower  failure 
deflection  value  than  recorded  from  the  tests.  However, the  energy  absorp¬ 
tion  contribution  to  the  total  crushing  energy  is  extremely  small  and  the 
small  discrepancy  in  this  region  will  not  introduce  any  significant  error 
nor  detract  from  the  analytical  representation  of  the  structure's  crashwor¬ 
thiness  capability.  Furthermore,  the  analytical  approach  presented  in  this 
study  implies  that  an  upper  bound  yield  stress  is  reached  at  failure  which, 
in  principle,  is  equivalent  to  obtaining  the  failure  deflection  from  pre- 
and  post-buckling  stiffness  calculations. 


The  literature  evaluation  indicates  that  the  calculation  of  pre-failure 
stiffnesses  and,  consequently,  deflections  is  a  formidible  task  due  to  the 
effect  of  initial  structural  imperfections. 

The  analytical  method  is  modified  by  the  results  of  the  test  in  the  manner 
in  which  the  post-failure  curve  is  developed.  The  initial  assumption  of  a 
fully  warped  flange  resulting  in  a  conservative  (lower  bound)  plastic  hinge 
moment  is  modified  using  an  empirically  developed  exponential  function 
(Equation  51)  for  angle  stiffened  panels.  The  modification  results  in  im¬ 
proved  agreement  between  test  and  analytically  predicted  failure  deflec¬ 
tion. 

The  analysis  is  intended  to  apply  statically  predicted  load-deflection 
characteristics  to  crash  simulation  studies.  For  lightweight  structure  of 
the  type  used  in  helicopter  fuselage  design,  the  failures  due  to  crippling 
would  most  likely  be  caused  by  a 'lowest  mode,  in  which  case,  a  statically 
predicted  load-deflection  curve,  although  slightly  conservative,  would  suf¬ 
fice  in  lieu  of  a  dynamic  analysis.  The  test  results  show  that  the  dynamic 
tests  yield  slightly  higher  (  ~  5%)  failure  loads  when  compared  to  static 
tests  (specimens  7,  8  and  9)  for  the  type  of  structure  tested  in  this  pro¬ 
gram.  However,  there  are  situations  where  static  failure  modes  do  not 
occur  under  dynamic  loads.  Under  this  circumstance,  the  utlization  of  a 
statically  predicted  load-deflection  curve  for  a  crash  simulation  study 
may  not  provide  sufficiently  accurate  responses.  The  problem  of  structures 
exhibiting  deformation  modes  which  are  different  under  dynamic  loading  is 
discussed  in  References  86,  109  and  110. 

SUBSTRUCTURE  TEST  PROGRAM 


The  substructure  test  program  is  established  such  that: 

•  Typical  fuselage  structure  is  represented, 

•  The  effect  of  detail  design  changes  on  load-deflection  character¬ 
istics  is  investigated. 

•  A  representative  crash  loading  condition  is  reproduced. 

•  The  test  data  can  be  related  to  designer-oriented  analytical 
procedures . 

•  The  test  data  can  be  related  to  analytical  modeling  data  require¬ 
ments. 


The  specimens  selected  for  test  are  representative  of  the  transmission 
pylon  support  structure  in  the  region  bounded  by  BL  +  14  to  BL  -  l4,  WL  9 
to  22  and  FS  .102  to  l80.  Twelve  specimens  were  fabricated  for  the  test 
program.  The  specimen  configurations  ana  test  conditions  include: 
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static  load 

dynamic  intact:  l4  ft/sec  <  velocity  <  30  ft/ sec  5 

four-edge  support 
two -line  support 

skin  web  thickness:  .025  in.  to  .064  in. 
angle  stiffener  thickness:  .032  in.  to  .072  in. 

specimen  depth:  6.125  in.  to' 12. 125  in.  ; 

number  of  angle  stiffeners:  12  to  32 

•J 

specimen  weight:  13.6  lb  to  31*25  lb  I 

effect  of  lightening  holes  in  bulkheads 

A 

The  overall  specimen  dimensions  are:  : 

£ 

width  =  18  in.  ! 

length  =  46  in. 

depth  6.125  in.  or  12.125  in.  j 

The  test  results  indicate  that,  for  these  types  of  structural  elements, 
static  testing  to  determine  load-deflection  characteristics  should  yield  ; 

sufficiently  accurate  results  when  compared  to  dynamic  test  results,  but  1 

in  a  more  economical  manner.  A  static  test  can  provide  the  desired  infor-  ;*■ 

nation  with  regard  to: 

•  Shape  of  the  load-deflection  curve  | 

•  Peak  failure  load  (dynamically  obtained  failure  loads  are  between 
9$  and  24$  higher  than  the  loads  measured  from  static  test  and  the 
loads  from  predicted  static  calculations,  respectively) 

•  Deflection  at  which  failure  load  occurs 

•  Energy  absorption  capability  of  the  structure. 

Figure  125  shows  a  comparison  of  the  load -deflection  response  for  two 
dynamic  test  conditions  and  a  static  test,  all  performed  for  a  similar 
type  specimen. 
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Figure  125.  Comparison  of  Static  and  Dynamic  Test  Results. 

A  dynamic  test'  has  the  advantage  that  the  amount  of  springback  from  the 
maximum  deflected  value  is  immediately  evident.  With  static  test,  the 
structure  will  relax  slowly  (up  to  several  days)  to  its  permanently  de¬ 
formed  position.  The  springback  could  be  as  high  as  k0%  of  the  maximum 
deflected  value.  However,  under  a  static  loading  condition,  the  specimen 
experienced  more  deformation  than  under  a  dynamic  impact,  for  an  equal 
amount  of  input  energy.  Figure  126  depicts  a  possible  percentage  range  of 
structural  springback  as  a  function  of  impact  velocity  for  the  type  of 
specimens  tested.  The  curve  is  based  on  a  limited  number  of  tests. 

The  static  test,  in  addition  to  being  more  economical  to  run  than  a  dynamic 
test,  has  the  advantages  of  easier  and  more  precise  alignment  of  impact 
head  with  specimen,  a  more  rapid  test  setup,  and,  generally,  a  lesser  amount 
of  instrumentation  required. 

If  dynamic  tests  are  to  be  performed,  load  cells  should  be  installed  in  lieu 
of  accelerometers  to  obtain  a  better  definition  of  the  load  time  history. 

Hie  response  from  the  accelerometers  exhibits  oscillations  which  include 
the  frequency  contents  associated  with  the  impact  mass,  structural  specimen 
and  test  support  structure.  The  interpretation  of  the  accelerometer  data 
can  be  difficult  if  the  oscillations  are  severe.  Load  cells,  in  addition 
to  providing  a  smoother  continuous  time  history,  can  be  used  to  indicate 
the  alignment  of  the  intact  her'’  with  the  specimen.  In  the  tests  conducted 
in  this  program,  the  load  cells  were  placed  below  the  support  along  each  of 
the  four  edges,  thus  providing  data  regarding  slight  misalignment  at  impact. 
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Figure  126.  Percentage  of  Structural  Springback  Versus  Impact  Velocity. 

The  typical  response  of  the  type  of  structure  tested  in- this  study  is  a 
relatively  sharp  load  buildup  in  the  elastic  region  up  to  the  peak  failure 
load  at  approximately  .25  inch  of  deflection.  After  failure,  the  load  re¬ 
duces  at  a  rate  less  severe  than  the  initial  buildup  until  the  deflection 
reaches  approximately  .75  to  1.0  inch.  Thereafter,  the  load  remains  con¬ 
stant  until  the  structure  crushes  to  the  point  that  it  again  stiffens. 

This  latter  point  is  approximately  half  the  depth. 

The  test  program  demonstrated  that  representative  aircraft  structure  can 
be  fabricated  and  tested  to  either  obtain  or  verify  load-deflection  char¬ 
acteristics  within  a  reasonable  budget,  particularly  if  static  tests  are 
performed.  Hie  average  specimen  cost  in  this  program,  including  dynamic 
testing,  specimen  fabrication,  instrumentation,  test  support  and  data  re¬ 
duction,  is  between  $2500  and  $3000. 

For  each  specimen,  photographs  of  the  post-test  condition  are  included  in 
the  Test  Program  section.  In  addition,  a  motion-picture  film  combining 
footage  from  each  of  the  tests  and  depicting  the  test  setup,  crushing  of 
the  structure,  and  the  post-test  condition  of  selected  specimens  has  been 
prepared  as  part  of  the  study  described  herein. 
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CORRELATION  BETWEEN  OEST  AND  ANALYSIS 


A  comparison  of  the  predicted  crushing  characteristics  versus  test  results 
is  presented  for  all  12  specimens.  Reasonable  agreement  between  test  and 
analysis  is  obtained.  For  the  specimens  supported  on  four  edges  (Speci¬ 
mens  1  through  11),  the  comparison  of  test  and  analysis  using  load  cell  date 
is: 

•  Peak  failure  load  (percent  error  between  -  10.6$  and  +  8.0$) 

•  Energy  absorption  (percent  error  between  -  16.1 %  and  -  1.5%) 

For  the  dynamic  tests  (specimens  U,  5  and  6)  using  only  accelerometer  data 
the  difference  between  analysis  and  test  for  the  peak  load  varied  from 
-  2k%  to  +  6.8$  and  for  energy  absorption  from  -  1$  to  +  26$.  For  the 
line-supported  specimen  (No. 12),  the  comparison  of  analysis  and  test 
showed  a  difference  of  approximately  +1.3$  and  -  27$  for  peak  failure 
load  and  energy  absorption,  respectively. 

The  results  of  the  correlation  indicate  than  an  analysis  to  determine  load- 
deflection  characteristics  of  multiweb  constructed  angle  stiffened  panels, 
typical  of  the  UH-lh  lower  luseiage,  is  feasible  using  relatively  simple 
techniques.  The  overall  crushing  behavior  of  the  substructure  can  be 
predicted  within  reasonable  accuracy  although  the  failure  mode  of  each 
individual  stiffener  cannot  be  consistently  accounted  for.  The  analyt¬ 
ical  procedure  described  herein  is  capable  of  satisfactorily  determining 
the  two  most  significant  parameters  for  crushable  structure,  i.e.,  peak 
failure  load  and  energy  absorption  capability.  The  approach,  as  devel¬ 
oped  in  this  study,  has  limitations  with  regard  to  accurately  predicting 
the  deflection  at  which  failure  occurs  and  the  deflection  at  which 
post-failure  stiffening  occurs.  Although  the  predicted  failure  deflection 
is  a  consistently  lower  value  than  that  obtained  from  tests,  it  represents 
an  extremely  small  percentage  of  the  total  crushing  energy  and  this  dis¬ 
crepancy  is  not  significant  as  long  as  the  peak  failure  load  and 
load-deflection  characteristics  are  reliably  represented.  The  analytical 
assumptions  regarding  the  mode  of  failure  in  the  post-failure  region  and 
the  hinge  formulation  at  failure  are  generally  considered  valid  in  light 
of  the  test  results.  For  example,  although  failures  in  the  form  of  an 
asymmetrical  mode  shape  as  well  as  symmetric  mode  failures  (as  initially 
assumed)  occur,  the  predicted  energy  absorption  capability  of  the  struc¬ 
ture  is,on  the  average,  still  very  close  (  =  +_  lU$)  to  the  test  results. 
Assuming  the  asymmetric  hinge  formation  sensitivity  factor,  shown  in 
Figure  95,  was  applied  to  all  specimens,  the  energy  absorption  would 
be  raised  approximately  15$.  While  improving  some  comparisons  with  test 
results,  it  would  create  a  larger  disparity  for  other  comparisons. 

Although  the  failure  modes  experienced  during  the  tests  show  both  types 
(symmetric  and  asymmetric)  of  failures,  the  initial  assumption  of  a 
symmetric  failure  is  "on  an  average"  reasonable. 
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The  initial  assumption  of  a  conservative  plastic  hinge  moment  formation 
(fully  warped  angle  stiffener  flange)  results  in  a  lower  predicted  fail¬ 
ure  deflection  than  is  observed  during  the  test  portion  of  the  program. 
Therefore,  a  modified  post-failure  load-deflection  curve  is  formulated. 
Although  this  improves  the  comparison  between  the  test  and  analysis 
obtained  failure  deflections,  the  analytical  values  still  tend  to  be 
consistently  lower  than  the  corresponding  test-  value. 

For  the  prediction  of  the  pre-failure  stiffness  (pre-  and  post-buckling) 
in  the  analysis,  it  is  assumed  that  the  yield  strain  is  reached  at  fail¬ 
ure.  Due  to  uncertainties  regarding  the  initial  structural  Imperfections 
and  their  effect  on  the  post-buckled  stiffness,  the  overall  pre-failure 
stiffness  is  difficult  to  predict.  The  approach  described  herein,  although 
approximate  and  resulting  in  a  slightly  higner  elastic  pre-failure  stiff¬ 
ness,  is  considered  adequate  because:  (l)  the  discrepancy  i&  minor  in  re¬ 
lation  to  the  total  crushing  energy  and,  (2)  the  analysis  represents  the 
proper  failure  load. 

Figure  127  shows  a  summary  of  the  comparison  between  test  and  analysis 
results  for  load  and  energy  absorbed.  A  detailed  discussion  of  the 
correlation  between  analysis  and  test  is  presented  in  the  Correlation 
section. 

In  particular  the  validity  of  the  following  analytical  assumptions  is 
discussed  in  detail: 

•  Plastic  hinge  formation 

•  Empirical  corrections  to  the  post-failure  load-deflection  curve 

•  Determination  of  the  deflection  value  at  failure 

•  Strain  data  and  failure  modes 


Program  "KRASH"  is  refined  to  facilitate  its  usage  by  designers.  The  input 
data  format  is  changed  for  ease  of  data  input  and  subsequent  data  changes, 
such  as  is  normally  required  for  parameter  tradeoff  studies.  The  input 
format  changes  are  divided  into  three  categories: 

•  reordered  data 

•  standardization  of  certain  inputs 

•  allowance  for  more  general  KR  curves 

As  a  result  of  the  format  and  associated  programming  changes,  the  program 
capability  has  been  expanded  to  at  least  80  masses  and  100  internal  beam 
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elements  and  120  Stiffness  Reduction  Curves  simultaneously.  The  program 
capacity  for  lumped  masses  and  internal  beam  elements  was  previously  40 
and  80,  respectively.  However,  program  run  time  is  still  the  same  for 
problems  of  equal  size. 

The  refined  program  KRASH  has  been  exercised  to  show  that  results  are  in 
agreement  with  previously  obtained  results  using  the  original  version  of 
the  program  (Reference  l).  Case  31-52,  which  corresponds  to  the  test  con¬ 
dition  described  in  Reference  1,  was  used  as  the  verification  case.  An 
updated  program  source  deck*  listing, sample  input  and  output  are  being 
submitted  to  the  Fort  Eustis  Directorate. 


Program  capability  for  special  impact  conditions,  not  previously  performed, 
is  demonstrated  using  the  refined  program  for  the  following: 

•  Three-directional  velocity  impact  (  ~  1+0  ft/sec  combined  impact  - 
27  ft/sec  longitudinal,  23  ft/sec  vertical  and  18.5  ft/sec  lateral) 

•  ’  Upper  mass  penetration  into  a  specified  occupiable  volume 

•  Simplified  blade  contact 

•  Utilization  of  load-deflection  data  obtained  from  specimen  testing 
during  this  program 

Energy  balance  equations  are  included  in  Volume  II  along  with  a  sample 
energy  balance  data  format.  Volume  II  also  contains  a  revised  User's 
Guide  and  associated  format  revisions. 


DESIGN  PROCEDURES 


Three  means  by  which  crashworthiness  capability  can  be  specified  are  pre¬ 
sented,  including: 

•  Buyer-specified  velocity  envelopes  and/or  peak  acclerations 

•  Manufacturer-specified  capability  including  procedures  by  which 
compliance  is  shown 

•  Specified  procedures  which  must  be  adhered  to 

A  discussion  is  presented  in  which  the  general  approach  for  performing 
crashworthiness  analysis,  particularly  during  preliminary  design,  is  de¬ 
scribed.  Included  in  the  general  approach  is: 

•  Identification  of  overall  failure  mode 

•  Identification  of  individual  element  failure  modes 
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•  Prediction  of  load-deflection  curve  for  each  element 

•  Prediction  of  substructure  total  load-deflection  curve 

•  Simplified  representation  of  load-deflection  curves  for  typical 
structural  elements 

A  detailed  design  procedure  is  included  by  which  load-deflection  charac¬ 
teristics  of  multiweb  composite  substructure,  typical  of  fuselage  design, 
can  be  formulated  by  designers.  The  procedure  includes  tables,  figures, 
nomographs  and  equations  and  is  presented  in  a  step-by-step  sequence.  In¬ 
cluded  in  the  sequence  is: 

•  Failure  load  analysis 

•  Beam  bending  analysis 

•  Composite  load-deflection  curve  formulation 

In  conjunction  with  the  detailed  design  procedure  a  comprehensive  sample 
problem  of  a  typical  segment  of  the  UH-1H  lower  fuselage  is  presented  which 
illustrates  the  manner  in  which  design  procedure  is  applied. 

Included  in  the  design  procedures  section  is  a  subsection  which  presents 
structural  element  modeling  guidelines.  Approximate  load-deflection  curves 
are  presented  for  some  typical  Type  I,  II  and  III  structural  elements.  In 
particular  available  pertinent  data  for  the  following  structures  are  pre¬ 
sented  : 

•  Solid  cross-section  column 

•  3Mn  walled  column 

•  Honeycomb  panel 

•  Frangible  tube 


CONCLUSIONS 


1.  Currently  available  methods  for  analyzing  the  nonlinear  behavior  of 
structures  are  oriented  toward  emphasizing  the  analysis  methods  or 
solving  classical  types  of  problems  rscher  than  providing  tools  to  be 
used  in  practical  design  applications. 

2.  A  simplified,  experimentally  verified  analytical  technique  for  pre¬ 
dicting  crushing  characteristics  of  one  type  of  typical  helicopter 
fuselage  airframe  structure  has  been  developed.  The  technique  is 
oriented  toward  designer  usage.  The  technique  is  used  to  predict  the 
peak  failure  load  and  energy  absorption  characteristics  of  the  struc¬ 
ture  with  satisfactory  accuracy  for  preliminary  design. 

3*  Tests  of  representative  helicopter  airframe  structure  have  demonstrated 
that  static  tests  provide  satisfactory  load-deflection  data  suitable 
for  preliminary  design  such  as  peak  failure  load,  failure  deflection, 
energy  absorption  and  curve  shape  when  compared  to  similar  data  obtained 
from  dynamic  tests.  However,  dynamic  tests  pr oxide  immediate  informa¬ 
tion  regarding  springback  effects,  which  in  a  static  test  are  not  evi¬ 
dent  until  relaxation  takes  place. 

4.  Load  sensitivity  studies  have  shown  that  the  energy  absorption  capa¬ 
bility  of  structural  elements,  as  well  as  their  particular  location  in 
the  vehicle,  influences  the  degree  of  accuracy  with  which  they  should 
be  modeled.  Modeling  of  post-buckling  behavior  for  noncritical  energy¬ 
absorbing  regions  can  be  approximated  with  simplified  load-deflection 
curve  representation. 

5.  Program  KPASH, developed  and  verified  under  a  previous  Army  study,  now 
has  the  capability  for  80  lumped  masses  and  100  internal  beam  elements. 
The  program  has  been  demonstrated  "qualitatively"  to  be  operational  for 
analyses  of  three-directional  velocity  impacts,  mass  penetration  into 
an  occupiable  volume,  and  simplified  blade  contact. 
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SUBSTRUCTURE  ANALYSIS 


a 

A.. 


w 


s 


w 

C 

c 

d 

E 

e’eo 

f 

f 

cr 

f 

w 

G 

g 

H 

I 

Jo 

K 

KE 

K 

w 

L,l 

L’ 
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area;  for  stiffened  panel,  area  of  stiffener  plus  area 
of  sheet  corresponding  to  stiffener  spacing,  in. 

plate  length,  in. 

2 

area  of  stiffener,  in. 

2 

area  of  web,  in. 

width  of  skin,  in. 

width  of  stiffener  web,  in. 

constant  (=  3.0 ) 
number  of  cuts 
rivet  diameter,  in. 
modulus  of  elasticity,  psi 
distance  to  neutral  axis ,  in. 

number  of  flanges 
crush  strength,  lb 

effective  rivet  off-set,  in. 

acceleration  level 

2 

gravity  tern  =  386.2  in. /sec 
distance,  in. 

k 

moment  of  inertia,  in. 

4 

moment  of  inertia  about  c.g.,  in. 

eccentricity  factor 
kinetic  energy,  in. -lb 
wrinkling  coefficient 

one-half  column  length,  in. 
effective  length  =  L/  \fc,  in. 
effective  moment,  in. -lb 


plastic  hinge  moment,  in. -lb 
neutral  axis 

effective  axial  force,  lb 
plastic  hinge  axial  force,  lb. 

rivet  pitch,  in. 
applied  axial  force,  lb 
monolithic  failure  load,  lb 

wrinkling  failure  load,  lb 
distributed  load,  lb/in. 

stopping  distance,  core  radius  of  given  cross  section,  in. 

distance  from  'midcross-section  of  skin  to  neutral  axis 
of  stiffener,  in. 

ratio  of  o/a 
1  cy 

ratio  of  critical  stress  (cr  )  to  yield  stress  (cr  ) 

cr  cy 

minimum  core  thickness,  in. 

skin  thickness,  in. 

stiffener  thickness,  in. 

impact  velocity,  in. /sec 
weight,  lb 

effective  half  width  of  skin,  in. 
lateral  deflection 

function  of  column  length  and  end  shortening,  in. 

defined  in  equation  51 

plasticity  reduction  factor 

cladding  reduction  factor 

Pi  =  3.14 

tan"1  Z 

radius  of  gyration  =  v/l/A,  in, 
partial  derivative 


SUBSTRUCTURE  ANALYSIS  (Contirued) 


6u 

6,  A,  u 
v 

X 

a 


co 


cr 


cy 


cy 


cy 


w 


°y« 


fr 


w 


20 

Subscripts: 


b 

f 

i 


incremental  displacement,  in. 
displacement,  end  shortening,  in. 

Poisson’s  ratio 
distance  to  c.g.,  in. 
compressive  stress,  psi 

stress  at  outermost  fiber  of  sheet  skin,  psi 

column  buckling  stress,  psi 

buckling  stress,  psi 

compressive  yield  stress,  psi 

effective  compressive  yield  stress,  psi 

effective  compressive  yield  stress  for  stiffener  web,  psi 

effective  compressive  yield  stress  for  skin,  psi 

Euler  stress,  psi 

crippling  (or  failure)  stress,  psi 

failure  stress  of  riveted  panel,  psi 

failure  stress  due  to  wrinkling  mode,  psi 

buckling  stress  at  L'/p  =  20 


bending 

filler 

I*111  segment 


n&.*. 

min. 

P 

s 

w 


maximum 
minimum 
plastic  state 
skin  or  sheet 
stiffener  web 
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SUBSTRUCTURE  ANALYSIS  (Continued) 

'Constanta  (function  of  stiffened  panel  configuration) 


g 

m 

n 


PROGRAM  "KRASH" 

N 

rotational  transformation  matrix  from  body  axes  to 
ground  axes 

k 

T 

transpose  of  matrix 

CE 

total  crash  spring  ( external  spring)  energy  absorbed 

DE 

total  damping  energy  dissipated 

DRI 

dynamic  response  index 

^TOT 

total  system  energy 

^ijtf 

till 

internal  beam  damping  force  (or  moment)  for  the  ij 
beam  in  the  fob  direction 

FSPijk 

Vi 

crash  spring  forces;  spring  ij  in  the  k  direction 

Fijj? 

force  (or  moment)  at  point  j  due  to  beam  ij,  in  the 
£th  direction 

T  T 

’  Vi*  Zi 

moments  of  inertia  of  lumped  mass  m^  about  i  body 
fixed  axe 8 

Ixyi>  Iyzi*  Pzxi 

*t/h 

product  of  inertia  of  lumped  mass  about  i  body 

axes 

KE 

total  kinetic  energy 

■^ik 

length  of  vector  from  m^  to  ground  contact  point  C"  ^ 

mi 

4th  1 

i  lumped  mass 

N 

number  of  lumped  masses 

PE 

total  potential  energy 

SE 

total  strain  energy  absorbed 

TERM^j  term  in  expression  for  crash  spring  energy 
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f1  (Continued) 


xi*  yi*  Zi 


i*  i*  i 


P4  >  <1«  >  r4 


Xci»  Yci*  Zci 
Lci’  Mci’  Nci 


Axi  » Ayi»  Azi 


weight  of  i  lumped  mass 
ground  coordinates  of  m^ 

i^h  body  axes  component  of  absolute  translational 
velocity  vector  of  mass  i 

Xtk 

i  body  axes  components  of  absolute  angular 
velocity  vector  of  mass  i 

crash  (external)  forces  and  moments,  ith  body  axes 
=  [Avc,  , 

(  1  th 

incremental  displacement  of  point  i,  i  body  axes 


Ainp^,  Ainq^,  Ainri  incremental  rotation  of  point  i,  ita  body  axes 

Avb. incremental  displacement  vector  of  point  j  with 

respect  to  point  i,  due  to  deformation  of  beam  ij, 
in  the  direction 


USER'S  GUIDE  AND  REFERENCE  1  TERMS 


Terms  used  In  Euler* 


s  equations  of  motion 


M 

iq 

W 

MJ 


Rotation  transformation  matrix  from  body  exes  to 
ground  axes 

Matrix  relating  (fa,  to  (p,,  q1#  rj  in 

equation  (92) 

Time  derivative  of  W 

Rotation  transformation  matrix  from  ith  body  axes 

to  c.g.  axes 


Rotation  transformation  matrix  from  beam  ij  axes  to 
body  1  axes 


M  Rotation  transformation  matrix  from  c.g.  axes  to 

ground  axes 

M  Matrix  relating  (0*,  ■»*,  0^*}  to  (p*,  q*,  r*)  in 

equation  (88) 


avetjk 

d^cijk 

M 

O’] 

^ijkl 

^Ljkl 

^pijk 


End  point  of  kth  spring  on  ith  macB 

Ground  contact  po^t  of  kth  spring  on  ith  mass 

Ground  axes  components  of  rector  from  mi  to 

Ground  axes  components  of  vector  for  to  Clk 

Derivative  matrix 

Derivative  matrix 

Running  time  sum  of  AfM^j^ 

Value  of  FM^j^  time  of  loading  reversal 

Body  1  axes  components  of  spring  force  at  ground 
contact  point 


FSPOik  Axial  compressive  force  in  kth  spring  on  lJ;h  mass 

Value  of  FSPOjj^  at  time  of  loading  reversal 
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'•'m» 


,&> 

f-t. 
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USER'S  GUIDE  AM)  REFERENCE  1  TERMS  (Continued) 


VC*iJk 

•vc' 

•p*' 

V* 


Wi 

wT0T 

xtij>  y^iji 

■  Aij 

XG# 

yg» 

*G 

(»«j) 

XG> 

8g 

If 

XG# 

it 

yw 

*G 

xi> 

n* 

*i 

t 

xi. 

YU 

i 

*i 

(viPij) 

v" 

xi> 

Y'i, 

-It 

zi 

(▼iPPij) 

xiJ>  yij*  ziJ 


»  I  I 

*ij>  *u 


/x3i>  Yji*  zjA 

\Lji#  Mji>  Nji/ 


Ground  axes  components f  of  absolute  velocity  of 
ground  contact  point 

Velocity  vector  of  c|k  with  respect  to  m* 

Velocity  vector  of  Cjk  with  respect  to  ground 
Velocity  vector  of  mi  with  respect  to  ground 
Weight  of  ith  lumped  mass 
Total  vehicle  weight 
Beam  ij  coordinates 

Ground  coordinates  of  initial  (t  »  o)  c.g.  position 

Ground  axes  components  of  initial  (t  •  o)  c.g. 

▼air*  tty  vector 

Helicopter  axes  coordinates  of  vehicle  c.g. 

(point  G) 

Ground  coordinates  of  % 

Coordinates  of  »i  in  center-of -gravity  coordinate 

system 

Coordinates  of  mi  in  helicopter  coordinate  system 

Ground  coordinates  oi"  vector  from  point  i  to 
point  j 

ith  body  coordinates  of  vector  from  point  i  to 
Point  J 

Total  (summed  over  time)  internal  forces  and  moments 
at  poi nt  i  due  to  beam  i  j ,  i^b  tody  axes 

Total  (summed  over  time)  internal  forces  and  moments 
at  point  J  due  to  beam  ij,  Jth  body  axes 


2:3 


USER’S  GUIDE  AND  REFERENCE  1  TERMS  (Continued) 


J 


Yi*  Zi\ 
\h*  *v  Ni/ 

XM  YM>  *Ai 

/xci*  Yctt#  zcA 
V  Lci>  M&*  NcJ 


X^,  Y^,  Gravity  forces,  1°“  body  axes 

xIi>  YIi»  2jA  Internal  forces  and  moments,  ith  body  axes 

LH*  Nii/ 

^^ijk  Ground  axes  components  of  spiring  force  at  ground 

contact  point  C^,  positive  up,  left  and  aft 

7jCMAX  Vertical  distance  from  c.g.  to  lowest  cJk 

Determinate  expression  used  in  equation  (68) 

^Fijk  Incremental  forces  and  moments  at  point  ,5  due  to 

beam  ij 

A^iikl  incremental  load  due  to  1th  incremental  deflec¬ 

tion  for  beam  lj 

Ajtfbij ,  Atb^j,  Ap>j«  Incremental  rotations  of  point  J  with  respect  to 

point  i,  in  beam  ij  axes 

Incremental  change  in  ith  mass  Euler  angles 

At  Numerical  integration  time  interval 

Avb . ,  Six  element  vector  made  up  of  Axb, ,,  Ayb. ,,  Azb., .. 

1J  A(Sb1J(  Afcij,  £f*1}  1J  1J  1J 

A^ij  Incremental  displacement  vector  of  point  J  with 

respect  to  point  i,  due  to  deformation  of  beam  ij 

Incremental  displacement  vector  of  point  j  with 
respect  to  point  i 


/xIi>  Yn,  ZjA 
\LIL>  NIi' 


Total  forces  and  moments  on  mass  1,  in  ith  body  axes 


Aerodynamic  forces,  i  “  body  axes 

Crash  (external)  forces  and  moments,,  1th  body  axes 


Ai 

AFijk 


A^j 

—  bm1 


i 
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USER’S  GUIDE  AND  REFERENCE  1  TERMS  (Continued) 

Final  value  of  FSPO^  in  input  table  of  s^  vs.  FSPOik 
Center-of-gravity  of  total  vehicle 


G 

H 


Origin  of  helicopter  coordinate  system  (F.S.O, 
B.I..Q,  W.L.O) 


Kexi,  Heyj,  Hesi 

Jxi*  ^i »  hi 
Ixyi t  *yzi# 


k*ik 

hj 

OuJ 

hx. 

^xi>  ^yi>  ^zl^ik) 

^ci 

UFT1 

“1 

®“lk 


N 

“ik 

V 

0 


Angular  momenta  of  ml  due  to  rotation  of  masses  in¬ 
ternal  to  m^ 

Moments  of  inertia  of  lumped  mass  mi,  about  ith 
body  fixed  axes 

Products  of  inertia  of  lumped  mass  m^,  about  ith 
body  fixed  axes 

Linear  unloading  stiffness  for  kth  spring 

Six  by  six  linear  stiffness  matrix  for  beam  ij 

Six  by  six  diagonal  stiffness  reduction  matrix  for 
beam  lj 

Length  of  vector  from  in^  to  ground  contact  point  c{k 
Free  length  of  kth  spring  on  ith  mass 
Aerodynamic  lift  constant 

Aerodynamic  lift  on  ra^,  positive  up,  in  ground  axes 
ith  lumped  mass 

Ground-spring  friction  coefficient  for  kth  spring 
on  ith  mass 

Total  number  of  lumped  masses 

Unit  vector  triad  fixed  in  ith  body  coordinate  system 
Unit  vector  triad  fixed  in  ground  coordinate  system 
Origin  of  ground  coordinate  system 


■) 

i 
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.USER'S  GUIDE  AND  REFERENCE  1  TERMS  (Continued) 


Pi>  ri 


PS  r* 


*1X 

»ik 

•?ik 

#lk 

®ik 

t 

ftj 

«i>  vi*  wi 

™ij 


vb 


1J 


*>1,}1 

▼ik 

^ijl 


vb 


IJ1 


VC 


ijk 


ith  body  axes  components  of  absolute  angular  velo¬ 
city  vector  of  aass  1 

c.g.  axes  components  of  initial(t«0)  vehicle  angular 
velocity  vector 

Contact  point  velocity  matrix  used  in  equation  (60) 

Axial  external  spring  compression,  kth  spring  on 
1th  mass 

Value  of  s^  at  time  of  loading  reversal 

final  value  of  s^  in  input  table  of  •lk  vs.  ISPO^ 

kth  spring  axial  compression  measured  relative  to 
current  load  stroke  curve  origin 

Haris sntal  shift  of  *&  coordinates  with  respect  to 
Six  coordinates 

Time 

Static  balance  matrix  used  in  equation  (30b) 

Body  1  axes  components  of  absolute  translational 
velocity  vector  of  point  mi 

*i*  *±»  *i 

Running  time  sum  of  Avb^ 

Value  of  vbijj.  at  time  of  loading  reversal 

Magnitude  of  ground  plane  contact  point  velocity 

1th  total  beam  deflection  measured  relative  to 
current  load-stroke  curve  origin 

Horizontal  shift  of  Vb^  coordinates  with  re¬ 
spect  to  vbijj  coordinates 

Ground  coordinates  of  point 
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■1 

V* 

.ft 


% 


"«S 

■*s 


| 


I 

| 
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USER'S  GUIDE  AND  REFERENCE  1  TERMS  (Continued) 


Incremental  displacement  vector  of  point  J  with 
respect  to  point  1,  due  to  rotation  of  mass  i 

Coordinates  of  in  beam  ij  axes 

&*i>  6yi,^zi 

Incremental  displacement  of  point  i ,  ground  axes 

Incremental  displacement  of  point  J  with  respect 
to  point  i  in  ground  axes 

/toy,  £11^,  &  %j\ 

V  6Liy  &:,ij 

Incremental  internal  forces  and  moments  at  point 
d  due  to  beam  ij,  in  beam  ij  axes  (elements  of 
AFij  vector) 

\  AI1J»  / 

Incremental  internal  forces  and  moments  at  point  J 
due  to  beam  ij,  ground  axes 

/  Axfj>  \ 

p?j  -p?/ 

Incremental  internal  forces  and  moments  at  point  i 
due  to  beam  ij  ,  ground  axes 

/n*Jj»4lty,  AZy\ 

( ^Lij» 

Incremental  internal  forces  and  moments  at  point  i 
due  to  beam  i  j ,  ith  body  axes 

/A*^i,  Avjj,  AzjA 

Incremental  internal  forces  and  moments  at  point  J 
due  to  beam  ij ,  j*“  body  axes 

$i* 

Euler  angles  from  ground  axes  to  body  axes  (time 
varying) 

Euler  angles  from  ith  body  axes  to  beam  ij  axes 
(constant) 

0',  ' 

Euler  angl.es  from  ground  axes  to  c.g.  axes  ( c di¬ 
stant );  initial  (t*o)  attitude  of  vehicle 

0;, fi 

Euler  angles  from  c.g.  axes  to  ith  body  axes  (con¬ 
stant  ) 
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